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Summary. — A method is established to solve Einstein’s equations by 
successive approximations, uniqueness and compatibility at each step 
being secured by imposing de Donder’s condition on the field. This turns 
out to restrict the freedom of choice of the sources in a way that deter- 
mines their laws of motion, directly from the field quantities, without 
any use of integration around singularities. If the approximation para- 
meter is the reciprocal velocity of light, the first non trivial approximation 
results in the constancy of the rest mass, and in Newton’s law of attrac- 
tion. The next step shows that there is no first-order correction to these 
laws. The third approximation then gives a small correction to the rest 
mass, which is found to be increased by the presence of other bodies. 
There is also a second-order correction to the acceleration (not explicitly 
evaluated here) which causes the perihelion advance in the two-body 
problem. Te fourth approximation involves a contribution from the 
radiation field. It is found that in order to fulfil de Donder’s condition, 
we must introduce into the field a term which is linear in the sources, 
although it does not appear in the so-called « linear approximation ». 
This result may cast some doubts on the validity of the latter. In a 
further paper, we shall solve still higher approximations by means of 
an improved technique, and find that the fifth-order correction to the 
acceleration involves a non-conservative term. This last result may be 
taken as evidence for the reality of gravitational radiation. 


40 - Jl Nuovo Cimento. 
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1. Introduction. 


Exact solutions of Einstein’s equations are known only in an exceedingly 
small number of cases. In most physical applications, the problem must by 
tackled by successive approximations, in which we assume that the solution 
is an analytical (1) function of some small parameter (e.g. the reciprocal velo- 
city of light, or the gravitational constant) and thus can be expanded into 
a power series such as 


gl” = que que gl + Po (a) 
0 1 2 


where each term is by one order of magnitude smaller than the previous one. 
It is then found that the equations to be solved are not, in general, com- 
patible. Indeed, let us write Einstein’s equations in the form 


(1) (qe gn = Gr Oe = 249 T tà Ver + doi can Cae gio a dali Ce ; 


where 


Qe = gg [9 (— 4020.05, + 10° 0.0, FIT Ino) 0,000 — 
m0:0 000 210,0 000 00,0 80 Cea) Cea Cat ile. o) 


and 9,5 — 9 9 

When we compute a”, the right-hand side of (1) depends only g‘, m<n, 
and thus can be considered as a known quantity. But the ten equations for 
g” are not independent: indeed 


n 


[(g”” 9 de. gue sa =) 


and therefore, equations (1) are compatible only if the four-divergence of their 
right-hand side is null. This in turn imposes four additional conditions on 
the previously found approximate solution, which can be fulfilled only for 
special values of ZT". In such a way, compatibility of Einstein’s field equa- 
tions determines the laws of motion. 


(1) B. Berrorrr: Nuovo Cimento, 4, 898 (1956); 6, 755 (1957). 

(*) Greek indices run from 0 to 3; Latin ones, from 1 to 3. It is convenient (2) to 
use as basic functions the contravariant density 9g = gg”, where g = V — (Det 94) 
Pseudoquantities (i.e., those of which the transformation laws involve some power of 
the Jacobian) are usually denoted by gothic letters. 

(?) P. T. Hoane: Compt. Rend., 246, 61 (1958). 
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2. — General covariance. 


When the laws of motion are satisfied, we remain with only six independent 
equations for the ten gu: this corresponds to the possibility of giving to four 
of them arbitrary values, by an infinitesimal co-ordinate transformation of 
the n-th order. In order to secure uniqueness of the solution and to get rid 
of apparent fields, we must impose on gf four additional conditions. These 
can be made fully covariant by introducing, side by side with g“”, an auxiliary 
flat (but not necessarily Minkowskian) metric y (8) which we shall identify 
with the «zeroth» approsimation gi. It is then convenient to choose, as 
additional conditions, de Donder’s relations 


RIP RZ 
9 yO, 


where the comma now denotes a covariant derivative with respect to the affi- 
nity belonging to y"”. It is easily shown, by generalizing a theorem of Fook (*), 
that de Donder’s relation prevents any infinitesimal co-ordinate transformation, 
modifying only g””, (n > 0), besides those generated by harmonic functions (5), 
and thus, under suitable conditions at infinity (°), secures the uniqueness of 
the solution, once y has been chosen (*). 

We shall assume, in this paper, that the velocities are small quantities of 
the first order. (This is equivalent to the choice of e-1 as the expansion para- 
meter. However, it is more convenient to take c=-G= 1). Accelerations, 
and masses are therefore of the second order. It is easily seen that the only 
explicit dependence of the field quantities on time is through the positions 
and velocities of the sources. It follows that a time derivative of. a field quan- 
tity is, in general, by one order of magnitude smaller than a space derivative. 
This, however, is due only to the choice of the expansion parameter (which 
is Lorentz invariant) and does not impair at all the general covariance of the theory 
as stated by some authors (!*). 


(3) N. Rosen: Phys. Rev., 57, 147 (1940). 1 
(4) V. A. Fock: Theory of space, time and gravitation (in Russian, Moscow, 1955), 
. 441. 

È (5) L. Lanpav and E. Lirsuirz: The Classical Theory of Fields (English translation, 
Cambridge, Mass., 1951), p. 324. 

(6) A. TRAUTMAN: Bull. Acad. Pol. Sc., 6, 407 (1958). 

(*) Without the introduction of y”, de Donder’s condition is not covariant, thus 
seriously impairing the generality of the theory (7). 

(7) L. InrELD: Rev. Mod. Phys., 29, 398 (1957). 

(8) P. Havas: Phys. Rev., 108, 1351 (1957). 
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3. — Approximation procedure. 


We shall assume that matter is present in the form of N particles of rest 


A A A 
mass M, and co-ordinates é*= é*(t), A=1,..., N. We have (5-1): 


(2) To 


ly, 
A=N A dé“ dé do | A 
Sant SS ole), 


O zi 
ey dc do 
where 


Ao > Ap 
do* = y, dé" dé". 


The Lorentz-invariance of the rest mass follows from the fact that T“’dtdV 


is a genuine tensor. 
It is easily seen that the first non-trivial terms of g’” are g°°, g* and g”. 
2 3 4 


We shall therefore call n-th approximation the computation of 9°, g° and g”. 
ne nee n+3 


It will turn out that 
(3) CLG + ae = 0 à 
m+ n+2 


will give the (n —1)-th order corrections to the masses, and 


ko kl 
(4) gue g 1 0, 
n +2 n+3 
the (n —1)-th order corrections to the accelerations. 
As we are dealing with a problem having no special symmetry, it is conve- 
nient to chose y” = 7”, where 7” is the Minkowskian metric: 7° — 1, 7°*= 0, 


yi =—o™, With the help of de Donder’s conditions, Einstein’s equations 
now read 

(5) VA ge a 1679T" ate Ee” ; 

where 


Qu — qu a (ge 7 D) A 


The ten equations (5) are independent, and therefore always compatible. 
If we neglect 0’, we get the so-called «linear approximation » which is the 


basis of usual computations on gravitational radiation (11). It is therefore 


(°) J. LUBANSKI: Acta Phys. Pol., 6, 356 (1937). 
(0) W. TuLezyEw: Bull. Acad. Pol. Sc., Ol. IIT, 5, 279 (1957). 
(4) L. Lanpav and E. LirsHITZ: loc. cit., p. 330. 
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interesting to examine whether Einstein’s equations (5) have exact solutions, 
in which the part of gf” which is linear in the masses is identical with the «linear 
approximation ». 

We shall therefore write 


ge = CE pp 3 2 uv 
where 
V7"? Sd he CATA 1679T/ 
is the linear approximation, and 
V2847 — ELLE Ou ; 


The equation for 6“ can be solved by the use of Lienard-Wiechert potentials: 


M,9 dés dé”. 


TEA, ab ; ie = 
) > \R,(d&/do)| do do dai 
where R? = 2° &, and we have omitted the superscripts A. 
Defining the velocity 
dé’ 
re : o=1, 
and the « effective source » 
dt 
M = M, do g , 
this can be written as 
Moto È ate, 
p= 4 DE — dn computed at time (tF R), 


where R—4RERE and the + sign is to be taken according to the choice of 
retarded or advanced potential. For more generality, we shall take a linear 
combination of these potentials, such as 


1 pie € 
bY” = (D + D) SEI), 
dd 


where e is an arbitrary constant. It may be shown (!2) that 
(co) 

=> == say R2n— zs 

5 | Déav) Tr Gaia > wR 


(12) A. S. EDDINGTON: The Mathematical Theory of Relativity (Cambridge, 1924), 
p. 249. 

(13) F. R. CROWNFIELD jr. and P. Havas: Phys. Rev., 94, 471 (1954). 

(14) B. BERTOTTI: Nuovo Cimento, 2, 231 (1955). 
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Il d2r+1 
+)! dit 


ib 14 2 = v n 
— 3 (Hi — Hie) = 4D oe PEARLS” 


where all quantities are now to be computed at the field time 4. 

In these expressions, we suppose that the co-ordinates &* and velocities v* 
are known (initial conditions), and that the masses and accelerations have 
to be expanded 


M=m+Mm+M+.., 
vw at ata ak... (*). 


4. — First approximation. 


We have 


m 
lo 00 2= h® = 4 — 
ee eae: 
mot 
qu pr A TS —_, 
MU Der 
Application of (3) gives 
m 
a= 0, 


which can be satisfied if and only if m-=0. The rest mass is therefore a 
constant at this approximation. 
Further, we have 


mory! 
pti 4 : 
4 2 Ro 


but in order to get 8%, we first have to compute 0”: 
4 4 


00 T 
— — 00 00 
9 Pal 9 n 9 Hone) 


a oy 


(or = 0, 
4 


Qui 20 00 00 1 00 ( 
cl — —jk ne — ‘0 È 
; 809 nd ng 9 à 


(*) Note that a* and m are small quantities of order (n+2). 
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Fr i 15 
rom this, we deduce (5): 


MOST OU 


AB 
= eee ee (loc 4 DSK San 4 Jkt 
9 > m pi jet pei (108 À) 2 mm (26% S 484), 


where terms with A= B are excluded from Di and 


gr SA Lente (R + R a) I; 
og* 0€! 


D being the distance between particles A and B. 
Application of (4) gives 


at SM pr ; 
B 
whence 
AB 
‘ 1 
I mm (= DI Set, ae = su) =) 
A,B R 4 
The coefficient of inm is 
AB i AB BA 
It k Fr D*|/D3 Fe— De! D3 
= NIE IE ari + 3 l AUS 
Ik R R R 


A B 
where Dt — &*— & , 
Newton’s law of attraction directly follows. 


(CS WW PA SE OCR GY, Cilio, pato 80% 


623 


624 A. PERES 


5. — Second approximation. 


We have 


because of 7 — 0 and Y ma*= 0. (This latter relation follows from the just 
found value of a*). Application of (3) gives > (n/R) = 0 whence it follows 
that # is a constant, which can be introduced into m. We thus may take, 
without loss of generality, #= 0, whence g°—0 and g°*= 0. 


3 4 
The next step necessitates the computation of 0”. We get 
a 


A” — 0 ; 


Aor — go (qe — qon ) = o a. qo 
5 5 sn 3 sn 9 sk 4 3 ,0 9 ie) 
Orr 0. 


5 


Thus 
gS FE de D m(vrat+ vat) 4 
5 5 


This is a function of time only: g*,= 0. It does not follow, however, 
Ho 


that (4) holds identically, because in the course of the first approximation we 


have made 
il ve it ak 
We —— 4 ax (È / mr | — = 
go > m [o Jero (a) +I: 


\ 


It therefore remains, for the second approximation 4 ¥ (ma*/R) == ()s 
whence a* = 0. 


6. — Third approximation. 


We have 
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and by a method similar to Fock’s (5) we get 


m* 48. 
80 — 7 $S— +145 mmS", 
4 ER? AB 


0% 


oF CO, 2 SI bee 
ica sn (los n) +16 Y mm (san + ok SU — fau). 
L 5 ; 


eo 
22 


1 IY 


Application of (3) gives 


™m "7 
2 ARE RS À) Qu Be B 3 
De + D' mm 5 SUL SDS, + 4urg,, — Avis) — 0. 


S 
& 


We tentatively try a solution such as 


na A B . . 
and equate to zero the coefficient of mm. This gives 


AB A È BA E “ 
HÉROS NON) A) 
I B Pa O 


A 
R Vir 


This result, together with Newton’s law of attraction which has just been 


derived, leads to 


This is the second-order correction to the «effective source » M = M,(dt/do)g. 


The correction to the rest mass M, is then found to be (*) 


B 
A A , My 
Mo = Mi. 
2 B 


S 


(*) At this approximation g=1+49%=1+42) (m/R). As this is computed at the 
A 2 
point R= 0, we have to omit the infinity self-energy of the particle A. This can be 


done in a consistent way by the use of « good » 6 functions (1). 
(6) L. InreLp and J. PLEBANSKI: Bull. Acad. Pol. Sc., Cl. III. 4, 689 (1956); 


5, 51 (1957). 
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The rest mass is thus increased by the presence of other bodies, a result 
that has already been found, on quite other grounds, by EINSTEIN (7). 

The next step is to compute OF in order to get ghi We find that OO, 
but O° and O*! are extremely cumbersome expressions We shall therefore 
postpone the evaluation of @* to the following paper, in which we develop a more 
powerful method, enabling us to obtain the laws of motion with a much smaller 
amount of computational labour. This will be, however, at the expense of 
mathematical elegance and logical simplicity. 


7. — Fourth approximation. 


As O0— @*=— 0, we should have 
5 6 


00 00 a € pi 
3 =!) = 4Y pt em + ay mh) ’ 


mo . € 506 
qo" Ss por ae fl Di È ne E(MU) + 31 m(v*R?2) + ome o 


However, application of (3) leads to 


which implies both m= const — 0, and Y im = 0. This last relation is not, in 
general, satisfied. We thus have to add to h* a term such as 
6 


go =—¥e)> (mR*) . 
(Note that Visa 0.) That is, we have to introduce a term which is linear 
in the sources, although it does not belong to the «linear approximation » 
This result may cast some doubt on the validity of the latter. 

As most computations on gravitational radiation by a system of bodies 
are based on this linear approximation (11), it seems to us necessary to pursue 


(7) A. Einstein: The Meaning of Relativity (Princeton, 1953), p. 102. 
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the approximation procedure until we find the correction to the acceleration 
due to radiation reaction. This will be the subject of a later publication. 


I am indebted to Professor NATHAN ROSEN for stimulating discussions. 


RIASSUNTO (*) 


Si elabora un metodo per risolvere le equazioni di Einstein per approssimazioni 
successive, l’unicità e la compatibilità essendo assicurate ad ogni passaggio imponendo 
al campo la condizione di de Donder. Ciò equivale a restringere la libertà di scelta 
delle sorgenti in modo da determinare le loro leggi di moto direttamente dalle srandezze 
di campo senza far alcun uso di integrazioni intorno a singolarità. Se il parametro 
d’approssimazione è il reciproco della velocità della luce, dalla prima approssimazione 
non banale risulta la costanza della massa a riposo e la legge d’attrazione di 
Newton. Il passo successivo dimostra che non esistono correzioni del prim’ordine di 
queste lesgi. La terza approssimazione dà poi una piccola correzione della massa a 
riposo, che si trova aumentata della presenza di altri corpi. C’è anche una correzione 
del second’ordine dell’accelerazione (qui non calcolata esplicitamente) che provoca la 
precessione del perielio nel problema dei due corpi. La quarta approssimazione com- 
prende un contributo del campo radiativo. Si trova che per soddisfare la condizione 
di de Donder dobbiamo introdurre nell’espressione del campo un termine lineare nelle 
sorgenti per quanto esso non appaia nella cosidetta « approssimazione lineare ». Que- 
st’ultimo risultato può far sorgere qualche dubbio sulla validità di tale approssimazione. 
In un ulteriore lavoro risolveremo approssimazioni superiori per mezzo di una tecnica 
perfezionata e troveremo che la correzione del quint’ordine dell’accelerazione com- 
prende un termine non conservativo. Quest'ultimo risultato può essere ritenuto una 
prova della realtà della radiazione gravitazionale. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 15 Ottobre 1958) 


Summary. — Matrix elements of the potential energy of interaction are 
calculated for the case of two spherically symmetrie, localized and non- 
overlapping wave functions. The result is expressed as a power series, 
and is applicable to a wide class of central force potential energies, and 
non-overlapping wave functions. As an illustration, a discussion is given 
of the diagonal matrix elements of the potential energy corresponding 
to a Van der Waals law of force. 


1. — Introduction. 


In many problems one needs to evaluate potential energy matrix elements, 
ap between two wave functions, p,, p;, localized around points P and Q, i.e., 


(1) VP, Q) = | i Pala — P)pi(x — Q) v(x, 22) px (x — P)p;(x, —Q)d°r, dr, . 


In most cases of interest, the potential energy function, v(x;, x»), depends only 
on the scalar |x,— x,|, 7.e., depends only on the distance between the two 


() This research was supported by the U.S.A.F. through the Air Force Office of 
Scientific Research of the Air Research and Development Command, under Contract 
no. AF 18(603)-108. 

(”) On leave from Instituto Tecnologico de Aeronautica, Sio José dos Campos, 
S. Paulo. 


628 


POTENTIAL ENERGY MATRIX ELEMENTS BETWEEN NON-OVERLAPPING ETC. 629 


points, and the gs are spherically symmetric functions of their arguments. 
In this case we may write: 


©) Va(P, Q) =| Jt (rs) fora) 0(0) den drs 
where 


ie 1, =|%.—Q|, 1 == |X, — 2 : 


prs) = P (13) PAT) 9 


and f,(r,) = pylra) @;(12). 

In the following we discuss a general method, suitable for numerical cal- 
culations, to evaluate the six dimensional integrals for V,;(P, Q) for a wide 
class of potential energies of interaction and angle-independent, localized non- 
overlapping wave functions, è.e., to evaluate these integrals with: f,(7,) — 0 
ee Fa, ne =0 when r,>b, and (a+b)<|P—Q|=R. Our results for 

,(P, Q), Eqs. (6a), (6b) below, will be expressed as a power series in terms 
È Da and na and will be found to involve the potential energy function 
v(r) and its derivatives, and the moments of the probability densities f,(7,), 
f,(r,). We shall also discuss the rate of convergence of the series, which depends 
on the functional form of v(r) and on the ratios a/R and b/R. 


2. — Integration. 
The integral yielding the matrix element under discussion is 
(3) PAG, a, 0) =| | flrdflra)e(r) dr dre 
CLIO 


where 2, and ©, are respectively two non-overlapping spheres of radii a and b, 
with centers at points P and Q. 
We may write: 


(4) Pa De fier) PE a. 
(Da 


Vi(e, a) = fieno) dr: 2 


Qy 


the meaning of the symbols 
being given in Fig. 1. 

In order to calculate 
V,(z, a) we first integrate 
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over the rim of the cone: r, — const., 0, — const., and obtain: 


Vila, a) = 2a| ff) sin 0, d0, dr, . 


0 0 


Next we change the variable of integration from 0, to r by means of the re- 
lation: 


r? = 22+ ri — 2er, cos 0, , 


so that 
27 : no 27 : 
Vi (2, a) = n fran fre dr — = [ve a Ty) ae (era US) 
0 2-7, 0 


where by definition: 
w(r) = rv(r) 


and 


W(a) = frein dr =| row) dr, 


t being an arbitrary number (t+) whose value does not have any effect 
on V(z, a). 

Any v(r) of physical interest is such that W(x) is a function which can be 
expanded into a uniformly convergent power series and hence can be integrated 
term by term. These two steps lead to 


eh a DS n ce 
where 
wa = [TT Mi =, 
dy bee 
and 
200 a) = fréftr) dm; Aa) =1, 
(CA 


is the (2n)-th moment of the probability density around point P. 
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Eq. (5) is a series of even powers of a, and represents the potential energy 
produced at a point + exterior to 2, by a particle whose wave function y, is 
effectively confined to the interior of Q,. 

Combining Eqs. (5) and (4) and changing the variable of integration in 
Eq. (3) from 6, to e, yields: 


ome (a 
Vup(R, a,b) = R. > Ont AS, feu + ra) — U(R — 12) | rofa(r2) dre , 
where 


x 


U(x) = fu") d2% 


t 


so that expanding U(x) in a power series, and integrating term by term, we 
finally get: 


1 (co) Aa WAGED ) 
6 Ve, (Ria ae 
(6a) «pl , a, b) DEI (2n +1)! (2 a (R), 
where 
antètn+m)l D d2t+m) 
w AC ae = WY) DEE 
and 


25” 0) = 13" faa) d°r , 


Qi 


is the (2m)-th moment of the probability density around point Q. 
Eq. (6a) represents the desired matrix element. Its meaning becomes more 
transparent if we separate out the first term in the series, é.e., the term with 
—0=n. We obtain 


to ERO) 
(60) Vap(B, a D = BITS 2 Gn+iCm+DI wR) 


n+m+-~0 


The first term, v(R)-1, in Eq. (60) represents the classical potential energy 
of interaction of two point sources located at points P and Q. The remainder 
term, v(R)-Y > ..., is either a polynomial or a power series in ar, b/R, which 


n m 


in either case vanishes in the limit a > 0, b + 0. The coefficients in this term 
can be evaluated in a straightforward way once v(r) and the wave functions 
Palr1); dp(r.) are known. 


_ 
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In the particular case of diagonal matrix elements, «=, and if a=b 
we can rearrange the terms in Eq. (60) in the following way 


(7) Via) (RES MET ES (AR) SI Si 
a (al) ETES here) TE De] 2 uti 


Raf Age a) 22 (a) WEAR) 
Ar 2(s+t) _ E AR A 
(3) SEI) 


3. — Illustrative examples. 


The terms inside the bracket in Eq. (6b) may contain a finite or an infinite 
number of terms, depending on the functional form of v(r): 


a) For a Coulomb inverse square law of force, v(r) — 7-4, and w(r)=const.; 
the remainder term then vanishes identically and V,(R, a, 9) = v(R) in agree- 
ment with the well known result of classical electrostatics. 


b) If v(r) is a polynomial of (p)-th degree in r, p> 0, then our series 
in Eq. (6b) or (7), reduces to a polynomial of the same degree. In this case 
the evaluation of the right hand side of Eq. (6b) or (8) is in general less tedious 
than a direct evaluation of the integral in Eq. (2). 


c) Finally, for all other v(r), the series in Eq. (6b) does not terminate 
but none-the less converges for (a+b) < R. The rate of convergence depends 
on the values of the space derivatives of the potential v(r), evaluated at r= R. 


In order to see the usefulness of our expansion for a v(r) of type c), we, 
in a particular case, compare approximate values of V.,(R,a) with the cor- 
responding exact values. The integral in Eq. (2) for V.,(R, a) has been eva- 
luated in closed form by HADAMAKER (1) for a Van der Waals potential energy 
v(r) = Kr, (K = const). and with probability densities 


f,Ar) = f,(r) = (4maÿ/3) 1, Lorgini sta, 


and 


f,.(m=f;@=0 for r> a. 


(*) H. C. HAMAKER: Physica, 4, 1058 (1937). 
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With these probability densities the moments y©"(a) = 3a2"/(2n+3), so that 
Eqs. (7) and (8) yield: 


(9) V.(R, a) = KRS{1 + (6/4)? 


= 
Le 
Pea 
© 
at 
S 
3 
US) 
CES 
> 
pu 
_ 
ss 
PS 
= 


3) a? + Lato 


where x = 2a/R. On the other hand, the exact result found by HADAMAKER is 


= 3K | 2a? 2a? 4a? 
10 ih) LIn{l ; 
vy sal, €) = 3500 LR 4e 7 RF n | =} 
In Fig. 2 we have plotted, for va- 1 ae 
rious values of 2a/R, the exact value a 
of V (R,a) as given by Eq. (10), R Va 


(curve H), and our result for V.,(R, a) 40 
in Eq. (9) including terms up to 2° 
(curve A), and «8 (curve B). Table I 
illustrates the accuracy of our approx- 
imation. 

From Fig. 2 and Table I we can 
conclude that our series for V,,(R,a) 
converges rapidly for 2a/R <0.5 in 99 
the case of Van der Waals forces. For 
potential energies of the form v(r) — r-? 
(p>1), this series converges more 
rapidly for p< 6, than for p> 6. In 0 0? OZ 06 08 
a paper to appear elsewhere (?), the 
author has applied the present method, 
using Lennard-Jones «12-6» potential 
energy functions and various probability density functions, to study the binding 
energy of crystals at low temperatures. 


30F 


Jr D 


IPABLE TI 
Î 
2a/R AA (%) | AB (%) 
0.2 | 0.01 <0.01 
0.4 0.25 0.20 
0.6 2 1 
0.8 25 16 


(2) N. BERNARDES: to appear in Phys. Rev. (1958). 
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RAS SI UNERONt(S) 


Si calcolano gli elementi di matrice dell’energia potenziale di interazione pel caso 
di due funzioni d’onda a simmetria sferica, localizzate e non sovrapposte. Il risultato 
si esprime in serie di potenze ed è applicabile ad una vasta classe di energie potenziali 
di forze centrali e di funzioni d’onda non sovrapposte. A titolo illustrativo, si dà una 
discussione degli elementi di matrici diagonali dell'energia potenziale corrispondenti 
a una forza di Van der Waals. 


(*) Traduzione a cura della Redazione. 
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The Partition Function 
and the Subsidiary Conditions in the Bohm-Pines 
and Migdal-Galitsky’s Supplementary Variables Method. 


A. PAWLIKOWSKI and W. SZCZUROWNA 


Polska Akademia Nauk, Instytut Fizyki, Pracownia Statystycena - Wroctaw 


(ricevuto il 31 Ottobre 1958) 


Summary. — The supplementary variables method is applied to the 
approximate calculation of the partition function of N interacting fer- 
mions. It is shown that the supplementary variables method in its original 
formulation due to Bohm-Pines and Migdal-Galitsky leads to wrong 
results. One obtains correct results if one assumes a modification, due 
to Kanazawa, of the subsidiary conditions in this method. 


1. — Introduction. 


In this work (*) we calculate approximately the partition function of a 
N-body system using the supplementary (collective) variables method of Bohm- 
Pines (1) and Migdal-Galitsky (?) (hereafter quoted as BPMG). In this method 
the subsidiary conditions which are imposed on the state vectors of the system 
play an important role. Recently KANAZAWA (*) has pointed out that the 
subsidiary conditions of BPMG must be modified analogously to the modification 
of the Lorentz condition in quantum electrodynamics due to GUPTA and 
BLEULER (4). The necessity of this modification results also from our work: 

(*) Some of results obtained here were published as a short note in Dokl. Akad. Nauk 


SSSR, 118, 61 (1958). 
(1) D. Pines: Solid state physics, 1, 368 (1955). 
(2) D. N. ZUBAREV: Candidate's degree dissertation (Moskow, 1953). 
(3) H. Kanazawa: Progr. Theor. Phys., 18, 287 (1957). 
(4) S. N. Gupta: Proc. Phys. Soc., A 63, 681 (1950); K. BLEULER: Helv. Phys. 


Acta, 23, 567 (1950). 
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only the partition function calculated with the use of the modified subsidiary 
conditions fulfils the correspondence principle for # — 0. 

In tbe calculations we use a discrete parameter x which takes two values: 
0 and 1. It enables us to consider simultaneously both formulations of the 
subsidiary conditions: the original and the modified one. 


2. — The supplementary variables method. 


In this Section we follow the papers (7°). 

Let us consider a system of N fermions (the case of bosons can be treated 
quite analogously) contained in a volume v interacting through a two-body 
central potential W,,—W{|x;—x,|). The Hamiltonian of the system is 


(1) HT Ta). 
Here 
N DI 
2 d IE = DI 
(2) or 
Na INS SEE: HER Se N 
(3) De YW.,= D Wont Wo— ge D Wa (*), 
Dei 20% DI DU 
where 


Wi pr x|)exp[—1kx]|dx, 


1 2 
== fs (ODE — ikx, . 
eva > exp [ ] 


dele 


By è; and 9; we shall denote the hermitian and the antihermitian parts of 
the operator V20,, respectively: 


1 
V2 


(4) te (x + 10%) . 


We now extend the above system by adding some new degrees of freedom 
described by variables Q, (the supplementary variables), (k #0 and |k|< ky, 


(‘) B denotes an operator and B his eigenvalue. 
(*) D’ denotes sum over k with omitting the term k = 0. 
k 
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ky is some constant). For these variables we have, similarly as for 6,, 


(5) Q_ 
and 

i A Je _ 
(6) UM O 


We define the extended Hamiltonian as follows 


(7) AL=T+Y7,, (<=, 0), 
where 
a. N Ae N 
(8) Voc >Wodx+ = Wi(0x+ On) (02 + 0-3) + 4e) 
AV |k|>xko AV |k|>%o 
and 
INR N x 
(9) A, = - 0 5 > Wi—— > ho, ; 
2a VE EX 
one 
10 D lle ee ie... 
ae i | mv 


(ab) U = exp 1-3 na Più] ; 


where P, = (1/V2)(P% + PI) is the momentum canonically conjugate to On 


(12) ansi ep 


After performing the transformation (11) we have 


(13) i = 08H 02=7+48,,4+4,4+4,, 


where 


lk] <k 


fee N NES 1 Cas 
(14) oot = ZH ose k Re (sn! P,P _, 4 + Mida) = > (Het Ha) 


is the Hamiltonian of many linear oscillators vibrating with the frequencies @,, 
and È, contains terms describing the interaction between the fermions, between 
the oscillators, and between the fermions and the oscillators (1°). 
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The extended system with the Hamiltonian (7) must be equivalent to the 
original system. Therefore, it is necessary to impose some subsidiary con- 
ditions on the state vectors of the extended system: 


(15) LD = 0 |k|< ko, 
where 
a / NW 
DR A a - k 
(16) Lia Wan at ke ve me Ta 


After the transformation (11) the subsidiary conditions become 


(17) Lgl 0 |E|<ko; 
where 

(18) IO ce 

and 

(19) EUR. 


It is easy to verify that the substitution x=0 in the above expressions yields 
the original formulation of the supplementary variables theory due to BPMG. 
and the substitution x=1 yields the modified formulation due to KANAZAWA 
(the method b) of (*)). 


3. — The calculation of the partition function. 
The partition function Z, of the system of N fermions is 
(20) Z, = ©r(exp[—8#]))=Tr(exp[— PH.) , 


where 6 =1/kT and Tr denotes the trace over a complete set of states of the 
original system, and Tr that over a complete set of states of the extended 
system satisfying the subsidiary conditions (15). The trace is an invariant 
of the unitary transformations, and, therefore, we have also 


(21) Zy pa Tr'(exp [= PÙ) ’ 


where Tr‘ denotes the trace over a complete set of states of the extended system 
satisfying the transformed subsidiary conditions (17). 

4 We shall calculate Z, in an approximation in which we neglect the term 
H, in the exponent. 
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We now build up the complete set of the states over which we sum in (20). 
Let us note that Le depends only on the supplementary variables. There- 
fore, the states of the system satisfying the conditions (15) may be assumed 
as follows 


(22) |ab>, =|x>|b>, , 


‘ 


where |x) is the common eigenvector of the position operators of N particles 


NA REA 


(23) (2724 EEE 
jel 
and 
(24) >>. = Ip = TD od, lb ©) 
lkl<k Ikl <i 


is a vector describing the oscillators. The vectors |b,>, satisfy the conditions 
(29) Li 


We take such vectors |x) in (22) because they are convenient with respect to 
the form of the operator È. 

In further calculations we shall use the common eigenvectors of the ope- 
rators 0, and we assume for them the notation 


(26) Ox [Ven — eer en 


lkl<ky lkl <iq 


We perform now the transformation (11) and after simple calculations we 
obtain 


(27) bX, = È |ad>, = Jin 1Q) dQ <Q — ole) |b». , 


where 


dQ = TT ageag;. 


Ikl<ky 
We see then that 
(28) Tr'(..)= far K0b]... ob, 


() []' denotes the product over k(k # 0) lying in a hemisphere (because CELSO 


k 
QUi= — Qi): 
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where 
N 
dg = [] dx, : 
i=1 


Therefore 


(29) 2%, = Tr’ (exp[— (7+ Hu + A,)))= exp[— PAg] | da<x| exp [87]: 


: | dQ d9'6*(0' — ole) @(Q — o(@)) <Q’ |exp [— BH] |Q>, 


where 
(30) CQ) = TI Ga (0 IO) è 
IEl<ko 
(31) FrodQt) = Order  JralQr) = <Q ld - 
From (25) we get 
; SI BIO 
(32) Jrx(q) = Sta exp | 9a (4) | ’ 
where 
3 1 1 
(33) IZ =, 
(2,70) VaCr 1/7 
Vr 
Thus, we have 
(35) God) = 9(9), (*) 
and 
aye 1 lg)? 
(36) dad <= ov |—5 (4) | 
de 0. 2\C, 


For the mean values of the operator exp [— BT] we have 


Dre SA SA, al 
(37) <ælexp[—BT]|x) = (i v > ee exp [— 2 &,]|det (exp[—if,x,]) |? , 


() We treat gy) as lim Ins assuming here that « is a continuous parameter. 
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where 


hf? 


Oa 


The matrix elements <Q’ \exp[— AH.,.]|Q> have the following structure 


‘ 


(38) <O' exp |= pH,.\|0> = 


= ll On jexp[— BE Oe (0 


exp [— PH 4] |0i 


After short calculations in which we have made use of the summation formula 
for the Hermite polynomials (5) we obtain 


OF LH NOD — 108. Gn) 
(39) | 
ERRO 
where 
(40) h(a. q) = La EE de pet BG? 200") 00] 
delie bon lt + de dd]; 


ba = exp [— phox | . 


Substituting (32), (37), (39) into the integrand in (29) and performing the 
integration over dQ dQ’, we obtain 


(41) Linn = exp[— BA, | 1, x io D. pp S ZA 
fi Sy j=1 


. far idet (exp [— ifix]) |? []” Yio, 


lkl< 
where 
(42) = ee Oe te ae da 
10,9 


Ceo i) 


The final expressions for the partition function for both values of « are as 


(5) E. C. TircHMARSH: Introduction to the Theory of Fourier Integrals (Oxford, 1948), 
p. 77, theorem 53. 
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follows 


N 
(43) Zino = eXP[— BA] (i il p À D Se > exp[— 8 > sg] 
ie Sy JE 


LAIT I= 2 1l—t, 
cp 


È | = oar, 12 Sp use 
Jariaer (exp [— ifix])| ANE ER 


kl<ko 


GARA |, 


IR as 
(44) Bruni = exp[— BA] Gi) vas 2 ae 2. exp[— fp È CAN 


far det (exp[— if,x,]) 2: II" t exp (1 ae : 


lkl< ko 


-a 


In the next Section we shall show that the correct partition funetion is 
given by (44). 


4. — The classical limit #0. 


We shall now consider the transition to the classical statistics, and, there- 
fore, we shall neglect the fact that the quantum particles are indistinguishable 


assuming 
_ i 2 Ù 
(45) WI |det (exp [— ifix]) 2 ~ 
Thus we have 
(46) hee == > erate > exp [— B > Saale 3 
fi in fea J 
where 


(47) Ryo = exp[—fBA,]v TI’ RA = i 
0 es Xp dol ox(æ)0_r(x)|, 


inizia, O ty 
and 
(48) linge ue cee El i fa exp a (1 — tn es(o) 01 (0) : 


It is easy to verify that the classical limit # — 0 with the assumption k, = co 
gives for R,, the expression 


n 


OE N 
(49) Resp E( 5 È Wa) farexp | 


N 
È D > one 5 
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which is exactly the classical configurational integral « das Konfigurations- 
integral »: 


(50) fexp | P > Wil dz. 


The classical limit # > 0 gives for R,, a‘divergent expression. 

We see, therefore, that the correspondence principle requires the modifi- 
cation of the subsidiary conditions of the method of BPMG in the form pro- 
posed by KANAZAWA because the subsidiary conditions in the original form 
lead to meaningless results. 


5. — Further approximation. 


The evaluation of the integrals (49) is very complicated and, therefore, we 
shall perform them approximately. ZUBAREV (°) has given a method for the 
approximate calculation of the classical configurational integral (50) based on 
the expansion of the integrand in powers of 1/VN. This method can be 
applied to the evaluation of R, (we write E, for R,,). We have calculated R, 
by this method up to terms ~1/N with the following result: 


N° N | 
(51) [Rex | b > Wo Di =: w TI exp 
a ZU ‘E 


lkl< ky 


») 


di 


bho i 


Il È ’ \—1 
{ n (1 exp [ bhox]) += i) 7 


RAA SIENS ONE) 


Il metodo delle variabili supplementari si applica al calcolo approssimato della 
funzione di ripartizione di N fermioni interagenti. Si dimostra che il metodo delle 
variabili supplementari nella sua formulazione originale dovuta a Bohm-Pines e a 
Migdal-Galitsky conduce a risultati errati. Si ottengono risultati corretti se si intro- 
duce una modificazione delle condizioni sussidiarie del metodo dovuta a Kanazawa. 


(*) Traduzione a cura della Redazione. 
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Gravitational Motion and Radiation - II. 


A. PBRES 


Department of Physics, Israel Institute of Technology - Haifa 


(ricevuto il 6 Novembre 1958) 


Summary. — General relativistie laws of motion for slow particles are 
deduced from successive approximate solutions of Einstein’s equations, 
subject to de Donder’s conditions. The first three approximations give 
results identical with those of the previous paper () (Newtonian and 
post-Newtonian equations of motion). Higher approximations involve 
contributions from the radiation field. It is found that they can be 
treated independently of other relativistic effects, and cause fifth-order 
corrections to the rest-masses and accelerations of the particles. In par- 
ticular, it is shown that as a consequence of gravitational radiation, a 
system of particles gains energy, thus implying that gravitational waves 
carry away a negative amount of energy. This is made plausible by 
showing that the energy density of the gravitational field is also a nega- 
tive quantity. An electromagnetic analogy of these results is proposed. 


1. — Equations of motion. 


In a previous paper (*) we chose de Donder’s relations as uniqueness and 
compatibility conditions, in order to solve successive approximations of the 
g’” field. It was shown that the laws of motion could then be deduced directly 
from the field quantities without any use of integration around singularities. 

While the principle of this method is extremely simple, its application to 
higher approximations necessitates an enormous amount of computational labour. 
If we are interested only in the laws of motion (and not in the g'” field per se), 
an equivalent (?) and much more speedy method is based on the conservation 


() A. PERES: Nuovo Cimento, 11, 617 (1959). Hereafter referred to as I. All nota- 
tions throughout the present paper are those of I, unless otherwise stated. 
(?) PHAM TAN Hoane: Compt. Rend., 246, 1497 (1958); Nuovo Cimento, 9, 647 (1958). 
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UV 


laws T/°,= 0. As we are dealing with point particles, it is convenient to 
introduce Infeld’s «good » 6 function (*), having the properties: 


[fede gar re Pie aol “n lav=0, 


where f(x) is regular at «= &. 
We have shown, in I, that the «effective source » of the field is the second 
order tensor density 


gT = Y Mov" via —E) 


quo 


rather than the simple tensor density T 
by writing 


Indeed, we get simpler formulae 


(1) (gt), = (T°), + (GT), = 


where 


ap — È ap 
ZI UV 1 “SY U SY To 
ARIE E 
Now 


ji CN J À . = 
(SMS i (Mot) — Move vr sex? = > (Moro, 


Thus (1) reads 


(Mv") + Mv Fr, = 0, 


where the bar over ne means that we have to compute #, at the position 
of the particle which we are dealing with, and to cut-otf all infinite self 
Teams (5) 

Equation (2) gives laws of motion that are perfectly equivalent to those 
obtained with /%, instead of F" 3 (°°). However, as we are working with g”, 
FF"; is more readily computed o Li 

Let us now remember that 


M= mt m + ME. ; DE = gk + Qh + ak + ..., 


(3) L. InreLD and J. PLEBANSKI: Bull. Acad. Pol. Sc., Cl. III, 4, 689 (1956); 
5, SL (1957): 

(4) L. INFELD: Acta Phys. Pol., 16, 177 (1957); Rev. Mod. Phys., 29, 398 (1957). 

(5) E. SCHRODINGER: Space-Time Structure (Cambridge, 1954), p. 55. 
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where m and a are small quantities of the (n+2)-nd order. Let us define 
n n 


DE at + Fh+20 H+ vo Fi. 


m m 
n+2 n+1 m 


Equation (2) is then equivalent to 


pan-1 
; OSE 
(3) n 
=) 
p=n 
(4) mort > m br — 0. 
n n_- PD 
p=0 


The problem of motion has thus been reduced to the determination of 
the db“. 


2. — Newtonian approximation. 


We obtain the Newtonian approximation by solving (3) and (4) for n= 0. 
From (3) we get m= 0, whence, from (4) br 0. It follows 


e 0 
ak == — ne — QU x = 


,m 
x > ja? 
da R 


î.e., Newton’s law of attraction. 


3. — Post-Newtonian equations of motion. 


For n=1, (3) gives 


Va sari oa 10 
m= mb — m È 0, 


AUD == 
whence, we deduce from (4), mde — ORC: 
R (en Ie Oy Wis Cone 
Ges Fy — 20' Fy = 0. 


There are no first-order corrections, either to mass or to acceleration. For 
n= 2, (3) gives 


e as ANDE om sm 
= — 0 — __ 0 2 ) D) = Va 
Di mb M( Fo + vf Fox) m (3 ot DI R 40% dal DI a a 
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Now 

ae 6 (m\ d_m Cl 

re BES gece yk = yk - 5 

RE 2! ce Ri UÈ2R "n2R 
Thus 

m = m oka? + 3- Sse 
di— RI” 

Whence 


1 
/ — 52 
m = gu + 3m > 3 


This result has been discussed in I. 
We then obtain from (4) Mo + mde — 0, whence 


m 
Da = Ce x sn Dal Tk 
ak — pe ae Fi yy” 4 Qo Pi AE 
4 m di 
where 
fa © 1 00, re = 1 ol Ok 
SEA Ep oh E 
2 2 3 3 3 


RER Ok ( ll 
Foo = 39°99, — ato — 219” + gx - 
2 5 


All the needed quantities have already been computed in I. 


4. — The perihelion advance. 
If we explicitely work out the previous formula for the case of a body of 
mass m and velocity v* revolving around another body of mass M, taking 


for simplicity the center of mass at rest, we get, by a tedious but straight- 
forward computation: 


ae MR! e 4m | 2m? È a 5m\ M 3m? (Rw? 
MIE ' M})R er 3 

ne mm, Im? Mo Rv! 

| COMM INR IR EE 
where R* is the vector from M to m. 


If m< M,we have, in ordinary vector notation 


a (1 3 4 M ; 4AMv(R-v 
i A RIA Ra) 
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This equation can now be handley by the usual technique (5). In polar co- 


ordinates, it reads 


2 


M | M 
1 — 4 
ipe 


| 7 — 2 = 


| ve — 4M È Tr. 


x 
This last equation is readily integrated 


726 = lexp[— 4M/r], 


where / is a constant. 
We pass to a new variable wu =1/r, and to 0 (primed) derivatives, instead 


of time (dotted) derivatives (*) and get: 


M ; 
(5) y’ u A + Mu2+ Mw’? + E (exp[8Mu] —1—4Mu exp[8Mu)) . 


| | 


Till now, there have been no approximations, except those explicitely 
stated. 

We now suppose that Mu<1, and replace «', in the right-hand member, 
by its value deduced from w= M(1+ecos0)/P, which is the appreximate 
solution of (5). 

We obtain 


which shows that the perihelion advance, per revolution, is 


670°  67M 
2 Sal See 


where a is the semi-major axis of the «ellipse », and e its excentricity. This 
result is in agreement with that computed by elementary methods (7), and 
with experimental data. This is a good test for our theory, as some other 
approximation methods lead to wrong results in this point (*°). 


(°) H. GoLpstEIn: Classical Mechanics (Cambridge, Mass., 1951), Ch. 3. 

() H. WEyL: Space Time Matter (English translation, New York, 1922), p. 258. 
(8) F. J. BELINFANTE: Phys. Rev., 89, 914 (1953). 

(°) P. Havas: Phys. Rev., 108, 1351 (1957). 
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5. — Radiation field. 


We call radiation field those terms which depend on e. As e is an arbitrary 
constant, they can be treated quite independently of others. We shall denote 
by stars (e.g. g'”) the quantities that are proportional to e; (quantities pro- 
portional to e? appear only in the eighth ‘order, and we shall not have to deal 
with them). 

We have seen in I that ql” vanishes in the 3-rd order as a consequence of 
the law of conservation of mass, and in the 4-th order, because of the con- 
servation of linear momentum. Radiation terms thus first appear only in the 
5-th order of the metric: 


n = i 1 Ra 
que ele > m + — m(rR?)|, 
5 5 2 3! 


* 
* n E fa 
e NOS 
5 5 


gr cas je SA > m{vta!+ var). 
5 


5 


These expressions depend only on time, because, by Newton’s law of at- 
traction 


DI m(R?) = 2 D mo? — Y' DI 


is a space-independent quantity, at this order of approximation. 
It follows that gl” contributes to the curvature tensor R” 
5 


M Only in the 
seventh order, and thus can be annihilated by an infinitesimal co-ordinate 
transformation (1:11), The consequences of this fact will be examined in a 
later publication. 


The next approximation is: 


gee =0, 

6 

* 1 

gue de > |(m mo") "a mar + À = (OER?) SA (mk*)|, 
6 

ge! 0 

6 


(19) L. InreLD and A. E. ScHEIDEGGER: Can. Journ. Math., 3, 195 (1951). 
(11) A. E. SCHBIDEGGER: Phys. Rev., 82, 883 (1951); Rev. Mod. Phys., 25, 451 
(1953). 
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(The origin of the last term of g° is explained in I.) We see that go is 
6 6 


space-dependent, but 


g*,— gt, = ie > mo Re — vR*)=0, 

Ci. A à 
by the law of conservation of angular momentum, which holds in this ap- 
proximation, as a consequence of Newton’s law of motion. It follows that 
CIRO, È 7 
g® also contributes to RZ. only in the seventh order, and therefore can be 


È ‘vot 
annihilated. 

Moreover, as the laws of motion can be written in an invariant form with 
the help of RY (1213), it also follows that radative effects can appear only 
at the seventh order (*) (14). i 

Indeed, it can be directly verified, with the help of the equations of mo- 
tion (3) and (4), that 


Further, we have 


(mR?) meo: corrective terms 
31 ae BI (RES from do ’ 


5 


: m 
= 5 ) 
E 


fox 

0% — 
LIU 
7 


De 4e > m(v* a? + vl a’) = (more!) L xi (vey! R2) ; 


. * A 
where the «corrective terms from g°° » result from 6* = at + ar, 
5 
We must add to this 


which results from the integration of 
* * * 
V2800 — (9% — qk qo 
D Mu ie 


(!2) J. L. SyNGE and A. ScHILD: Tensor Calculus (Toronto, 1949), p. 93. 

(3) F. A. E. PIRANI: Acta Phys. Pol., 15, 389 (1956); Bull. Acad. Pol. al Olle THU E., 
5, 143 (1957); Phys. Rev., 405, 1089 (1957). 

(*) L.e., as fifth-order corrections to masses and accelerations. This corresponds 
to Infeld’s ninth-order equation of motion: 

(4) L. InreLD: Phys. Rev., 53, 836 (1938)- 
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and 


Bh Ap > fi Seno v'R*)] PRA (a RI + a'R*) + i (mR*R?) es (m2) 
to 4 ‘2 2 One 30 asa he 


which is needed in order to fulfil de Donder’s condition do 10, — 0. Note 
that V?3*— 0. 


7 


6. — Radiation reaction. 


From (3), we get 


m mb° fre Lex 00 Gi? 
= — my = —_ m = — È 
5 4 RIO 4. 9 9 )o 


It follows 


* 


d 
Mez em — (31 —17V), 


5 dt 
where 
1 7 A B 
mm 
SS, SE DE FR METTE a 
Gl! == 9 > mu? and i = 5 : 


(Note that 7+V is a constant of the motion, at this approximation). 
Next, we have, from (4) 


* 


‘a Us TES * 
Qk = — © pk 2 Be Fr, 
= m À 7 00 
where 
* * * * * 
(e Me) 1(q0% ol ee Sey 
De AU eg 
6 5 6 5 
and 


All the needed quantities have already been computed. 


7. — Radiated energy. 


We define the radiated energy, per unit time, as minus the rate of work 
of radiative forces 


(6) U=— > motar. 


Lou 
1D 
© 
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di 


We can also compute it with the help of the conservation law (1°) 


[EM = 07 
Here 


i 


UR 


where use has been made of de Donder’s condition, and the explicit form of 
O” is given in I. Thus 


Poe —5, | {[se"+ ar — + [oto : 


Such a computation has already been proposed by LANDAU and LIFSHITZ (1°) 
on the basis of the «linear approximation » of Einstein’s equations, that is, 
by keeping in gt only terms which are quadratic in §’,. This method will 
be correct only if 8“ falls faster than h”’, at spatial infinity. 

That this is not the case can be seen in the following way. Let us sup- 
pose that at spatial infinity gl +h” falls like R-1 The same will be true 
of 6”, (unless the sources are time-independent) and thus ©‘ will fall like R-?. 
It then follows from 


2apv Suv _ uv 
Vr: Sa — OM , 


that 3” will not vanish for large R, but rather behave like log R unless O! 
has a null time-average. This shows that the whole problem of boundary 
conditions in General Relativity should be reconsidered. At any rate, this 
argument may cast serious doubts on the validity of the linear approximation. 


8. — Two-body circular motion. 


We shall explicitely evaluate U in the case of two bodies of masses m 
and M, revolving on circular orbits with constant angular. velocity © = 
—V (m + M)/D*, D being the distance between them. As the final result may 
look paradoxical, and the computations are rather short, we shall quote them 
with some detail. 

First we note that as the kinetic and potential energies are constant (at 
the needed approximation), m, m, ji and qu vanish. Next we choose, for 


(5) L. Lanpav and E. LirsHITz: The Classical Theory of Fields (English trans- 
lation; Cambridge, Mass., 1951), p. 317. 
(16) L. Lanpav and E. LirHsirz: loc. cit., p. 330. 
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convenience, a system of co-ordinates where the center of mass is at rest: 
mor + MV*= 0. 


As a consequence, we can neglect, in Pi all spatial constants which give 


equal contributions to d” and À, as they cal cancel in the computation OLA US 
5 


Thus the only relevant term of po is 4e > m(R*)/5! and gU gives no con- 


ri iti 
tribution. 
Moreover 
I dba mR'R" one MERE" _ mR? 
#00 LE n SIA qu 
HA oi >( RS. i, 3g 3" 2 Ri guy 


The first term vanishes, as a result of the motion being circular, because 
D'D"(a'v + av!) = 0. The second term of 3°, exactly cancels the third term 
of Fe. It remains still — qu, > AOL cara the only relevant term is 

7 
— 4e > m(vtR?)/3! 

Thus 


a ee ae > Mm = (vik) er na (R?°R*)|. 


2 


A straightforward computation, making use of at=—MD"/D* and omit- 


ting all spacial constants of Y, gives 


> TELA 
ie — © m(m + M) zi (221M — 269m), 
5 


30 
whence 
osi 3 m?M?(m + M) 
cb DS i 
For outgoing waves (retarded potential, £= — 1), we obtain that the ra- 


diated energy is negative, i.e., the mechanical system gains energy (D will 
slowly increase). This result is in qualitative agreement with that found 
bit; 


9. — Four-body circular motion. 


If we consider four equal masses placed at the corners of a uniformly 
rotating square, we can expect that there are no radiative effects at this order, 
because the quadrupole moments are constant. This can indeed be easily 


(1?) N. Hu: Proc. Roy. Irish Acad., A 51, 87 (1947). 
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verified by noting that hi k1— 0, and therefore d* reduces to 
5 5 


a = =e>m 2 (© ER?) + 55 (ht R*)|. 


The explicit evaluation of this quantity shows that the contribution of 
the farthest particle is exactly cancelled by that of the two nearer ones, 
De U= 0; 

It may be concluded that the hypothesis that gravitational waves are, 
at the lowest approximation, quadrupole ones, is a reasonable assumption. 
Tf a general proof could be found, the radiated energy would be, in ordinary 
units (7°): 

G 


= kt D 
E 18065 (D ); 


where 


Dr = Si m(3é k El 29! ONE) 2 


This result is extremely striking, because it is exactly the same as for 
electromagnetic quadrupole radiation (1°), if we merely replace m by e and (— G) 
by (1/4xeo), te., if we replace Newton’s law by Coulomb’s law! 


10. — Negative energy. 


We still have to discuss the «unbelievable result » (11) that, as a conse- 
quence of gravitational radiation, our system gains energy. This implies that 
gravitational waves (generated by a system of freely gravitating pole par- 
ticles) carry away a negative flux of energy. As a matter of fact, this result 
becomes quite reasonable if we recall that the energy density of the gravi- 
tational field 


7 


+00 = — (joo — SESIA 00 D0: x 
9 167 0" = 1287 9 9 


2 


is also negative. We can write it, in ordinary units 


oc np 
5 8168 
where g is the gravitational field. Apart from the factor 7, this is the exact 
equivalent of the electrostatic field energy density +(e,/2)E?. If we integrate 


(18) L. Lanpv and E. LirsHITZ: loc. cit., p. 206. 
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the latter throughout the whole space, it is equal to the total potential energy (19). 

In a similar fashion faire 7V. If we add this result to gT®, computed up 
4 

to the fourth order, included, we find 


fonti AOC III METE 


i.e., the total (rest + kinetic + Newtonian potential) energy. 
The electrostatic analogy still holds, although it is less trivial than in the 
previous section, because rest-mass is not conserved, while charge is. 


I am indebted to Professor N. Rosen for having suggested this problem 
to me, and for stimulating discussions. 


(9) W. Panorsxy and M. PaILLIPS: Classical Electricity and Magnetism (Cambridge, 
Mass., 1955), p. 88. 


135 ASS AN IPO) 6) 


Dalle successive soluzioni approssimate delle equazioni di Einstein vincolate alle 
condizioni di de Donder si deducono leggi di moto per particelle lente secondo la rela- 
tivita generale. Le prime tre approssimazioni danno risultati identici a quelli del pre- 
cedente lavoro (1) (equazioni di moto newtoniane e postnewtoniane). Nelle approssima- 
zioni superiori compaiono contributi dal campo radiativo. Si trova che questi possono 
trattarsi indipendentemente da altri effetti relativistici e producono correzioni delle 
masse a riposo e delle accelerazioni delle particelle. In particolare si dimostra che in 
conseguenza di radiazione gravitazionale un sistema di particelle acquista energia, sot- 
tintendendo così che le onde gravitazionali asportino una quantità negativa di energia. 
Il che si rende plausibile dimostrando che anche la densità d’energia del campo gravi- 
tazionale è una grandezza negativa. Si propone per questi risultati un’analogia elettro- 
magnetica. 


(*) Traduzione «a cura della Redazione. 
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The Construction and Working 
of a Wide Line Nuclear Magnetic Resonance Spectrometer 
and the Measurements of Some Chemical Shifts. 


S. S. DHARMATTI, K. J. SUNDARA RAO and R. VIJAYARAGHAVAN 


Tata Institute of Fundamental Research - Bombay 


(ricevuto il 10 Novembre 1958) 


Summary. — A wide line nuclear magnetic resonance spectrometer of high 
sensitivity operating between 4 and 10 MHz, has been described together 
with a high stability electromagnet system. Using the spectrometer, the 
chemical shifts of 170 resonance in several compounds have been measured. 
The results have been interpreted in terms of CO bond character in 
various groups of organic compounds. 


1. — Introduction. 


The nuclear magnetic resonance spectroscopy is based on the fact that 
the various isotopes of the elements can be distinctly identified because of their 
differing nuclear gyromagnetic ratios. Each isotope, with non-zero spin, pos- 
Sesses a magnetic moment, and behaves as though it were a tiny magnet with 
a unique strength. So, if matter be placed in a steady magnetic field H,, say, 
in the Z-direction, the nuclear magnetic moments will tend to orient them- 
selves parallel to the field. After the establishment of thermal equilibrium, 
there will be a net paramagnetic polarization in the Z-direction given by 


(1) My = 7H, , 


where y is the nuclear magnetic susceptibility. If now, we superimpose an 
oscillating magnetic field in the «-direction such that 


(2) H, = 2H, cos wt, 
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the polarization vector starts deviating from the Z-direction and reaches the 
maximum deviation when the Z-field H, approaches the resonance value H 
given by Larmor’s relation 


(3) y|H=@, 


4 


where y is the gyromagnetie ratio of the nucleus under consideration and © 
is the precessional frequency. It is well known that the resonance phenomenon 
was simultaneously discovered in 1946 by BLocH and his collaborators at 
Stanford (1) using the « crossed-coil » method, and PURCELL and his co-workers 
at Harvard (?) employing the bridge method. We have essentially followed 
the technique developed by BLocH and our spectrometer is similar in prin- 
ciple to those of PRocToR (*) and WEAVER (*). 


2. — Description of the spectrometer. 


The spectrometer consists of an R.F. transmitter, a probe, an R.F. pre- 
amplifier, an R.F. amplifier, a detector, an audio frequency amplifier and 
an oscilloscope. To record very weak signals, a narrow band amplifier with 
a phase-sensitive detector and a graphic recorder have also been included. 
The steady field of several kilogauss is produced by an electromagnet with its 
associated power supply and current regulator. 

The transmitter (4) consists of an R.F. oscillator operating between 4 and 
10 MHz and a power amplifier which produces enough current for an ade- 
quate R.F. field in the transmitter coil. A 6SN7 double triode serves as an 
R.F. oscillator which drives enough current at the grids of the 832A power am- 
plifier. With — 45 V as grid bias and 300 V on the plates of the 832A tube, 
by varying the screen-grid voltage from —10 V to +30 V, one can get R.F. 
fields in the transmitter coil ranging from 100 milligauss to 2 gauss. The 
power amplifier serves to isolate the transmitter coil from the load variations 
as also to supply adequate and variable R.F. current to the coil. A twin lead 
co-axial line feeds the R.F. power to the probe which is kept inside the gap 
of the polefaces of the electromagnet. As the transmitter and the receiver 
are located close to each other, adequate precaution is taken to minimize the 


(1) F. BLocH: Phys. Rev., 70, 460 (1946). 

(2) E. M. PurcELL, H. C. Torrey and R. V. Pounp: Phys. Rev., 69, 37 (1946). 
(3) W. G. Proctor: Phys. Rev., 79, 35 (1950). 

(4) H. E. WEAVER: Phys. Rev., 89, 923 (1953). 
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unwanted direct coupling between the two units which will result in an ex- 
cessive leakage voltage at the detector. The transmitter is mounted in à 
rectangular brass box which is fixed rigidly to a platform to minimize micro- 
phonics. The regulated power supplies, which serve the two units, are ade- 
quately decoupled to reduce interference between them. 

It is needless to say that the probe or the « Head » is the most delicate 
and sensitive part of the spectrometer. The construction of the « Head » is 
similar to that of Weaver’s (4). It contains two coils mounted orthogonally 
and rigidly on a perspex frame. The transmitter coil is of the Helmholtz type 
and one half of the coil can be moved relative to the other coil known as re- 
ceiver coil. The receiver coil holds the test tube containing the sample and 
is electrically shielded by a Faraday cage from the transmitter coil. In ad- 
dition to the two usual controls known as U and V mode controls which 
effectively decouple the two coils and introduce a leakage of the desired mode 
at the detector, there is also an additional arrangement by which the receiver 
coil can be moved vertically up and down relative to the transmitter coil. 
The reasons for such an elaborate system of controls are apparent from the 
following considerations. 

For an R.F. field of the order of 0.5 G, about 50 V exist across the trans- 
mitter coil. If the two coils can be mounted exactly perpendicular to each 
other, then there will be no induced voltage in the receiver coil. But there 
is always a practical limit to which this accuracy can be achieved and usually 
the tuned receiver coil may pick up voltage as high as 5 V. But, the maximum 
allowable voltage across this coil is of the order of a few millivolts, as other- 
wise, this voltage, after amplification, will saturate the other stages of the 
R.F. amplifier, resulting in the loss of the efficiency of the spectrometer. By 
moving either the transmitter or the receiver coil relative to each other, with 
the help of a differential screw arrangement, considerable decoupling can be 
achieved. ‘The remaining voltage can be removed by the operation of the 
U and V mode paddles. The V mode paddle (5) consists of a semicircular 
copper dise of about 0.5 mm thickness while the U mode (*) control consists 
of a resistance and an inductance in series in the form of a loop. By a judi- 
cious operation of the two controls, complete decoupling can be effected, and 
any desired mode of leakage can be introduced. This leakage serves to make 
the 1N34 crystal diode detector work efficiently and also determines the shape 
of the nuclear induction signal. The whole assembly is mounted rigidly in 
a brass box on whose outer covers are pasted two ten-turn coils of Helmholtz 
type which are supplied with an audio-frequency current to modulate the 
nuclear signal. 


(5) F. BLocH, W. W. HANSEN and M. E. PACKARD: Phys. Rev., 70, 474 (1946). 
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or 
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The receiver, which consists of an R.F. preamplifier, R.F. amplifier, a de- 
tector and an A.F. amplifier, is shown in Fig. 1. The detected nuclear signal 
is enhanced by a factor of Q (about 80 in our case) of the receiver coil using 
a tuned circuit. This voltage is further amplified by a low noise amplifier 
employing Wallman’s circuit which has a gain of about six. The amplified 
signal is fed into a six stage R.F. broadband amplifier having a total gain of 
fifty. The modulated R.F. signal is finally detected with an 1N34 crystal diode. 


Preamplifier RF Amplifier (total six such units) 
6BQ7 64K5 Demodulator 
06 27k& 
PAR ane 1N34 
| five 1 
RF amp! 
0.01 | 001 
L BAKE ye pF 
470k92 dr ie 
ee 33k 
Audio an fier Cathode -follower | ok: 
+ 300V 2 6SL7 
470 
kS2 
== 
To cho and 
narrow -band 
amplifier 


8 

IM& $ 10003 = 
ls n) UF 
Fig. 1. — Receiver. 


The detected audio signal is further amplified and applied to the Y plates 
of a 304A DuMont oscilloscope. The X-axis plates are supplied with a voltage 
derived from the modulating coils. 

Strong signals, which have a large signal to noise ratio, can be displayed 
thus on the oscilloscope. However, many signals cannot be displayed on a 
scope as they will be considerably less than the noise of the system itself. 
Such signals can still be observed by using a « lock-in-amplifier » with a twin- 
T-filter which rejects all other than the fundamental frequency of modulation. 
We have adopted a circuit similar to the one used by BLOEMBERGEN (°) as 


(6) N. BLoeMBERGEN: Nuclear Magnetic Relaxation (The Hague, 1948). 
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shown in the diagram of Fig. 2. The component of the nuclear signal is am- 
plified in the first two stages and is applied to the grids of the two pentodes, 


used as a mixer with a mixing voltage of proper phase and the same modu- 


65/7 6SL? 65/7 (2) 6SN 7 


16 pF <100K 


470K 22M 


|—e Mixer 10V RMS 37~ 


” F- Filter (Twin-Tee) 


Fig. 2. — Narrow band amplifier with phase sensitive detector. 


lation frequency. This eliminates all spurious signals and noises. The detected 
voltage is amplified and smoothed by an R.C. filter. Then it is applied to 
the grids of the twin triode and the signal is recorded in the D.C. milliammeter 
which is connected across the cathodes of the twin triode instead of across 
the plates. We found that this modification resulted in a great increase in the 
efficiency in detecting weak signals. The magnitude of the current is propor- 
tional to the amplitude of the nuclear signal and its sign depends upon the 
phase of the input voltage. By utilizing this method within practical limits 
a large improvement in the signal to noise ratio (between four and five in the 
case of a 2D signal in 0.5m MnSO, solution) can be realized. 


3. — The electromagnet system. 


The electromagnet (*) used with the spectrometer has pole pieces about 
9.5 in. in diameter and an air gap of 1.75in. For a typical sample volume 
the homogeneity is better than 1 in 10° at the most favourable spots. In 
Figs. 3 and 4 we have shown the signals of ?D in a40% D,O sample and of 


() We are greatly indebted to Prof. Hans H. Sraug of Physikalisches Institut der 
Universität, Zürich, for providing us with an electromagnet constructed by Wild 
Muri Co., Switzerland, under his personal supervision. 
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Fig. 3. — 2D signal in 40%, D,0 sample. 


Fig. 4. — ?*Na signal in saturated NaOH solution. 
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23Na in saturated NaOH solution with the wiggles taken at the fields of 8 kG 
and 5 KG respectively. There are two water cooled windings around each pole 
piece. Each section consists of about 20000 turns and the magnetization 
arises from a large number of turns rather than from a large current. This 
feature makes it possible to regulate the current by hard tubes, specifically 
six 304 TL in parallel. The diagram showing the, rectifier, the current regu- 
lator and the D.C. amplifier is given in Fig. 5. A maximum field of 15 kG 


304 TL 
(6 Valves) 


Magnet coil 1 


400V to 4000V 


transformer 872A (6 Tubes) 


1508 


standard 
resistance 


+4 kV 
304 7L di 
cathode DPDT switch 45V 
A45V ti 
150 $100V= 25V We 
i 30 2 
304 7L 


12AX7 6AU6  VRIOS 


Fig. 5. — 4 kV magnet power supply and current regulator. 


is obtained for a current of 2A at 4000 V. This voltage is derived by a 
step-up transformer by full wave bridge rectification from 400 V 3 phase lines 
using six 872 mercury rectifiers in parallel. To obtain current regulation, the 
magnet current is passed through a fixed resistor of high wattage which is 
immersed in an oil container. The difference in voltage across the resistor 
and that of a standard battery is amplified by a D.C. amplifier of two stages 
with a total gain of 2000. The output is applied to the grids of 304 TL 
tubes, thus providing a feed back capable of reducing the line voltage varia- 
tions of about 10% to current variations of 0.001% approximately. Thus 
an overall stability of about 1 in 10° is achieved in a typical operating range 
of 3 KG to 11kG. The magnetic field, after three hours of warming up, does 
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not drift by more than 500 mG in an hour. The steady field is modulated 
by passing an audio frequency current of 37 Hz to the coils cemented on the 
plates of the nuclear induction head. This current is obtained using a Hewlett 
Packard 206A model audio oscillator with a power amplifier, similar to the 
one used by BLOEMBERGEN (°). For a current of 0.5 A in the sweeping coils, 
the modulating field is of the order of 5 G which facilitates the presentation 
of the signals on the oscilloscope. It is possible to vary the sweep field from 
5G to 20mG. The mixer voltage, which is derived from the unit, is applied 
to the phase sensitive detector of the narrow band amplifier. 

The power supplies for all the units of the spectrometer are well reguiated 
and have a ripple content less than one millivolt r.m.s. The D.C. filament 
supply is obtained from a bridge rectifier. Though it has a ripple of 10 mV 
r.m.s., it does not contribute much to any pickup. This unit is more con- 
venient than the batteries which require constant attention. 


4. — The working of the spectrometer. 


Usually liquid samples are preferred because they often lead to narrow 
lines, and the relaxation times can be controlled by the addition of para- 
magnetic salts like manganese sulphate or ferric nitrate. To obtain suitable 
conditions for the signal observation, the magnitudes of the sweep and the 
R.F. fields must be chosen carefully. Usually the R.F. field is kept about 
0.5 G (2H,) and the sweep field is varied from 5G to 0.1 G depending on 
whether the signal is to be observed on the oscilloscope or recorded on the 
charts of the D.C. milliammeter. The total receiver gain is normally kept 
around 100 and the leakage at the detector is about 0.2 V for the efficient 
operation of the diode. The twin-T-filter in the narrow band amplifier has 
a band width of about 1 Hz and we have used total output time constants 
of 10, 20, 40 and 60s. The frequency of the R.F. oscillator can be varied 
from 4 to 10 MHz and the frequency is measured very accurately using a 
Hewlett Packard frequency meter. The magnetic field is changed with the 
help of a 50 kQ helipot either by hand for a quick change or by a motor at 
the rate of about 1 G per minute through a range of 2000 G for a slow and 
steady change. A somewhat lower rate of field scanning can be achieved by 
employing a step down gear arrangement for the motor drive. 


5. — The efficiency of the spectrometer. 
There are two nuclei whose signals serve as severe tests of the optimum 


efficiency of the spectrometer. They are 2D and 170 which occur in nature 
in such low abundances as 150 and 350 parts per million respectively in or- 
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dinary water. Both the nuclei possess small quadrupole moments, which due 
to interaction with the molecular fields, should attain favourable relaxation 
times. However, we found it necessary to prepare à 0.5 molar MnSO, solution 
to record a good signal of ?D in 3 em? of water at 8.5 KG. For the case of the 
170 nucleus, we used ordinary water with- 
out any paramagnetic catalyst and obtained 
a signal‘of 170 with signal to noise ratio of 
about 20 in a field of 9 kG. The two signals 


= vie 


a) b) 


Fig. 6. — a) The derivative of w-mode signal of ?D in water containing 0.5 M MnSo, solu- 

tion. R.F. field of 0.3 G (H;) and sweep field of 0.4 G (H,) are used. b) 2D signal 

with audio gain reduced to + of that used in the recording of Fig. 6a. The other para- 
meters are the same. 


lef 


are Shown in Figs. 6 and 7. An equally strong signal of 170 was recorded 
using concentrated nitric acid as sample. The absorption mode signal of 170 
in concentrated HNO, is shown in Fig. 8. 


== 


Fig. 7. — 170 signal, the derivative of 
u-mode, obtained in 3 em° of ordinary 


Fig. 8. — 170 signal, the derivative of 
v-mode, obtained using concentrated 


water. R.F. field of about 0.25 G (H,) 
are used. Audio gain is about + of that 
used in the tracing of Fig. 6a. 


HNO; as sample. The other parameters 
are the same as those used in the re- 
cording of Fig. 7. 


To test the overall efficiency of the spectrometer, we studied the tempe- 
rature and the field dependence of *Co resonance shifts and compared the 
results with those obtained by another group in our laboratory using the 
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Varian spectrometer (7). PROCTOR and YU (*) studied the resonance of 5Co 
in several cobalt complexes and found the chemical shifts as large as 1.3% 
between two compounds such as potassium cobaltioxalate and potassium 
cobalticyanide. It was explained by RAMSEY (°) that there exist several ex- 
cited electronic energy levels close to the ground state and that there should 
therefore be a slight temperature dependence of the chemical shift in 5°Co 
resonance. We studied four compounds, viz., CO(NH,).Cl;, [Co(NH3);CO;]N03, 
Na,;CO(NO,), and [Co(NH,);ClCl, at a resonance field of 5.7 kG for the tem- 
perature dependence of the chemical shifts. These samples were heated to 
85 °C and the signals were observed on the calibrated screen of the oscillo- 
scope at this temperature. The resonance field of each sample at room tem- 
perature (27 °C) was observed immediately on the oscilloscope screen by re- 
moving the hot sample and replacing it with the ordinary sample. The dif- 
ference in field for each sample is given below: 


TABLE I. 
Sample Difference in field 
(Hy, [Coram Ag; °c) 
| Na;CO(NO:)s + 0.8 gauss 
[Co(NH;);CO,INO; + 0.7 » 
[Co(NH,);C1JCI, + 0.75 » 
Co(N'H3)gCl (8) + 0.7 » 
(*) PRocToR and Yu (°) reported a shift in frequency of 
about 1 kHz for the same difference in temperature at the 
resonance frequency of about 7000 kHz for Co(NH3)s Cl 


solution. 


As the difference in the field for each sample was observed in the calibrated 
oscilloscope, the error is of the order of 100 mG. Within this experimental 
error, our results are in good agreement with those of PROCTOR and YU (8). 

The field dependence of the resonance was checked in a few compounds at 
different fields such as 5.7 kG and 7.8 KG. The measurements at the latter 
field were done using the Varian V4200B spectrometer. The results of the 
two experiments are given in terms of the dimensionless quantity o defined 
in Table II. They are in excellent agreement showing thereby an extremely 
satisfactory working and the high fidelity of our spectrometer. | 


(?) S. S. DHARMATTI, C. R. KANEKAR and S. C. MATHUR: Proceedings of Interna- 
tional Conference on the Peaceful Uses of Atomic Energy, 1958 (to be published). 

(8) W. G. Proctor and F. C. Yu: Phys. Rev., 81, 20 (1951). 

(9) N. F. Ramsey: Phys. Rev., 86, 243 (1952). 
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TABLE II. — Chemical shifts of 5°Co in cobalt compounds. 
lea Shiftin G|ShiftinG] gx 108 
No. Compound at 5700 G | at 7800 G Remarks 
1 2 1 2 
1 | [Co(NH,),(NO,),JNO,| + 6.8 + 9.35 | +11.9 | + 12.0| Cis-isomer 
2 | [Co(NH,),(NO,).|NO3 + 5.20 + 7.25 | + 9.1 | + 9.3 | Trans-isomer 
3 | [Co(NH,),;(NO,)]CL + 3.5 4.8 + 6.1 | + 6.2 | Yellow solution 
| 4 | [Co(NH;);ONO]CI, + 3.5 + 4.8 + 6.1 | + 6.3 | Orange solution 
x + 0.35 — = 0:6 — Weak signal 
5 | NasCo(NO,); +415 | +57 | + 7.3] + 7.3| Strong signal 
6 | Co(NH;);Cl 0 0 0 0 Reference sample 
7. | [Co(NHx);CO;JN0; also It Ser ost, to 
8 | [Co(NH;);CIICL — 5.6 — 7.7 — 9.8 | — 9.9 
0 — 5.7 — 7.7 —10.0 | — 9.9 | Strong signal 
9 |[Co(NH;),COsINOs == 8.7 vis =12.0 153 05.4 | Weak sienal 


MATTI et al. (*). 


1) Chemical shifts measured using our spectrometer. 
2) Chemical shifts measured, using Varian’s V-4200 B wide line spectrometer by DHAR- 


6. — Chemical shifts of YO resonance. 


o= (H;— H,)/H,; where H; is the resonance field for a given compound and H, is the 
resonance field for the reference sample. 


The chemical shifts of 170 resonance were first reported by WEAVER et al. (1) 
who investigated several compounds in a field of 10 kG using the Varian 4200B 
wide line spectrometer. They found the values of chemical shifts ranging from 
200 mG to 6.9 G in a field of 10 KG. Barring water and concentrated nitric 
acid, the signals of YO in all the other compounds are very weak and the 
detection of such weak signals offered us an excellent opportunity for testing 
the sensitivity and performance of the spectrometer. We extended the study 
of 170 resonance shifts to many organic compounds in order to interpret some 


of the shifts in terms of the bonding of oxygen in these compounds. 


The 


(2) H. E. Weaver, B. M. Totpert and R. LaForce: Journ. Chem. Phys., 28, 


1956 (1955). 
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shifts in the compounds studied thus far in a field of 9 KG, are given in terms 


of o in Table III. The organic compounds used in the experiment are the 
products of either B.D.H. (AnalaR) or Merck (extra pure). 


TABLE III. — Chemical shifts of O0 resonance. 


No. Compound alo Remarks 
IR EM eh y]talcoho oa.) sa + 0.400 | = | 
OE: ya lcoho NN 0 98% pure 
SO VV ACT RI o 0 Reference sample 
Fe RNA ER O 0 — 
5 RFOPIO DIC, GLAS SE — 2.1 — 
6 | AGG: gol 45 RE — 2.2 — 
7 HOENILCRA GIO SR e a 6 6 — 2.7 | — 
RARE ro plonie anhydride yas a. — 2.4 Stronger signal 
| — 4.2 Weaker signal 
9 Acetyl acetone ...... — 2.35 | — 
10 EthylMacetate CN — 3.1 — 
11 Concentrated nitric acid . . — 4.6 | Very strong signal 
12 lNecAcetoney i a cet eer — 6.0 — 
13 | Methylethyl ketone . . . . — 6.0 — 
14 Aceto acetic ester. . . . . — 6.2 | Broad signal 
15 | Ethyl malonate . . . . . . — 6.25 — 


The resonance of 170 in methyl and ethyl alcohols and diethyl ether is 
obtained close to that of water with the signal of 170 in methyl alcohol oc- 
curring at a slightly higher field (about 350 mG at 9kG). There is no measur- 
able chemical shift of the 170 the resonance in the other three compounds. These 
results indicate that, in compounds, in which oxygen forms a single bond with 
carbon or hydrogen, the resonance of 170 occurs at very nearly the same field. 

It is interesting to note that the resonance of 170 in the two ketones studied 
(namely CH,COCH,; and CH,COC,H;) occurs on the other hand at a consi- 
derably lower field than those observed for the aliphatic alcohols. In the 
ketones, the oxygen forms a double bond with carbon in which p-orbitals are 
used for the formation of both x and o components. This would result in an 
increase in the effect due to the second order paramagnetism for the double 
bonded CO group as compared to that for the single bonded CO group in the 
alcohols. Hence one would expect the shielding of 170 to be lower in the 
ketones than in the alcohols. No measurable chemical shift has been found 
between acetone and methyl ethyl ketone. This is probably because, the sub- 
stitution of the ethyl group in place of the methyl group in the latter ketone, 
does not affect the electronic environment of the strong C = O bond appreciably. 
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We have observed only one 170 signal for the carboxylic group in formic 
acid, acetic acid and propionic acid. This indicates that, although this group 
contains two oxygens, they are structurally equivalent, owing to the possi- 
bility of a resonance between the following structures (11) (neglecting the di- 
merization of acids due to hydrogen bonding), 


AO O- 
BR UC R CCE 
OH OH 


This resonance gives the CO bond in the carboxylic group a partial double bond 
character. The observed bond distance for CO in the formic acid dimer (1?) 
(1.29 À), which lies between that of a single and double bond (1.42 À and 1.22 / 
respectively), confirms this. It may be noted here that the !0 resonance 
for this group also occurs at fields intermediate between those for the 
alcohols and the ketones. 

Resonance similar to the carboxylic group is also possible in esters of these 
acids, giving a partial double bond character to the CO group in esters. The 
chemical shift in ethyl acetate confirms the above reasoning. 

We observed two signals of 170 in propionic anhydride, with an intensity 
ratio of about 2:1, the stronger signal lying close to the fields of the acids 
while the weaker occurs at a field intermediate between those corresponding 
to the acids and the ketones. In the anydride, the following type of reso- 
nance is possible (11). 


40 /O 10, 
CH;_—CH,—CK CH,—CH,—C¢ CH,—CH,—0Z 
O Ot 


Oe: 
(neu, a CH,—CH.—OC 


So it can be seen that there are two types of oxygens in this anhydride 
with one oxygen exhibiting a stronger double bond character than the other two. 

The resonance of 17O in acetylacetone is obtained close to the fields of 
resonance for the acids. It has been found by the method of bromine titration, (1?) 
that about 80% of this compound is in the enol form 


| | and the rest in the keto form | | 
OM O O 


(11) G. W. WHELAND: The Theory of Resonance (New York, 1945). 

(12) L. PAULING: The Nature of Chemical Bond (New York, 1940), p. 306. 

(18) E. A. BRAUDE and F. C. NacHOD: Determination of Organic Structures by Phy- 
sical Methods (New York, 1955). 
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ScHOOLERY et al. (14) using the proton resonance estimated that about 
85% is constituted in the enol form with 15% in the keto form. Our result 
clearly indicates that a substantial part of this compound is present in the 
enol form. We could not, however, detect the signal due to the other form 
as it would be very weak. 

A preliminary investigation shows, on the other hand, that the resonance 
of 17O in aceto-acetic ester and ethyl malonate are obtained at fields close to 
the fields of resonance for the ketones. This is understandable since the two 
compounds exist mostly in the keto form (14°) 


O 


CH;—_C—CH,—C—0C0,H; C,H,0—0—CH,—C—OC,H, 
| | and | | 
O O 


respectively with a slight percentage of enol form 


CS ha ead ine A GIA se BUSH: 
an 
OFy. HO O PPS O. 


respectively. 

However, it is possible to change the relative strength of each form and 
then observe the signal intensities of 170 due to both the forms. A check, on 
the amount of conversion, can then easily be made either by observing the 
proton resonance using a high resolution spectrometer, or by the conventional 
method of bromine titration. Investigation along these lines is being pursued. 

It can be concluded that the study of the 170 resonance yields valuable infor- 
mation about the bond character between carbon and oxygen in the organic 
compounds. In the case of single bond, it appears that the contribution to 
the chemical shift is predominantly due to diamagnetic effects of nearly the 
same magnitudes thus placing the resonance of 170 at the highest field in the 
compounds so far studied. In the case of pure double bond, the second order 
paramagnetic effect is quite pronounced (due to the p character involved) and 
sometimes even overrides the overall diamagnetic effect; in the partial double 
bond case, the net shielding is an average of the two effects. 

While we are still continuing the study of O resonance shifts in other 
organic and inorganic compounds, we are also interested in their line widths 
which may throw more light on the structure of these compounds. 


(4) J. N. ScHooLeRY, H. S. JANETT and M. 8. CADLER: Journ. Chem. Phys., 
21, 2092 (1953). 
(5) H. B. Watson: Modern Theories of Organic Chemistry (London, 1946), 2nd. ed. 
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RASS S URNERON() 


Si descrive uno spettrometro ad ampia risonanza magnetica nucleare ed alta sen- 
sibilita funzionante tra 4 e 10 MHz, ed inoltre un sistema di elettromagneti di alta 
stabilità. Coll’aiuto dello spettrometro si sono misurati in alcuni composti gli sposta- 
menti chimici della risonanza del 170. I risultati si interpretano in termini del legame CO 
in vari gruppi di composti organici. 


(*) Traduzione a cura della Redazione. 
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K*-Deuterium Scattering - I (*). 


M. GourDIN and A. MARTIN 


Laboratoire de Physique, Ecole Normale Supérieure - Paris 


(ricevuto il 6 Dicembre 1958) 


Résumé. — Nous avons étudié les diffusions élastique, inélastique et 
d’échange méson K+-deutérium à l’aide de l’approximation «impulse ». 
Dans ce calcul préliminaire, nous avons négligé l’interaction entre les 
nucléons dans l’état final. Des calculs numériques sur machine électro- 
nique 704 IBM ont été effectués dans le cas particulier où les seules 
amplitudes S méson K+-nucléon interviennent. Nous trouvons que la 
distribution angulaire de la diffusion d’échange et le rapport des sections 
efficaces d’échange aux sections efficaces sans échange sont les seules 
quantités qui peuvent permettre de déterminer facilement la longueur 
de diffusion dans l’état de spin isotopique [= 0 pour les réactions méson 
K+-neutron, quantité actuellement mal connue alors que la diffusion 
méson K+-hydrogène a permis déjà d’évaluer l’importance de l’état de 
spin isotopique J = 1. 


1. Introduction. 


The K*-nucleon interaction appears in three basic processes: elastic scat- 
tering on protons and neutrons, exchange scattering on neutrons: strangeness 
conservation forbids processes involving hyperons. Since the isotopic spin of 
K° is ;, the three corresponding scattering amplitudes can be expressed in 
terms of two scattering amplitudes in the total isotopic spin states J = 1 
and J=0, and of the Coulomb scattering amplitude f, (1): 


(I) Kt +p +>Kt4+p a, + fe 
(II) K*+n>K*+n (Sini 
(III) KP sen Keep (Ce See 


(") Supported in part by the United States Air Force through the European Office, 
Air Research and Development Command. 


(*) It can be easily seen that coulombic phase factors must be neglected. 
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At the present time, our knowledge of this interaction is based on two 
types of experiments: 


1) Scattering of K* by hydrogen which gives information on the J=1 
interaction. The results show a relatively constant total cross-section up to 
200 MeV corresponding to a scattering length of the order of 0.24 h/ue (?). 
The angular distribution shows the predominance of the S state in spite of 
an unexplained marked minimum in the backward direction. Obviously the 
experimental data can be improved and make the situation clearer. 


2) Elastic and inelastic scattering of K* in emulsions. — The description 
of elastic scattering by an optical potential shows that the average K* nucleon 
interaction is repulsive (*). The ratio of exchange to non exchange inelastic 
(or rather quasi-elastic) scattering of K* should directly give the ratio of the 
I=0 and 1 — 1 scattering amplitude if the nucleons, inside the nucleus, could 
be considered as perfectly free. This is not the case and we feel that the cor- 
rections to be done before giving the J = 0 amplitude are rather uncertain. 
Nevertheless it turns out that the 7 = 0 amplitude is small and seems to in- 
crease with energy. 


The need for another, more direct, way of investigation of the 1 —0 
interaction seems to be clear. In our opinion, the experimental data which 
can be handled by theoreticians are those given by scattering of K* on the 
simplest substance containing neutrons, i.e. deuteron. The experiment is quite 
feasible, since analogous experiments have recently been done by sending 
K” mesons in a deuterium bubble chamber (¢). 

The object of the present paper is to give a preliminary calculation of K*- 
deuteron scattering, both elastic and inelastic; this calculation uses impulse 
approximation; final state interactions have been provisionally neglected, and, 
in the numerical computation, done on IBM 704 digital computer, P-waves 
have been ignored. We have especially investigated the sensitivity of the 
results to the /= 0 scattering length. It turns out that the angular distri- 
butions for both elastic and inelastic non exchange scattering are dominated 
by form factors and Coulomb scattering; therefore it would be difficult to see 
a modification of this angular distribution due to P waves. The most inte- 
resting quantities are the differential and total exchange cross-sections. Their 
ratios to the non-exchange cross-sections are very sensitive to the 7 —0 


(2) J. E. LANNUTTI, S. GOLDHABER, G. GOLDHABER, W. CHUPP, 8. GIAMBUZZI, C. 
MarcHI, G. QuARENI and A. WATAGHIN: Phys. Rev., 109, 2121 (1958). 

(3) G. Ico, G. RAVENHALL, J. TIEMANN, W. CHUuPP, G. GOLDHABER, J. LANNUTTI 
and R. M. THALER: Phys. Rev., 109, 2133 (1958). 

(4) R. D. TRIPP: Rapport Conference Genéve (1958). 
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amplitude and the corrections to be done to extrapolate to the case of free 
nucleons are probably little affected by the approximations made, especially 
in the high energy region. 


2. — Differential cross-sections. 


We bombard the deuteron in its ground state with a given distribution of 
momentum for the nucleons, by a K* meson of momentum k in the labo- 
ratory system. We assume that the collisions can be treated by the impulse 

approximation where the nucleons 


k ki =k, k are considered as free. It is pos- 

x N N à N sible to characterize the scattering 
+, Sk ki in terms of K*-nucleon amplitudes. 
P ET QU PUR We neglect the binding energy 

ke k, kt of the deuteron and we have enti 

Initial state Final state 1 Final state 2 mated the multiple scattering 
Fig. 1. — K'-deuteron scattering in the effects and shown that they are 
impulse approximation. smaller than in the case of x 


meson-deuteron scattering (5). 
Let us now consider the elementary K* meson-nucleon processes. The 
amplitudes corresponding to (I), (II) and (III), can be respectively written: 


| . . . . / . 
Arno D A,+0,:B,; A,+0,:B,, 
where the spin of the involved nucleon explictly appears. We can exhibit 
three interesting cases: 


a) The elastic scattering where the deuteron is not dissociated by the 
collision. 


b) The inelastic scattering without exchange. 


c) The exchange scattering where the emerging meson is a K° meson. 
The final state is a p-p system obeying the Pauli exclusion principle. 


After summation over the final spin states of the nucleon and average over 
the three initial spin states of the deuteron, we obtain in terms of the three 
SE : : gi 2 = EST 4 : ee 
K'-nucleon amplitudes, expressions for the cross-sections in the laboratory 
System corresponding to processes a), b) and e). 


(5) S. DreLL and L. VERLET: Phys. Rev., 99, 849 (1955). 
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1) Elastic cross-section. 


do 


9 | A 
ee À fy 1| F@ 2d cos Oou 
(1) n cee ts | = |B, + B,|? Pp | 


d così J) F(0)| dcosì ’ 

q is the transfer moment k’— k of the K* meson; 0 the angle between k and k’ 
in the laboratory system and F(q) the form factor defined as the Fourier 
transform of the square of the modulus of the deuteron radial wave function. 
In practical calculations we use the Hulthén wave function to describe the 
deuteron (see Appendix). 


2) Inelastic cross-section without exchange. — Assuming no final state inter- 
action we obtain in the laboratory system: 


(ROUE 0) {AS +|A,/?+ |B, |? + |B, |?) (p, q) + 


=F AS Ag à ASA; ae 5 (B,B; * B,B;)iG(», È k'?dk'd(cos 0) 2 


The auxiliary quantities p and q are defined by the following equations: 


‘Pi ua 
(3) F 
gy ikke ee NE) 


The last equality is a consequence of momentum conservation. The form 
factors F(p,q) and G(p,q) are related to the radial wave function of the 
deuteron; they are explicitly given in the Appendix together with the nor- 
malization constant C. 


3) Exchange cross-section. 
(4)  do(k, k', cos 0) = C{[|A.? + |B.) F(p, 9) [4.2 + 41 B14, a} 
D pr ak! a (cos 0) . 


n 


3. - Numerical results when K nucleon interactions occur in S state only. 


In this section we restrict ourselves to spherically symmetrical energy 
independent elementary amplitudes a, and a, (6). With these assumptions the 
spin dependent scattering amplitudes B are zero and the spin independent 


(5) Experimental results on K+-Hydrogen scattering exhibit a nearly isotropic 
constant cross-section up to 200 MeV (?). 


oD 
n 
<© 
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amplitudes are given by: 


2M 1 
Aia x 3 gq?’ 
ne 
pers) 


where x is the fine structure constant and M the nucleon mass. 

We have calculated on IBM 704 digital computer the cross-sections for 
a,=— 0.24 fue as given by the K*-hydrogen experiment analysis (?). 
Several values of a, have been tried: a) = 0, 0.05, 0.06, 0.13 #/uc. Various 
incident energies have been considered: %k= 1, 2, 2.5, 3, 3.5u/fic, where mw in 
the x-meson mass. For the three processes involved we have computed both 
the differential and the total cross-sections restricted to cos 0 < 0.9 because 
of the Coulomb effect. In the case of inelastic and exchange scattering the 
angular distribution of K* mesons is obtained by integrating over k’. 

“2 The elastic and inelastic differen- 
do tial cross-sections are strongly domi- 
d(cos @) nated by two characteristics: 


— the Coulomb effects which 
are comparatively more important 
than in the zx-deuteron scattering 
because the K-nucleon forces are 
about ten times smaller than the 
r-nucleon forces; 


— the form factor F(q), F(p, q), 
G(p, 9); for 0— x the form factors 
considerably reduce the cross sections 
and the counting rate is too low. 


It appears very difficult to deduce 
useful information on elementary 
amplitudes from curves such as 
those given in Fig. 2. 


The exchange differential cross- 

section is more interesting because 

Coulomb effects are not present. The Pauli exclusion principle reduces the for- 
ward exchange scattering because we have assumed spin independent elementary 
amplitudes so that the final state of the two protons system would be a 
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triplet state. Consequently the angular distribution exhibits a maximum from 
which it is, in principle, possible to deduce the value of a. However, the 


do q=0.5 
ee fin) 
d(cos@) K-35 


cos@ 


lia St 


sensitivity of this maximum cross-section on the K-meson energy shows that 
a very good energy definition is required, which is not the case at the present 
time; Fig. 3 gives exchange. diffe- 

rential cross-sections for various 03 ne 
values of the energy and of a. 

In Fig. 4, we represent the va- 
riation of the total cross-sections 
with the energy. We note that 
the inelastic and exchange cross- 
sections are rather flat between 100 
and 200 MeV; but the ratio 0,,/0, ., 
increases with the energy and tends 
to an asymptotic value which can 
be computed by merely adding up 
the free cross-sections of K* mesons 
on the two nucleons of the deuterium 
(Coulomb repulsion is negligible in 
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the high energy limit): 


Tex (4, ic Ay)? 
Cc» TO) 239 
On ex 4ai ae (4, ca Mo) 


which permits an unambiguous determination of a) since it is known from 


experiment that |da)|<|a, 
Fig. 5 shows the energy dependence of the actual ratio; the asymptotic 
value is reached more rapidly by the differential ratio for 6 = x, but the expe- 


rimental uncertainty would be much larger because of the very low counting 


rate. We feel that the deduction of the above asymptotic limit from the actual 
ratio is not too dependent on the approximations made. On the other hand 
it can be seen that this measurement does not require very good energy de- 
finition. 


4. — Conclusion. 


The present calculation gives encouraging results about the possibility of 
a theoretical interpretation of a K* deuteron scattering experiment. In par- 
ticular we have shown that the exchange cross-section permits to draw con- 
clusions on the /= 0 elementary scattering amplitude. However this preli- 
minary work is unsatisfactory in two respects. First we should introduce pheno- 
menologically elementary P amplitudes which are certainly present in the high 
energy domain. Secondly it seems necessary to take into account the inter- 
action of the two nucleons in the final state which might enhance the cross- 
sections by an appreciable factor. A more careful calculation will be soon 
undertaken. 


APPENDIX 


The elastic form factor F(q) is related to the deuteron radial wave function by : 
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For the particular choice of the Hulthén function: 


yo(r) == {exp [— Br] exp[— rl}, 


we obtain an analytic expression for F(q): 


F(q) _14xB(x + 8) 16q(B — «)?(B + a) 


> arctg da si 
“ati 


In the inelastic case, we define the Fourier transform of the wave function: 


1 il 
k® + p? k2 + a2|" 


Polk) = [exp [tk - r| Wp ( Nol (= | 


We introduce the angle between the two vectors p and q. By definition: 


+1 +1 
Blp, 0) = | 2 (R:) d (cos 9) = | 3 (k5) d(c08 On) - 
-1 -1 


G(p, 9) = | AIA D 


The explicit expressions with the Hulthén function are: 


9 9 
FO D= Do +00 6-00) | PH +02) EH 0-07) 
4, +L + w/21)(6* + Lp — 027) 
pa(B? — &) ° (a? +[( ie }?) (62 +E + 0/29)” 
ne 2031048) 1 Le ape) a 
Pr PT palo + Be + (p? + @)/2) Le + (p? + @)/4 © +[(p — )/21° 
1 Pi cli crCUEARI 
pepe) °8 pa “i qg)2PJ : 


The normalization constant C is given by: 


aB(a + B)1 M 


Coi (2B)? 2.M.° 
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where M is the nucleon mass and M, the reduced mass of the K*-meson- 
nucleon system. 


RIASSUNTO (*) 


Abbiamo studiato gli scattering elastico, anelastico e di scambio K+-deuterio ser- 
vendoci dell’approssimazione impulsiva. In questo calcolo preliminare abbiamo trascu- 
rato l’interazione tra i nucleoni nello stato finale. Si sono effettuati calcoli numerici 
su macchina elettronica 704 [BM nel caso particolare che intervengano solo le ampiezze 
S del K+-nucleone. Si trova che la distribuzione angolare dello scattering di scambio 
e il rapporto delle sezioni di scambio efficaci alle sezioni efficaci senza scambio sono 
le sole grandezze che possano permettere di determinare facilmente la lunghezza di 
scattering nello stato di spin isotopico I=0 per le reazioni K+-neutrone, grandezza. 
attualmente mal nota, mentre lo scattering K+-idrogeno ha già permesso di valutare 
l’importanza dello spin isotopico 7 — 1. 


(*) Traduzione a cura della Redazione. 
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On the Elastic p-d Scattering at High Energies (’). 


L. FAVELLA and M. OLIVETTI 


Istituto di Fisica dell Universita - Torino 
Istituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto 18 Dicembre 1958) 


Summary. — We have investigated the validity of the impulse approx- 
imation in the elastic p-d scattering for protons of 310 MeV. The numer- 
ical results for the differential cross-sections and polarizations as funetions 
of the angle of scattering show clearly that this approximation is inade- 
quate for the prediction of the additivity of these observables as it has 
been proved experimentally. The reason of the failure is due to the 
large nuclear two-body interaction in the triplet states. 


Introduction. 


The scattering of high energy nucleons by complex nuclei raises the question 
whether the experimental values of cross-sections and polarizations can be 
explained by means of the data corresponding to the scattering of two free 
nucleons. 

The particular case of nucleon-deuteron scattering serves for deriving infor- 
mation about the continuum states of the neutron-proton system which are 
at present time not well known, because of the difficulty of producing a beam 
of monochromatic neutrons. This information is necessary to the scope 
mentioned before. 

There are measurements (!) for some types of quasi-elastic scattering which 


(*) This work is a continuation of a thesis in Theoretical Physics of M. OLIVETTI 
(Turin, Dec. 1957). 

(1) O. CHAMBERLAIN, E. SEGRÈ, R. D. Tripp, C. WieGanD and T. J. YPSILANTIS: 
Phys. Rev., 105, 288 (1957). 


D 


~ 


680 L. FAVELLA and M. OLIVETTI 


clearly show additivity of the differential cross-sections and polarizations of 
nucleons with energies greater than 300 MeV. 

At such energies one usually emploies the so called impulse approximation, 
whose validity can be checked by investigating whether it is able to predict the 
experimentally observed additivity. For this reason we have calculated the 
values of the p-d cross-sections and polarization using the impulse approxi- 
mation in the scattering of 315 MeV protons, which is known experimentally. 

For the effective calculation one needs of course to know the scattering 
matrix of two free nucleons. Since none of the various types ot phenomeno- 
logical or meson theoretical interactions is able to predict, with a good ac- 
curacy, the continuum states for the two-body problem, we have used for the 
proton-proton scattering a set of phase shifts calculated by STAPP (?) directly 
from experimental data. Since the data corresponding to neutron-proton 
scattering are lacking, we cannot evaluate what is left by using the impulse 
approximation besides the additive terms, so that we have evaluated only 
those terms corresponding to p-p scattering. 

However, as it will appear from a more detailed examination of the re- 
mainders, the corrections due to p-n scattering do not improve the situation. 
These remainders are functions of the scattering angle and are important al- 
ready at 15°. The situation is less favourable for the differential cross-section 
for which at 15° the additivity fails to be predicted with an error of nearly 
14% and this error increases to nearly 22% at 90°, while for the polarization 
it is less than 10% between 15° and 45° and increases rapidily for larger angles. 
Their values are shown in Fig. 2 and 3. 

For what concernes the total cross-section, it is easy to show (see Appendix) 
that there is a rigorous additivity if use is made of the well known theorem 
which gives the total cross-section as a function of the immaginary part of the 
scattering amplitude in the forward direction. 

In Sect. 1 we express the scattering amplitude operator of the three- 
body problem in the total spin space as a sum of the corresponding ope- 
rators of the two-body problem. This decomposition, supplied by the impulse 
approximation is used in Sect. 2 for evaluating the differential cross-section 
and the polarization in the p-d collision in terms of the two-body cross-sections 
and polarizations. 

The scattering amplitude operator for the incoming proton #(123) in the 
three body center of mass system is given by (#4) 


(1) H(125)= P'(12,3) A-LE"(1372)0 


(*) H. P. Srapp, T. J. YPSILANTIS and N. MErROPOLIS: Phys. Rev., 105, 302 (1957). 
(3) G. F. CHEW: Phys. Rev., 80, 196 (1950). 
(4) M. VERDE: Handb. d. Phys., 39, 143 (1957). 
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F' describes the scattering of particle 1 by 2 considering particle 3 as free; 
F" the scattering of particle 1 by 3 while nucleon 2 is free. #' and F” contain 
an average over the relative momenta of 
particles 2 and 3.-As we shall see later, 
the matrix elements of the operators F’ and 
fF” are linear combinations of those for the 
two-body scattering amplitude respectively. € 
We have calculated the differential cross- 

section o 


1 incoming proton 


; : 2-3 deuteron 
va and the polarization P,,, follow- 3 2 neutron 


ing the well known formalism of the density 


3 Biome 
matrix 


“ e EN 
(2) Opa = ae = ) LI oi Tr (SERE al PURE Ad IPS TRAI 7 


(3) GP = Tr [(1) PH] = Tr [o(1)(P'PT+ PP + PP + PP). 


S| 


The o(1) = (0,(1), o,(1), 0,(1)) in (3) are the usual Pauli spin matrices of 
the incoming proton. The factor } is the statistical weight of the six possible 
spin states (3 deuteron spin states x 2 proton spin states) which are all equally 
probable, since we assume no polarization in the initial states. 

Neglecting the mixed products of F’ and F”, as it is requested by our ap- 
proximation, (2) and (3) become 


Tr (EP FF"), 


(2) Opa = 


CHR 


ar Pr Tr{o(1) CE'F"++ P'P'H]} 
. n= PRA PET 


whence we shall deduce later 


16 
(2") Opa = 9 S(k, 0) fol AU ne Lee Ep} ’ 
(3) Opa Py wa On RA sa Onn E nas Non ? 


S(k, 0) being the deuteron’s form factor which is a funetion of the scattering 
angle and energy in the center of mass system of the three bodies. 

The remainders are all contained in the terms €,+e, for the scattering 
cross-section and in n, +, for the polarization. 


pn 


pon 44 - Il Nuovo Cimento. 
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1. — Reduction of the three-body scattering amplitude to the two-body one. 


Let 7” (J= 3,4) be a complete orthonormal set of spin eigenfunctions in 
the special case when particles 2 and 3 form a deuteron. They are quite deter- 
mined by the condition y*" 7” —0,,0,, and besides by imposing them to 
be simultaneous eigenfunctions of the operator $|o(1)-+e(2)+6(3)| with eigen- 
values 3 or } and of 3|6(2)+0(3)| with eigenvalue 1. Such eigenfunctions are 
obtained DE composition of the three triplet states of deuteron with those 
of particle 1. 


= = > x C(13J; um M) 2 (23) ay (1) - 


p=71 m=—# 


In order to reduce the three-body scattering problem to p-p and p-n scat- 
tering, it will be useful to write our eigenfunctions by means of those of the 
two bodies whose scattering is considered, for instance the bodies 1 and 2 
The y; are so given by 


j= A m=-1 


(4) gf = os > ALOG4I; M—m m M)" (12)a" (3) = > bi. (12)& (3) = 7, 
=1 
=1 


94) — he (3) T°(12) + pe) T1(19), 
3 V3 
J= > È 
= Ht = | Ormai, 
Gi xt = 47? = B(B)T- (2), 
3 (3) 3 
| | V=G = SS T2) A O1) sh (3) S(12) , 
o 
Li fg. (SOS à 3 
(A ge o ee WOR ae (12) -+ - B(3) S(12), 
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T"(12) and S(12) are the triplet and the singlet states of particles 1 and 2. 
We call: 


CRUE eZ 


k, and k, the momenta of the incoming and outgoing proton respectively; 
alto) the deuteron’s wave function normalized so that 


finto Pa al. 


y (123) the rigorous solution of the wave equation for the three-body problem 
which evolves from the initial state 


(5) exp[ik,-r]g;(0)x/(123) . 


Let V(12) and V(13) be the nuclear potentials which operate on space and 
spin co-ordinates, M the mass of the nucleon, P(13) the operator which ex- 
changes the space and the spin co-ordinates of particle 1 with those of 3. 

The rigorous expression of the scattering amplitude matrix of the incoming 
proton, in the three-body center of mass system is 


— 
DI 
ee, 


(Aa 
OES age (3 Dì 


[z°028) exp [— ik,-rlg*(o) [V(12) + V(13)](1 — P(13)) y/(123) do der. 


We shall now derive the expression of the operator F' in (1) as a function 
of the p-n scattering amplitude f(12). 
Let e"(12) (m=1...4} be the spin eigenfunctions of bodies 1 and 2 


(7) ge TA) et 92), EST (12) ei = (12). 


Let y7(12) be a complete set of rigorous solutions of the two-body wave 
a 
equation, corresponding to initial states labelled by the momentum p; and 
the spin state m. The rigorous expression of fi,(12) is thus 


eae | 2 8 eee 
(8) fom(12) == 5 fe" (12) exp [— ipa] V(12) yp, (12) dr, 
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where x=r,—r, and y, (12) evolves from the initial state 


exp (ip; x]2"(12) 


Let us define 
g(K) = [ex [iK-p]p(e) do , 


where K is the momentum of the third particle relative to the second. With 
the impulse approximation we expand w’(123) in terms of eigenfunctions Vir”) 
of bodies 1 and 2. We neglect during the scattering process the interaction 
with particle 3, multiplying by a plane wave travelling in direction 
w= ((r,+r;)/2) —r; and whose momentum is k, = K—(k,/2), so that the 
seattering amplitude operator depends upon the initial state, with momen- 
tum k,. 

w,x() is a function of the momentum k; of the incident nucleon and of K. 
The last che may assume all possible values according to the momentum distri- 
bution in the deuteron. Therefore 


(9) NL Is) ) Pr; (©) On1¢"(3) exp[ik, u|d'EÆ, 


m=1,2 
l=1.. 


BIS 


where «”(3) are the spin eigenfunctions of the third particle. The coefficients 
a, are defined by the condition that y'(123) approaches the initial state for 


r— 00, i.e. 


(9°) p”(123), — fan) exp (AK + wk;)(r2—1,) + 


5 k,; r+r 
i(K a rs; — 9 ") do dr dX - alk DI Dent (3)2' (12)] . 


4 M = 1,2 
Verdin’ 
One finds then 
A=3, u=i, 
therefore 
(10) Wir (0), Xp [19K -x+ 13k, -x]2(12) . 
Assuming as new variables x and p we have u = — (x/2)+e; r= x — (0/2) 
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and in account of (4) 


m,r=1,2 


1,8=1...4 


7 il M 
(IRD) Jes = SA, fed S*(k, 0) ea: Derr Ore À 
FIT \e ( 


-fexp [— tk,-x]2°(12) V (a) wr o(0) exp|ik;: la 3 


where 
rene 
Stk, 0) = fs" ( — D E La K | g(K) d*K ’ 
È | ZS | 


is the deuteron’s form factor and 0 the scattering angle in the three-body 
center of mass system. 

In eq. (11) we have taken into account the fact that i gives a contri- 
bution in the integral over X only in the neighbourhood of X= 0, and we 
have neglected the integral containing the exchange of coordinates of nu- 
cleon 1 with those of 3, which does not contribute in the forward direction 
and is unimportant at large energies and at small angles. 

From 


(e 


Dr, 0 go EXP [es k; ‘ale! 
it follows 


Pas > exp (t¢ h,-x] = \exp[ik;-a]a 


T—> © 


and in account of (7) 


We have now to evaluate the coefficients DT 


From (9’) it is clear that a/,=b/,, so that 


Ts” one 


Sia’, = Do, | 
ue T 


Since 


and 


we obtain 


(13) 2 bb; t= <4 123) |2°(12)> | <2" (12) | 771238)» . 
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The corresponding numerical values are thus derived from Eqs. (4) and (6°). 
The scattering process of 1 by 3, is evaluated in an analogous manner. 
The wave function y(123) becomes 


n 


(14) yp’ (123) = > o — r)bme"(3) exp [thy -u'] , 


where 


ey 


Li 


The term of the scattering amplitude which does not contain the permu- 
tation of 1 and 3 is reducible, as in the preceding case, to an integral of type (12), 
while the term containing P(13) gives the following integral 


m,r=1,2 


s,l=1...4 


4 DI Let 3 
Dy == L bibtudra 7 7a) [exp | Re) VO gia, 


using the new variables 


ALARE : : 
Yx,0(13) is the function which evolves from exp [i$k;, -x|2’(13), as follows from 
the condition at infinity for y”. 

The sum of the two integrals gives 


(15) F,(123) = 4 S*(k, 0) 


r 


rs rl 


DIANO 


i 


Î,1(13) are the matrix elements for the scattering amplitude of two nucleons, 
no matter if they are identical or not, whose initial and final momenta Pi and 
P; are related to k; and k, by (1) 


li, — k = aa 
(16) f= Pp Pr ; 


og ep) Pi- Pre 
2. — Evaluation of the differential cross section and polarization. 
We can now proceed to calculate the proton-deuteron cross-section and 


polarization as a function of those for p-p and p-n processes. 
The two-body scattering amplitude in order to satisfy the usual condition 
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of invariance by transformations of Galilei and space an time reflections, must 
have in the spin space the following expression (1) 


(17) Î = BS+# C(0,) ae Gu) "OF GT (Se ‘q)(0%) ‘q) + (6a) HAS VEDA ce 
yee [ (Ga) ‘q)(0 |G) — (Ga) P) (Ge -P)]T + N(6q)‘1)(6q):n)T, 


where : 
ee Lh del Paes De Pi Pe 
|PiAP.| | Ps — Pil IP; + pil’ 


T and S are the triplet and singlet projection operators 
T = t(Gqa) G+ 3), S = 7(1 —oq 64%) . 


The complex coefficients B, C, G, H, N, are scalar functions of the scat- 
tering angle (p;, p;) — 0 and of energy. 

Ordering the matrix f according to the spin eigenfunctions 2"(12) and 
choosing and z as y axes the directions of P and n respectively (5) we have 


1(G—H cos0) —vV2(i0+4Hsin0) 4(4+H cos6)—N 0 

| 4/2 (10 —4H sin0) N+H così —V2(i0—+Hsin0) 0 
Soana | 3(G+H cos0)—N V2(iC+4H sin0) 3(G— H cos0) 0 
0 0 0 B 


The operator f preserves the total spin, so that the elements f,, and f,, 
with « #4 vanish identically in both p-p and p-n cases. In the p-n case that 
is a consequence of the fact that the term (0,) — Go) ‘n is absent; its presence 
would imply a charge-dependence of nuclear forces. The ceofficients B, ©, 
G, H, N are different in the two cases, for the necessary antisymmetrization 
required bu Pauli’s principle. The matrix F' of (1) and (11), and so F”, is 
connected to f(12) by the following linear relations 


(19) GA |P" |) = YATA{0G 33"; M'—m m MCG 4d; M—mm M): 


-<j M'— m|f(12)|7 M — m) 


as one can easily deduce from (4), (12) and (13). 


(5) With this choice we have |Tr (onfF*)|=|Tr (0, af F*)| since the polarization 
must be normal to the scattering plane; see L. WoLrensTEIN and J. ASHKIN: Phys. 
Rev., 85, 947 (1952). 
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It is clear from (17) that f is a sum of terms each one expressible as a scalar 


product of a spin tensor of well defined rank with a space tensor of the same 
rank. Introducing these tensor components O (L= 0, 1,2), we have 


2 L 
(20) jm |i \jm>y = > > BE<im' | Or\im>. 


L=0 u=-£ 


From (17) one derives 


| BERRY, (Op), B°,=3(G+N)Y3(09), 
(20) } BY = CY;{n)(—1}, 
| Baz LG + N)Y,"(q) + L(G — H) Y,"(P) +iN te) (— 1)", 
<00|©3|00) = «1010, 110) = 4/4x, 
] AT 
10|02|10) = V = 
20" (10 | @2| 105 NE 


{01021105 = 2 WE 


Y* are the spherical harmonics normalized to one on the sphere of unity radius 


a (1) L+ ft (2h Je Gh = n)! 1 Q L+W 
vi + 
5 Da TC I DI | (sin 0)" (3 A si (sin 0)?" exp[iug] . 
Making use of Wigner-Eckart's theorem (6) for every matrix element one 
obtains 


(21) <jm' |f|jm) = > BY C(jLj; mum'Kj|0,}|j), 


Lyt 


and hence one derives the well known expression for the differential cross- 
section and for the polarization in the two-body scattering (see (2’) and (3’)) 


Le : I il Hy et 
(22) (12) Gipi CAVIA dro i Ti (St) BiB iG ODI > 


(°) See for instance, Rose: Elementary Theory of Angular Momentum, p. 85. 
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Il 
(229) 0(12) P(12) — Tr (o,(1) ff) = 


: il (2j +1)! pe ee ue 
ioe 2 Sori Bf W(1jLj; jL/)- 


*{ Br; OOL'L;1p'u'-1)4-1 oo he L'L;1u'u"+ 1)} jlo117)< 151915 HOME 
where the W(3,9:7js; 9j") are the Racah coefficients ("). 

Similarly using Wigner-Eckart’s theorem for the matrix elements of F 
one obtains 


(23) <J'M'IFIJM) = X Di 


MITI 


C(JLI'; MuM'XJ'|E,|J), 


where D, 


The differential cross-section o'(123) which is referred to the scattering 


and £, are quantities corresponding to BY, and O,. 


process by one of the nucleons in the deuteron is thus 


! 1 (e nies 1 \ lt LL — ll 2J' 1 jo D 
CARI a METRI 
JI LI 5 ni 


By comparison of (23) with (19) one obtains 


Ob) SOC Va Ty 


LIT" 


= D AA Bide We (ed ta) (2j +1)? W(Lid'4; 7J)<G|O,|9>. 


Consequently 0'(123) becomes 


D (42 BORA Aa eae 10); 
(26) 0'(123) =32 DI Aj A; À Ay ELLI 


:(—1)°Bi; Bd W(LjJ';jJ)W(Lj'Tk;I)G|0x]}G'|O|j>F, 


QI'+1)27+1)- 


and separating the contribution of the terms with j 7j from those with 
a; 
(SM, 


2 Ls. MRI Ent 
ORI ! 9 9 2 inal 72 ñJ +: ; J\2 Ja 
(26°) o (123) Go J +1)(2d > W? (Lid 5: jJ)7 (45) (A; ) Shae 


(Babol Ong 2 + 


3 n 43 
a ve W (0133; 14) W(0034; 03) Ref BS, B%<1|O,|1><0|O,|0>*} 


{7) See for instance, Rose: Hlementary Theory of Angular Momentum, p. 108. 
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(We have made use here of the triangular inequalities for the Racah coeffi- 
cients and of the values of A.) 

The first term at the right hand side of (26') contains, as one can easily 
see, just the same four terms with the same tensor ranks appearing in o(12), 
with different coefficients. 

By comparison of (26’) with (22) and denoting by C%(12) the coefficients 
of*|j|09,|j>| in 0(12), one finds 


(27) 0‘(123) = o(12) +5 Re {By,, Bre<1|@,|1}X0] ©, | 0>*} — 
$F 68(12) | 0 |, | 0> |? — zy C2 (12) |<1| D P— 
— ¥ Oi(12) |<1] O,[1> |? = $07 (412) |<1] @, [1 |? - 


| 117 | 

We can therefore assert that, apart from the first term, the others. are all 
comprised between 9% and 33% of the corresponding ones in o(12). 

It is noteworthy here to call attention on the preponderance of the 
term — À ,O7(12) |<1|@,|1>|? which practically dominates the remainder. This 
means font the impulse approximation is inadequate as soon as strong two- 
body interactions are effective. The amplitude {1 ©, 1) refers indeed to the 
triplet substate. Proceeding in the same way, the polarization results to be 


x (27 +1} (29 +1). 
LP PEER 
We i À, DA CIEL) 


RIPRODOTTI OR ELA TES ET) MC TETE UT) 
WT; fd") W(L'jot'%; jd") hi | Ox HD" fa | Oz 1e) <9" Jor | 
*[C(1LL';-1uu—1)Bi+ CALL';1u u+1) Bi \(=1) Be 


1 
(28) 7 Tr (o,(1)F'F*) = 
0) 


With reference to (2”) and (3”), we give here the expression of the cross- 
section and the polarization for the two-body scattering, and the remainders 


€ 1 2 (9 119 1 | 1 | il 
oI BRS 2h OP + A te Nee GaN 
6 1 il 1 dial 
29 3 | [9 Waar | lo T lo 
( ) re LA 16 (|B |? |G |?) SE 181) dr 


5 1 
+ 5; Re (GN*) +7 Re[(N+ @)B"], 


2 Re (ON*), 


2) = 
30 L 
tai = Im Ga (¢ + B—3N)0* — à HN* sin ol 


È 


690 


i lee en 


ON THE ELASTIC p-d SCATTERING AT HIGH ENERGIES 691 


In the numerical calculations of €, and 7, and their percentage with 
respect to the additive terms, we have used the set No. 4 of phase shifts in 
p-p scattering, as already mentioned (?). The choice of set No. 4 is due to 
practical reasons, since a large num- 
ber of tables and curves given by 
Stapp, are referred to this set. 


SRE ee Ee < 
QP 10° 20° 30° 40 50° 60° 70° 80° 90% 


Hig. 2. Rot Se 


However, as it is shown in Stapp’s work, the sets No. 1, No. 2, No. 3, No. 6 
are also good for what concerns the values of the cross-section and polariza- 
tion. The percentage value of the remainders e, and 7, are plotted as 
functions of the scattering angle 0 in Fig. 2 and 3. 


We wish to thank Prof. M. VERDE for suggesting this problem and for his 
kind encouragement and many useful discussions. 


APPENDIX 


Additivity of the total cross sections. 


We shall prove the rigorous additivity of the total cross section 0°, by 
help of the well known theorem relating the imaginary part of the scattering 
amplitude in the forward direction, to the total cross section 


(31) odi = va (es AE, . 
Since 3 tg 0/2 = tg 0/2 and, when 0 = 0, p= }k as one can see from eq. (16), 
it is easy to deduce that 


to) t ) 
(32) Ose == Om + On 
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It is just sufficient to bear in mind the expansion in terms of associated 
Legendre polynomials (?) of the fim. It is then clear that, in the forward 
direction, fis, fs1, fray fee, for, fas Vanish, since 


fis = fa = > a'Pi(9); (= Sh =e al Oe fu = — hres = 25 


where at, b!, ¢? depend upon the phase shifts and the energy. Moreover, it 
can be seen directly from the matrix (18) itself, that the sum of the last 
four elements vanishes. Since, 


Art ;j 
tot 
EIN re ; 
ites 


D 


and so for off), it follows 


ot = 
Opa = 


2 fai(o 
pete = Im SO 


q-e.d. 


RASS SIENNE) 


Si è studiata la validità dell’approssimazione impulsiva, nel caso dell’urto elastico p-d 
a 310 MeV. I valori numerici delle sezioni d’urto differenziali e delle polarizzazioni, 
come funzioni dell’angolo d’urto, mostrano chiaramente che questa approssimazione 
non è in grado di prevedere, per queste osservabili, l’additività verificata sperimen - 
talmente. Tale insufficienza deve attribuirsi alla forte interazione nucleare tra due corpi 
negli stati di tripletto. 
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Direct Decay x°->e' +e 


SD ADRE UT 


CERN - Geneva (*) 
Stanford University - Stanford, Cal. 


(ricevuto il 9 Dicembre 1958) 


Summary. — The direct decay 7°-—-e*-+e- is discussed and analyzed as 
a probe of «pion size ». 


The decay process 7°+>e++e7 is of interest as a probe of pion «size» in 
the following sense. Think of the decay proceeding via two intermediate pho- 
tons (!) according to the interaction 


(1) La if ag) E(x)- B(2). 


oy 


The photons then convert as shown in Feynman diagram 1. 


=P, 
Le 
Te. YG > 
TA 
vf 
Ter iL 


Aside from known and calculable electromagnetic factors describing the 
conversion of the photons to an e*e~ pair, the ratio of the decay rates 


Rate (7° > et + e-) 


2. ear Rate (7° — > + y) 7 


(*) Ford Foundation Fellow on leave from Stanford University, Stanford, Cali- 


fornia, summer, 1958. 
(1) In Eq. (1) ¢, is the pion-field amplitude, and #:B the gauge invariant pseudo- 


scalar formed from the electric and magnetic field strengths Æ and B. 


a 
Le 
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depends only upon the « structure » of the interaction as summarized, to indi- 
cate a special example, by an invariant structure factor /T(x — y)?] in 


(10) Bly = 2 deep n(0) dey Loe — 9) Ble) By). 


This is because the emission of virtual intermediate photons of large four mo- 
mentum is inhibited by the structure factor /. On the other hand the real 
photons in the 7° yy decay satisfy the Einstein condition k,k“ = 0 and 
thus have no invariant «length » with which to probe the structure in (1). 

Indeed, in the limit of a point interaction (1), the Feynman amplitude for 
diagram I diverges logarithmically (?). It is anticipated, however, that the 
«size» of this vertex will be & nucleon Compton wavelength (1/M) since 
nucleon-antinucleon pair annihilation into photons is usually held to be the 
mechanism underlying (1). If the Feynman amplitude for 7° + ete- is cal- 
culated with (1), and convergence is achieved with a Bopp-Feynman-Podolsky 
regulator corresponding to a mass M, it is found directly that the branching 
Taio (2) Istria 107. 

In order to proceed further and study more generally the dependence of r 
on the structure of (1’), it is advantageous to turn to a spectral representation 
for this vertex. The use of spectral representations for the vertex operator 
has been discussed in a number of recent communications (3); we shall only 
outline its application here. . 

We are interested in the matrix element and hence form factor defined by 


(3) <p, p= |i,(0)|0> = (4E, &, ) (us Fl(p,+p_) Jy; | bats) ‘ 


F is given by a dispersion relation 
(cel 


(4) (Qe) = + | DAS Q=p.+p-. 


o? — Q? — te 


0 


(*) Interaction (1) is not renormalizable because of the appearance of field strengths 
as required by electromagnetic gauge invariance. 

(*) G. F. CHew, R. Karpius, 8. GASIOROWICZ and F. ZACHARIASEN: Phys. Rev., 
110, 265 (1958); Y. NAMBU: Nuovo Cimento, 9, 610 (1958); R. OEHME: Phys. Rev., 111, 
1430 (1958); M. GOLDBERGER and S. TREIMAN: Phys. Rev., 110, 1178 (1958); 111, 
354 (1958); P. FepERBUSH, M. GOLDBERGER and S. TREIMAN: Phys. Rev., 112, 642 
(1958); S. D. DrELL and F. ZACHARIASEN: Phys. Rev., 111, 1727 (1958). The notation 
of this reference is used here. 
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in order to evaluate Im F(c?) we turn to the absorptive amplitude 
(5) Ape Pl We a Ve HG)? ) | Op) = 
i == (44/2) DA (27) d*(p, —Q)<p, lu, 9,(0) |n=><n* |9_(0) | 0 


n 


and find, with neglect of all but the two-photon contribution to Y and with 
n 

use of Born approximation for the electron-photon Compton scattering am- 

plitude, 


ah 40 
6 n ee e | OO? 
(6) Im F(Q?) (n m, In mi G(Q?) , 


The structure factor G((?) is defined by the matrix element to a two-photon 
state (momenta g, and g and polarization e, and e,; e? is the completely 
antisymmetric tensor of rank four): 


(7) CAE 1)-(0) 10) a ((4/2) (8m, 2) ) EA an EGG JE 2)?) o 


Equation (6) in (3) and (4) leads to the branching ratio (2) 


1 fo m.\2| f da : “ mi |? 
(8) VS (2 | | - G(amz=) ln et | 
= | 
0 


a—-1— die me 


which depends only on the space-time structure of the vertex (1'). The factor 
a? in this ratio is due to the additional electromagnetic vertices, and a factor 
of the form (m,//) In E/m, is typical for high-energy electrodynamie ampli- 
tudes; in this case £ corresponds to the pion rest energy. 

We give a few numerical results in order to show the sensitivity of ratio r 
to pion size (4). Choosing for simplicity 


lrn We 
G(aemi).= 
È (DE = E 


we integrate (8) to (°) 


1 fam, \? 212 272) 
= (= 2 | im p? In SI È In =| | 
À Mm | Me 


me | | i 


(4) Evaluation of G in a dispersion theoretic study of 7°--y+y decay is discussed 
in M. L. GOLDBERGER and 8. B. TREIMAN: Nuovo Cimento, 9, 451 (1958). 

(5) In a perturbation calculation which assumes that the 7° decays via a closed 
nucleon loop, the ratio r is finite and proportional to (In M/m_)?. The approximations 
here are not identical to this perturbation approximation since only the 2y-intermediate 
state is included in the absorption amplitude. In order to reproduce the perturbation 
calculation it is necessary to include in addition the two higher mass states (NN) 


and (NNy) where N(N) represent a (anti-)nucleon. 


(ox 
© 
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or r= 2,2:10-%, 1.2-107; and 0.3-107 for p=2M/m_, Mjm_., ard.i, res: 
pectively. 

It is not unreasonable to assume the possibility of detecting such a small 
branching ratio (9). We remark only that it is necessary to observe the decay 
2°->et+e with high resolution in order to reduce background from Dalitz 
pairs (7), 7°->y+et+e- which occur with the very much larger branching 
ratio 0.0118. The three-body decay into a Dalitz pair leads to a continuous 
distribution of final e- et pair energies and angles which is strongly peaked 
for the pair to emerge with small relative angle and with total energy & 5 m_ 
so that the one virtual intermediate photon is close to the mass shell. How- 
ever, a small fraction of Dalitz pair decays will yield final photons of very 
small energy so that the energetic e*e- pair will emerge essentially back-to- 
back. Such events will look like direct decays 7° ete. One can calcu- 
late directly from the papers of DALITZ (7) and of KRroLL and WADA (5) that 
the branching ratio for Dalitz pairs which leave an energy of up to dm_, with 
6<1, on the emitted photons is r, = (8x/37)ò'. For a resolution of 6% (de- 
ined by O— 0.06, or the photon energy ii S:4\MeV), 7 == 06-10) 7 ands 1s 
less than—or of the order of—the expected values of the branching ratio r 
for direct ete- decay. 

In this connection we note that detection of an accurate spectrum of Dalitz 
pairs may also provide information on the structure of the vertex 7° >y-+y. 
However, the direct decay ratio would appear to be cleaner from an experi- 
mental viewpoint and to be a more sensitive probe of high momenta and of 
small distances because of the momentum integration corresponding to the 
closed loop in Fig. 1. 

In resumé then, the branching ratio, 7, for m—+et+t+e-, is a measure of 
pion structure as operating at the vertex of Fig. 1. (We assume validity of 
quantum electrodynamics.) Data on this are extremely welcome, especially in 
view of the present unsatisfactory state of our understanding of elementary 
particle structure (°). It is also clear that observation of a larger branching 


(°) This is a much smaller branching ratio than discussed in P. LINDENFELD. 
A. Sacus and J. STEINBERGER: Phys. Rev., 89, 531 (1953). 

(7) R. H. Darrrz: Proc. Phys. Soc. (London), A 64, 667 (1951). 
(8) N. M. Krozz and W. Wapa: Phys. Rev., 98, 1355 (1955). 
(°) This has been particularly emphasized recently by the work of P. FEDER- 
BUSH, M. L. GOLDBERGER and 8. TREIMAN: Phys. Rev. 112, 642 (1958). See also: 
J. BERNSTEIN and M. L. GOLDBERGER: Rev. Mod. Phys., 30, 465 (1958); S. D. DRELL: 
VIII Annual International Conference on High-Energy Physics (CERN, 1958, to be 
published). It has been discussed here that pion structure may be a clue to a better 
understanding of nucleon structure. The «pion structure » operating in the nueleon 
structure problem corresponds to pion scattering with emission or absorption of a 
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ratio than predicted above could be interpreted as evidence for à direct non- 
charge exchange interaction of the form (pion) (lepton) (lepton). 


I would like to thank Dr. PETER HILLMAN for first calling my attention 
to this problem and for subsequent valuable discussion. I would also like to 
thank CERN, its Director, Professor C. J. BAKKER, and Professor B. FERRETTI, 
Director of the Theory Group, for their hospitality and for the award of the 
Ford Foundation Fellowship which made possible this visit. 


virtual photon and is different from vertex (1’). A future complete analysis of (1’) would 
relate its structure to «fundamental» ones such as the nucleon structure probed by 
virtual Compton scattering (4). 


RIASSUNTO (*) 


Si discute e si analizza il decadimento diretto 7°-e*+e- quale mezzo per saggiare 
la «mole del pione ». 


(*) Traduzione a cura della Redazione. 


45 - Il Nuovo Cimento. 
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Introduction of a Neutral Pseudoscalar Field 
and a Possible Connection between Strangeness and Parity. 


K. NISHIJIMA (*) 


The Institute for Advanced Study - Princeton, N.J. 


(ricevuto il 9 Dicembre 1958) 


Summary. — A neutral pseudoscalar field with vanishing rest mass is 
introduced. We first ask for the condition guaranteeing the vanishing 
rest mass, which turns out to be an invariance requirement under a 
certain kind of gauge transformation. Then assigning an appropriate 
transformation to each field one can restrict the possible forms of weak 
interactions. The main results obtained in this paper are as follows: (1) the 
change of strangeness necessarily requires the violation of both parity 
and charge conjugation; (2) the change of the number of electrons or 
of muons necessarily violates parity and charge conjugation; (3) stran- 
geness can change only by unity in weak interactions. 


1. — Introduction. 


The 7-0 puzzle led LEE and YANG to doubt of the validity of the conser- 
vation of parity (7) and various experiments that followed confirmed directly 
the violation of parity in weak interactions (*4). The violation of parity was 
then attributed to the two-component neutrino (5-7) but the original 7-0 


() Now at Osaka City University, Osaka, Japan. 

(1) T. D. LEE and C. N. YANG: Phys. Rev., 104, 254 (1956). 

(2) C. S. Wu, E. AMBLER, R. W. Haywarp, D. D. HopPes and R. P. Hupson: 
Phys. Rev., 105, 1413 (1957). 
. L. GARWIN, L. M. LEDERMAN and M. WEINRICH: Phys. Rev., 105, 1415 (1957). 
J. T. FrIED and V. L. TELEGDI: Phys. Rev., 105, 1618 (1957). 

A. SALAM: Nuovo Cimento, 5, 299 (1957). 
T. D. Lee and C. N. YANG: Phys. Rev., 105, 1671 (1957). 
L. LANDAU: Nucl. Phys., 8, 127 (1957). 
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puzzle was left behind the neutrino problem. It would, therefore, be worth- 
while to ask why the parity is violated also in neutrinoless weak interactions. 

In order to understand the nature of weak Fermi interactions several authors 
presented arguments in favor of the universal VA coupling such as two com- 
ponent theory (8), chirality invariance (*) and mass reversal invariance (1°). 
These arguments are, however, not fully satisfactory since, for instance, chi- 
rality invariance is valid in weak interactions only in the absence of strong 
interactions as well as rest mass terms, and also since mass reversal is not a 
quantum mechanical transformation. Consequently it seems likely that there 
is no room to accomodate a new invariance principle within the framework 
of the present theory. Hence the introduction of a new invariance principle 
will probably necessitate the introduction of a new degree of freedom into 
the theory, and we shall do so in what follows by introducing a new field. 
The introduction of a new field with vanishing rest mass seems to fit the 
present purpose, since the vanishing rest mass is generally guaranteed by an 
invariance principle, e.g. gauge invariance for the photon field and y;-invariance 
for the neutrino field. 

We have chosen a neutral pseudoscalar field denoted by B hereafter. The 
vanishing rest mass of this field is guaranteed by the invariance of the theory 
under the transformation B—-B-+/, where / is a real constant. In order to 
relate this transformation to other fields, we postulate that every Fermi field 
is simultaneously transformed as y» > exp[igd(y;-+a)|], where g and a are real 
constants chosen appropriately to each Fermion. We shall refer to the above 
transformation as « gauge transformation » hereafter. 

In Section 2 the general rule to write down a gauge invariant Lagrangian 
will be given, and also the relation of the gauge tansformation to the mass 
reversal will be discussed. In Section 3 the conservation law resulting from 
the gauge invariance will be given. In Section 4 we shall apply the theory 
to leptonic processes and determine the possible forms of Fermi interactions. 
In Section 5 the theory is extended to include strongly interacting particles. 
Assuming that all weak interactions are either of the Yukawa type or of the 
Fermi type and that every Fermi interaction comprises lepton wave functions, 
one can show that the strangeness selection rule AS = 0, +1 in weak inter- 
actions is an inevitable consequence of the theory and that the change of 
strangeness AS = + 1 necessarily requires the violation of parity and of charge 
conjugation. 


(8) R. P. Feynman and M. GELL-MANN: Phys. Rev., 109, 193 (1957). 

(°) E. C. G. SupARSHAN and R. E. MARSHAK: Phys. Rev., 109, 1860 (1958). 

(19) J. J. SAKURAI: Nuovo Cimento, 7, 646 (1958). See also J. Tromno: Nuovo 
Cimento, 1, 686 (1955); S. Hort and A. WAKASA: Nuovo Cimento, 6, 304 (1957); K. Fusit 
and K. Iwata: Progr. Theor. Phys., 19, 475 (1958); 20, 126 (1958). 
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2. Introduction of a neutral pseudoscalar field. 


We introduce a neutral pseudoscalar field B and ask for the condition 
guaranteeing the vanishing of the rest mass of the B quantum. The condition 
is given by the invariance of the theory under the « gauge transformation » 


(1) B—B+1, 


where À is a real constant. It is worthwhile to notice that the invariance re- 
quirement under the transformation (1) with an arbitrary real function (x) 
satisfying [JA(x) = 0 is superfluous, since there seems to be no theory that 
is invariant under this generalized transformation. Furthermore the gauge 


transformation (1) cannot be applied to charged fields. The reason for this 
is that this transformation is not compatible with the electromagnetic gauge 
transformation of the first kind. 

If the B field is coupled to other fields only through derivative couplings, 
the Lagrangian 


(2) he à 


OD \ wana E> 

A Jn OL 1 i 

is obviously gauge invariant. In order to derive useful results, however, it 
is necessary to assume that not only the B field but also other fields are simul- 
taneously transformed. As a possible gauge transformation of the first kind, 
we take 


(3) y>exp[ig4Ay;+ a)]y, 


where g and a are constants appropriately assigned to this spinor field, and À 
is the transformation parameter appearing in (1). It has to be noticed that 
the transformations (1) and (3) are given in the Heisenberg representation, 
and that the transformation of field operators in the interaction representation 
is quite different from the above ones. 

Then a possible gauge invariant Lagrangian is given by 


SE \2 
(4) L' = —Y(y,0,+ m exp[— 2igy;B]) le) è 
2 \0æ, 


A gauge invariant interaction is always obtained by the following substitution: 


(5) mpy > my exp [— 2igy;B]y . 
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This general method of constructing a gauge invariant interaction corresponds 
to Dirac’s substitution in quantum electrodynamics 


(6) VY OY > PY ple, — CAP. 


In addition to the above interaction one can add interactions of the derivative 
coupling type as appeared in (2). The gauge invariant interactions of this 
kind correspond to the gauge invariant Pauli type interaction in quantum electro- 
dynamics. Itis usually taken for granted that the Pauli interaction does not exist 
and that the electromagnetic interaction of the Dirac type (6) is the only one 
realized in nature, and this postulate is called the principle of minimal electro- 
magnetic interaction. Correspondingly we shall assume that the only inter- 
action of the B field with other fields is given by the substitution (5) in the rest 
mass term. Later in Section 5 we shall see that this postulate needs a minor 
modification concerning the pion-baryon interactions. This is indeed an es- 
sential assumption in determining the possible forms of weak interactions. 
There are two points of importance: 1) The electromagnetic interactions 
— tepy, pA, is invariant under «our» gauge transformation; 2) The inter- 
action of the B field with other fields in the absence of weak interactions con- 
serves parity and charge conjugation separately. The B field is odd under 
space reflection and even under charge conjugation. 

Next we shall discuss the relation of the gauge transformation to the so- 
called mass reversal. Under a special gauge transformation corresponding to 


(7) A= 1/29, 


the rest mass term is transformed as 


(8) my exp [— 2igy;B]y > — my exp [— 2igy;B]y , 


if only the B field is transformed. Since the B field appears only in the rest 
mass term this transformation can also be regarded as the transformation of 
the rest mass 


(9) M —— M. 


Under this discrete transformation à spinor field is transformed according 
to (3) as, e.g., 


yp > ty, explixa/2]y, 
(10) for a= +1. 
PT Ys » 
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The transformations (9) and (10) together form the so-called mass reversal 
transformation. As seen from (9) in which a c-number undergoes à trans- 
formation, mass reversal is by no means a quantum mechanical transformation, 
but we gave an equivalent quantum mechanical transformation in terms of 
a new field. 


3. — Conservation law. 


From the invariance of the Lagrangian (4) under the simultaneous trans- 
formations (1) and (3), one can derive a conservation law of the form 


od 
11 ia 
( ) OD 0 1 
where J, is given by 
OB bps 
(12) J, == On. = 19Yy » (Vs =F a)y È, 
IT, 


Hence the following quantity Q is conserved: 


(13) Q =| da. 


This conserved quantity Q corresponds to the electric charge in quantum 
electrodynamics. One can show that an operator of the form (12) generally 
leads to continuous eigenvalues. However, in our case so far as observable 
quantities are concerned @ has just discrete eigenvalues. This will be shown 
in Appendix A. 

The question to what extent the above conservation law restricts the pos- 
sible forms of the S matrix is not a simple subject, and we shall discuss it in 
a simple case in Appendix B. 


4. — Application to leptons. 


In this section we shall discuss the properties of leptons based upon the 
new invariance principle . The usefulness of this theory as applied to leptons 
rests fully upon the choice of a suitable set of transformations for leptons. 
Furthermore we shall assume the lepton conservation as well as the two compo- 
nent theory the of neutrino. 

When e-, 17 and two component v are particles, we choose the following 


N 
S 
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transformations among many other possibilities: 


y, > exp [19475 — 1)ly, , 


py Wa Py ? 
(14) 
and 


y. > XP [— igÀ(y5 — 1)]y, . 


Then the possible Fermi interactions invariant under the above transformations 
are given by 


1) 6-decay 

(15a) Gi Pp¥o(1 + ds) Pa VA + Ys)~y + he, 
2) p-capture 

(158) PA + 4275) Pa Pu%o(1 + ¥5)Py + h. 0. 


3) u-e decay 


(15¢) Gsp,Y A+ Ys) Pe Pol + vs), + h. c. , 


All Fermi interactions turn out to be of the VA combination just as other 
arguments give, and the u-e decay gives the maximum violation of parity in 
this unrenormalized form. However, unlike other theories we cannot put the 
ratio Gamow-Teller to Fermi a, or a, equal to unity, and further there is no 
relationship among the three coupling constants G,, G, and G;. The trans- 
formation (14) is so chosen as to give negative helicity to leptons and posi- 
tive helicity to antileptons in weak interactions (1:12), In this section baryon 
wave functions are treated as if they were not transformed, but we shall see 
in the next section that the interactions (15a) and (15d) are also invariant under 
the baryon transformations. 

When e-, u* and four-component v are particles as proposed by Kono- 
PINSKI and MAHMOUD (1%) several years ago, the vanishing of the rest mass 
of the neutrino is guaranteed by the invariance of the theory under the simul- 
taneous trausformation 


(16) PW, Voy » Yo > We eo Ye 


(11) H. FRAUENFELDER, R. BOBONE, E. VON GOELER, N. LEVINE, H. R. LEWIS, 
R. N. PEAcOCK, A. Rossi and G. DE PAsquartI: Phys. Rev., 106, 386 (1957). 

(12) M. GOLDHABER, L. GRODZINS and A. W. SUNYAR: Phys. Rev., 109, 1015 (1958). 

(13) E. J. Konopinsxi and H. M. Maxmoup: Phys. Rev., 92, 1045 (1953). 
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In this case we can again employ (14) provided that y, as well as y, always 
refer to the field operators for negative leptons. As has been discussed before 
one cannot distinguish between these two alternative choices so far as pre- 
sently existing experiments are concerned (14). 


5. - Application to strongly interacting particles. 


We shall now apply the requirement ot gauge invariance to strongly inter- 
acting particles as it was the first intention of this paper. Among many other 
possibilities we shall choose the following set of transformations for these 
particles in analogy to the case of leptons: 


1) baryons 
(EEE) y > exp [2 (C 1) aa vy , 
2) mesons 


(18) D > explifiU]®, 


where U is the sum of nucleon number N and strangeness S, and f is generally 
different from g for leptons. However, it seems likely that there is a simple 
relation between f and g. 

If the pion-baryon interaction is of the derivative coupling type, 2.e. 


AG 
(19) Ù n PTYVo%° Oc ; 


then this interaction is already invariant under the gauge transformation (17). 
If, on the contrary, it is of the form of pseudoscalar coupling, i.e. 


(20) IGPTY;4'® , 


then this interaction is no longer invariant, and in order to make it invariant 
we follow the pattern utilized in modifying the rest mass term by substituting 
the following interaction for (20): 


(21) iGpyst exp [— 2ify;B]y@ : 


(14) K. NISHIJIMA: Phys. Rev., 108, 907 (1957). 
() The factor (— 1)” is necessary in order to introduce leptonic d ecays of strange 
particles. 
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However, we need not modify the baryon-kaon interaction if it is strangeness 
conserving and of a non-derivative coupling type either scalar or pseudoscalar. 

The conservation law that follows from the gauge invariance is given again 
by (11) with J modified in the following way: 


(22) J, — JV nh : 
where 
tw oD oD 
23 J° = i U - lee a E 
GR rt DI OTT ar. gore t || 
0B 


Le x 
REATI 


baryons 


dx 


pe 


(24) Ji = 
Correspondingly the conserved quantity © can be split into two parts 


(25) G=9,+9, @=fU. 


In the above expressions the contributions from leptons are discarded since 
they are not essential in the following discussions. In strong interactions Q, 
and Q, are conserved separately since U=N-+S is conserved. Now let us 
emphasize that they have opposite transformation properties under space 
reflection and charge conjugation, 4.e. 


(26) POLE EREDE: 
and 
(27) CO Caan CT COO = TE 


In weak interactions only the sum Q is conserved since U is no longer con- 
served, and this means that quantities of opposite transformation properties 
are mixed up, or in other words, that both parity and charge conjugation are 
wiolated in weak interactions. To be more precise, we shall make further 
assumptions. Assumption 1): Weak interactions are either of the Yukawa 
type or of the Fermi type. Assumption 2): Every Fermi interaction comprises 
lepton wave functions, or in other words, there are no Fermi interactions con- 
sisting of baryon wave functions alone (*). 


(*) In this connection it might be suggestive to assume that lepton wave functions 
occur only in Fermi interactions but not in weak Yukava interactions. Indeed Yukawa 
interactions of the form Kyy, Key cannot be gauge invariant. 
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From these assumptions we can draw two important conclusions : 


1) In weak interactions the strangeness S obeys the selection rule 
AS= 0, +1 provided that we assign zero strangeness to leptons. 


2) When strangeness changes by unity, one can determine the unrenor- 
malized form of such an interaction uniquely. 


In order to prove the first statement we shall show that a bilinear form 
consisting of baryon wave functions and changing the strangeness by two 
cannot be gauge invariant. Under the transformation (17) y=0y, is trans- 
formed as 


(28) a0 ys > Hs exp if (— 7 +1) DI 


O exp in (- Ê 1) 


The invariance requires 


exp i (- mn + 1) O exp 


(5 +1) x 
but the only possible solution of this equation is given by 
(29) OUR 
This proves the selection rule AS= 0, +1 for Yukawa interactions of the 
form BBr and for Fermi interactions under the assumption 2). This can also 
be proved for Yukawa interactions of the form BBK. If the assumption 2) 
is dropped, AS = +2 cannot be excluded theoretically as given by the ex- 
ample An-An. 

The proof of the second statement is illustrated by considering 
(30) N > pErT : 
The bilinear form y,Oy, is invariant only when O assumes the form 


(31) O = AURAI 


which is the only solution of the equation requiring gauge invariance 


(ae. 


exp 


inf —2 ze i) 0 exp 
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Hence the unrenormalized form of the decay interaction for (30) is uniquely 
given by 


(32) Pool + 75) Pa “Oo” « 
Similarly for the process 


(33) ET = N+x., 
we get uniquely an invariant interaction of the form 
(34) Proll — Y5)Ya 0,9 - 


It is very suggestive that the signs of y; are different in A and & decays showing 
opposite polarizations for the decay products p and A. 

Finally it must be mentioned that when strangeness does not change as 
in 6 decay the form of interaction cannot be determined uniquely as illustrated 
by (15a) and (155). In this point the present theory differs from other theories 
basically. 

To conelude the author would like to remark that as far as the coupling 
constants f and g are sufficiently small the B quantum would hardly be ob- 
served. 


This work was initiated during the author’s stay at Brookhaven National 
Laboratory and completed at the Institute for Advanced Study. The author 
would like to express his sincere thanks to Professor S. A. GOUDSMIT and to 
Prefessor J. R. OPPENHEIMER for the hospitality. 


APPENDIX A 


Eigenvalue problem for the operator Q. 


From the Lagrangian (4) we see that Lis a periodic function of B, namely (4) 
is invariant under the transformation 


(A:1) B>B+alg. 


This periodicity requires that all physically observable quantities be invariant 
under this transformation. In order to determine the eigenvalues of Q, we 


D 


= 
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consider an eigenvalue equation of the form 
(A.2) ROSEN IAA 


If the eigenoperators Y form a complete set of operators, i.e. if the equation 
[0, X.]= 0 valid for every « requires that O be a c-number, the eigenvalues 
of Q are completely determined and given by 


(A.3) q= 04, (cy: integer), 


provided that QD,= 0, where ®, is the vacuum state. 


For the sake of definiteness, we take a=—1 in (12) 
0B RG 
(A.4) Ja = Ox + 190Y » (5 — 1)y ? 


then the eigenoperators and corresponding eigenvalues are given by Table I. 


a MATE 

Xx re do 
(V¥5—1)y 2g 
(75 + 1)w 0 
vil) — 29 
p*(Y5 + 1) 0 
exp [— 2ingB] K 2nq 
0B/0x, 0 


The invariance under (A.1) requires that » occurring in exp [— 2ingB] 
should be an integer so far as observable quantities are concerned, and hence 
we get all the eigenvalues of @ 


(A.5) q = (even number) g. 


For this reason it is sometimes convenient to introduce @/2g, since it has 
integer eigenvalues even after renormalization. In this connection it is desir- 
able to choose f equal to g or a simple multiple of g. 
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APPENDIx B 


A low energy limit thoerem for the emission or absorption of a B-quantum. 
Let us consider a gauge transformation U, which is defined by 


(B.1) UBI BE 
(B.2) U;yU;' = exp[igA(y; + a)]y. 


U, can be expressed in terms of Q by 

(B.3) U,= exp 101 

The statement that the theory is invariant under this gauge transformation 
is expressed by 

(B.4) UHU:=H, 

where H is the total Hamiltonian consisting of free part H, and interaction 
part H;,,. However, this does not always lead to 

(B.5) USS 

where S is the scattering matrix. The condition for (B.5) is more severe 
than (B.4) and is given by 

(B.6) CH =H, U;H UG ay 


(B.6) is true only for special choices of the value of 1, i.e. 

(B.7) Ne a (integer) . 

Let us denote U(z/g) simply by U, then we have (B.6) and also 
(B.8) UBI Be Tg = vin 


In general the transformation of incoming field operators is different from 
(B.1) and (B.2) and is very complicated. The simple transformation law (B.8) 
is valid only for special choices of the value of 4. 

The S matrix is a functional of B;, and for this special 2 we have 


NIT 


(B.9) U" S[ Bin] Cer => S Bin == LA = S[ Bin] e 


Then (B.9) gives us the following result: 
Let s(k,, ke, -.-1 Pis Pa, ---) be the S matrix element corresponding to a 


process in which B particles with four-momenta k,, k,, ... either incoming or 
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outgoing as well as other particles with four-momenta pi, Po, --. participate. 
Then (B.9) requires that s should vanish if any one of these momenta k,, ko, ..- 
is equal to zero. Thus we get a kind of low energy limit theorem. So long 
as S matrix elements on the energy shell are concerned we cannot expect any 
further relations. This result suggests that there is no difficulty corre- 
sponding to infrared divergences in quantum electrodynamics. However, the 
ultraviolet divergences get worse than in quantum electrodynamics as a general 
trend of weak interactions. 

In the present paper all results are referring to unrenormalized forms of 
interactions. However, it should be mentioned that the conservation law can 
explicitly be expressed in terms of renormalized operators alone using the 
transformation properties of field operators. 


RIAMSSIUNENONS) 


Si introduce un campo pseudoscalare neutro con massa a riposo tendente a zero. 
Cerchiamo dapprima la condizione che garantisce l'annullamento della massa a riposo, 
e si trova essere un’esigenza d’invarianza rispetto a una certa specie di trasformazione 
di gauge. Assegnando poi un’appropriata trasformazione a ogni campo si può restrin- 
gere il numero delle possibili forme dell’interazione debole. I principali risultati con- 
tenuti nel presente lavoro sono: 1) il cambiamento di stranezza richiede necessaria- 
mente la violazione della parità e della coniugazione della carica; 2) il cambiamento 
del numero di elettroni o di muoni viola necessariamente la parità e la coniugazione 
della carica; 3) nelle interazioni deboli la stranezza può cambiare solo di un’unità. 


(*) Traduzione a cura della Redazione. 
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On the Spin Dependence of the A-N Interaction (*)(**). 


L. E. Picasso and S. ROSATI 


Istituto di Fisica dell Università - Pisa 
Istituto Nazionale di Fisica Nucleare - Sezione di Pisa 


(ricevuto il 10 Dicembre 1958) 


Summary. — Rhe ratio wi between the lifetimes for the decays 5H, > 
>%He+n- and 8H,->D+p+n has been calculated as a function of 
the spin j of the 8H, (j=4 or j=3). The experimental value wp & 1 
favours j7=4, showing that the A-N potential is more attractive in the 
singlet than in the triplet state. This result does not depend critically 
on the choice of the wave functions, because wÎ/w? is independent of it. 


1. — Introduction. 


The available experimental results on the binding energies of the hyper- 
fragments, may be interpreted by assuming that the A-N interaction depends 
on the relative spin orientation of the two particles. The phenomenological 
analysis (1) on the binding energies of the lightest hypernuclear systems allows 
two possible solutions for the spin dependence, which lead to two different 
values for the hyperfragments spins. 

In this work it is shown how one can choose between these two possibilities 
by studying two decay processes of the *H, hyperfragment, the simplest one 
among the known hyperfragments. 


(*) This work has been presented at the XLIV Congress of the Italian Physical 
Society, 6-11 November 1958. 

(**) When this work was completed, we have received a preprint of a paper by 
M. Leon, who calculates the fraction of decays in two bodies of the 3H, and compares 
it with the experimental value; his conclusions are in agreement with the results we 


have here obtained. 
(1) R. H. Darirz and B. W. Downs: Phys. Rev., 111, 967 (1958). 
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The two decays we want to consider are: 


(1a) Ie DE OR ET, 
(1b) oie Hema 


Our purpose is to calculate the lifetimes 7) and 7, for the two channels as 
functions of the spin j of the hypertriton. 

It will be shown that the branching ratio w’= 1/1; depends strongly on 
the hypertriton spin j and on the ratio «= |g_|?/|/g, |? between the emission 
probabilities in s- and p-waves of the 77 meson in the free A decay. 

The experimental value (2?) w,,, < 1 singles out j= 3, showing that the 
A-N interaction is more attractive in the singlet than in the triplet state. 

This result is in good agreement with the conclusion reached by DALITZ 
in the study of the decay of the ‘H, and with the field-theoretic calcu- 
lations (5). 


2. — The interaction hamiltonian for the A decay may be written (*) (in 
the non-relativistic approximation) 


(2) H'= | dr y, (r) 


Di 


VI 
g-+ 9 PE valr) (1) , 


where 


y,(r), alr), g(r) are the field operators for the proton, lambda an x re- 


spectively 
© = (dz, Oy) 0.) are the three Pauli matrices, 
k, is the 7- momentum in the free A decay. 


If we neglect the influence of the Pauli principle in the final state for the 
process (1a), the lifetime 7? does not depend on the spin j of the *H a (indeed 
in this case the decay probability of the A does not depend on the orientation 
of its spin with respect to the deuteron spin). 


(*) For references to some *H, decay-events see, e.g., G. FRANZINETTI and G. Mor- 
PURGO: Suppl. Nuovo Cimeto, 2, 785 (1957). 

(3). R. Levi-SerTI, W. E. SLATER and V. L. TELEGDI: Suppl. Nuovo Cimento, 10, 
68 (1958). 

(4) D. B. LicareNnBERG and M. Ross: Phys. Rev., 103, 1131 (1956) and Phys. 
Rev., 109, 2163 (1958). 

(5) F. FERRARI and L. Fonpa: Nuovo Cimento, 5, 842 (1958). 

() We have everywhere used units #= mass (n-)= c= 1. 
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The square modulus of the matrix element for the transition (1a) summed 
on the final spin states and averaged with respect to the initial ones, is given by 


span 


A = Jan ar ar Polis — Pol) Pra, (ln Fm Ta) EXP | (ra 3 (Fat r))(&+p)]; 


| Hn 


= (+R) lo + la 


zh (O) 


where è 


r,, r, and r, are the neutron, proton and A co-ordinates, 


k and p are respectively the 7- and free proton momentum, 


Y, and ®,,, the space wave functions describing the deuteron and the 
hypertriton. 


The lifetime for the process (1b) on the other hand, has a very strong 
dependence on j, because for 7 = 3 only x in p-waves may be emitted in the 
decay. 

One obtains 


1 — ik È 
(3) ji Mali = do) 19 ge Ho È ||4 |, 
3 
(4) lio I A SG +R) 3] g; |? pid E ’ 
where 


A= = far, dr, dr, Pie (Ta, Tp) TA) Dorn, (a, Pp) TA) exp[— tk- ra], 


k2 
— = 1.28 
he i 
and Y,,, is the space wave function for the *He. 
To calculate the explicit values of the lifetimes we shall use: 


a) A Hulthén wave function to describe the final deuteron in the 
channel (1a) 


1 J 
Poe (exp[—ar]—exp[—br]),  Œ=ir—r,l), 
where 
(5) @ = 0.325, DI-1L:413t. 


e 46 - Il Nuovo Cimento, 


~ 
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b) Irving’s (9) wave function for *He 


Pr ex [Art (rat (rara; 


with À — 1.3. 

The hypertriton wave function is taken to be the product of a function ®,, 
describing the relative proton and neutron positions, and a function ©, 
describing the relative position of the A and the deuteron center of mass (7) 


(6) Dax, = D, (ra ral) B,t(ir, —àr, == 55) O 


Such a wave function disregards correlations among the positions of the 
three particles. 

The *H, is simply considered as consisting of a A bound to a nuclear core 
formed by a distorded deuteron. 

This may appear an oversimplification: indeed the correlations between 
the positions of the A and the nucleons are very important if one wants to 
predict the strength of the interaction A-N by using variational methods (*). 

In this case it is necessary to use a trial function for the hypertriton with 
sufficient flexibility in order to allow strong space correlations between the A 
and the nucleons. 

For the present problem, however, the use of a simple funetion of the 
type (6) is justified by the fact that, as will be shown, the results do not depend 
critically upon the choice of the wave functions. 

We give now to (6) the following expression 


— 
I 


1 
Dir, Ve (exp [— ar] — exp[— br]) exp[— yo], (eo =|ra—4(r.t+r,)|), 


where a and b are given by (5). 

The parameter y has been determined in the case that the deuteron is sup- 
posed undistorted, by a variational calculation for the A-N interaction strength 
necessary to give a binding energy for the *H, of 2.826 MeV. | 

If one takes à A-N interaction range which corresponds to an exchange 


(5) J. Irvine: Phil. Mag., 42, 338 (1951). 
(7) Various authors have used functions of this form for variational calculations 
on the *H, ; see, e.g.: G. H. DERRICK: Nuovo Cimento, 4, 565 (1956); J. T. JonEs and 


J. M. KELLER: Nuovo Cimento, 4, 1329 (1956); D. B. LICHTENBERG and M. Ross: Phys. 
Rev., 110, 737 (1958). 


(5) R. H. DaLitz and B. W. Downs: Phys. Rev., 110, 952 (1958). 
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of two x, one gets y= 0.75; for a range corresponding to an exchange of 
one K-meson, one would get y = 0.67. 
By using standard methods, one can easily get the lifetimes for the two 
processes (1a) and (1b) and therefore the values for w° and w*. 
The result is (for y = 0.75): 
0.85x + 0.407 , 0.850 + 0.407 
(3) wt = ——____--____; ee 
0.873x + 0.124 0.496 


where =|g-|?/|gs |?. 

The available experimental data are not sufficient to determine accu- 
rately x. 

They are compatible with 1 < x < 5 (1). 

In Table I, 3-rd and 4-th column, we have listed the values for w* and w* 
obtained from equations (8), for different values of x. 


TaBLE I. — The numerical results for R and for the w’s as functions of x. The 3" and 
the 4° column refer to the wave function (7); the 5° and the 6!" column refer to the wave 
function (9). 


| a R | we wi | | wi wi 
| 
SS DIE | ct _ 

Il | 2 1.26 2.53 | 1.75 3.5 
| 1.5 2.9 i Hz 3.4 1.62 4.7 

2 Salud 115 4.25 | 1.55 # 5.84 

3 5.56 1.08 | 6 i 1.48 8.2 
| 5 | 9 1.04 | 9.4 1.41 i 12.7 
| co co ; 0.97 fore) 1.32 co 

| 


We see that the experimental value w,,, < 1 greatly favours spin + for 


hypertriton. 
We may now ask how much this conclusion depends upon the choice for 


the wave functions which describe the states of *H,, *He and deuteron. 
To answer this question we first study the dependence of w! and w? upon 
the parameters of the wave function (7). 


is 
beam 
L= 
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In Fig. 1 we have plotted w* and w against y for 7-1 and #=2 with 
a and b given by equation (5). We see from the curves that within a fairly 
large range of y around 0.7 the situation is practically unchanged. 

In the same way, we could analyse the dependence of the w’s upon a and b: 
it can be shown (°) that in this 
case too, the results do not depend 
critically upon the choice of these 
parameters. 

We still would like to know 
how the above results would be 
modified by a different choice of 
the initial wave function. We may 
in general expect that the values 
for wt and w? would in this case 
be different. 

However we can easily see 
that the ratio 


2 
A0? 


Fig. 1. — wi and w? obtained from (7) re 


E) 
we 


against y. 


is independent of 7, #, and Da 

Indeed, as the lifetime for channel (1a) does not depend on the spin 7, 
one gets: 

4 (HAP 9 1 

ipa FA Ue da Ny i 


where | H,’~?|? and |H,/-?|? are given by (3) and (4) respectively. 

KR depends only on « and when x varies between 1 and 5, R varies cor- 
respondingly between 2 and 9. This shows that even if the wave functions 
we have used are incorrect, the conclusion on the spin is very probably still 
true. Indeed to get a spin 3 one would need to modify w? and w° by a factor 
of the order 1/R (in the more pessimistic case R — 2), which seems rather dif- 
ficult to believe. 

To verify this we have calculated w? and w? using a different wave function 
for ®H-,enamel;#(59) 


1 1 
(9) Din, — = (exp[— er] —exp[—dr]) > (exp[— ao] — exp[— Bol). 

(°) For more details on this question see the report of our work, to be published 
on Rendiconti del XLIV Congresso della Società Italiana di Fisica. 

(9) D. B. LicaTtENBERG and M. Ross: Phys. Rev., 110, 737 (1958). We are very 
much indebted to the Authors for communicating the values of the parameters which 
they used. 
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The values for w? and w° one gets in this case are given in Table I, 5-th 
and 6-th column. We see that in this case too w? is much more close to the 
experimental value of the branching ratio. 

In general we see from Table I that the agreement improves with increas- 
ing x, which could perhaps be taken as an indication that x is indeed fairly 
larger than 1. However this cannot be considered as a definite result due to 
the approximations involved in the calculations and the uncertainty Of w,,,- 

In conclusion we may say that our results are in favour of spin 4 for *H, 
and therefore for a stronger attraction in the singlet than in the triplet state 
of the A-N system. 


The authors wish to express their thanks to Prof. L. A. RADICATI for his 
constant interest in this work and for many helpful and stimulating discus- 
sions, and to Prof. D. B. LICHTENBERG for a valuable conversation. 


RIASSUNTO 


Si calcola in funzione dello spin j dell’*H, (j = + o j = à) il rapporto w’ fra le 
vite medie dei due decadimenti 3H, +?He + x e 8H, >D + p-+7-. Il miglior accordo 
col valore sperimentale w,,, < 1 si ha per j = 3, e ciò mostra che il potenziale A-N 
è maggiormente attivo nello stato di singoletto che in quello di tripletto. Questo risul- 
tato non dipende criticamente dalla scelta delle funzioni, poichè w?/w? ne è indipendente. 


NOTE TECNICHE 


Camera a bolle a propano. 


A. Loria, P. MITTNER, I. SCOTONI e G. ZAGO 


Istituto di Fisica dell Università - Padova 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 4 Dicembre 1958) 


Riassunto. — Si descrive una camera a bolle a propano di circa due litri 
di volume e si indicano scopo ed alcune modalità dell’impiego fattone 
al sinerociclotrone da 600 MeV del CERN. Il meccanismo di espansione, 
a membrana, permette una ricompressione rapida: la proporzione di 
tempo sensibile può quindi essere notevole. Sistema di riscaldamento 
assai pronto e flessibile e particolare struttura delle finestre determinano, 
senza inconvenienti di altra natura, buone caratteristiche termiche. Un 
condensatore di luce conferisce maggiore efficienza al sistema usato di 
illuminazione a campo oscuro. 


1. — Introduzione. 


La struttura generale dello strumento, in vetro e metallo con guarnizioni, 
ha costituito il naturale sviluppo di precedente attività di questo Laboratorio 
in materia (*). Sia per renderne praticamente certo il regolare funzionamento 
ed eliminare quasi del tutto il rischio legato all’impiego di fluido infiammabile, 
sia per ridurre le distorsioni così reali come ottiche (2) ed ottenere un’illumi- 
nazione più intensa a parità di energia dissipata nel flash, sono stati però 
impiegati alcuni nuovi accorgimenti la cui conoscenza si ritiene possa essere 
utile per chi si trovi a dover impiegare camere a bolle, particolarmente a pro- 
pano, in analoghe circostanze; essi risultano, coi vantaggi che ne derivano 


(1) P. Bassi, A. Loria, J. A. MEYER, P. MITTNER e I. Scoront: Nuovo Cimento, 
4, 491 (1956). | 

(*) P. Bassi, R. Cano, S. FOCARDI, A. MICHELINI e F. SAPORETTI: Report on 
« Distortion of tracks in bubble chambers », presented to the 2nd International Con- 
ference « Atoms for Peace » in Geneva (to be published). 
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dalla dettagliata descrizione che segue, ma si ritiene opportuno segnalare 
particolarmente: 


1) la struttura composta delle finestre (Fig. 1); 


2) il sistema di riscaldamento a sorgenti indipendenti situate nel blocco 
metallico; 


3) il meccanismo rapido di espansione, a membrana mobile fra due dia- 
frammi metallici (Fig. 2); 


4) il manometro a ponte, preciso e sensibile, per seguire dal locale di 
comando e controllo l'andamento della pressione nel propano; 


5) lo specchio sferico condensatore usato nel sistema di illuminazione. 


La Fig. 3 in cui la camera appare montata sul suo telaio con dispositivi 
sussidiari, fra cui la macchina fotografica stereoscopica, mostra come sia stata 
realizzata una struttura compatta, assai robusta, facilmente trasportabile e 
mobile. 

Lo strumento è stato costruito e provato nell'Istituto di Fisica dell’ Uni- 
versità di Padova per essere impiegato al sincrociclotrone da 600 MeV del 
CERN. Nella prima esposizione, compiuta di recente, non si è verificato alcun 
inconveniente degno di rilievo. In quasi otto ore di funzionamento continuato 
sono state fatte circa sedicimila fotografie al beam di x+ da 130 MeV; la loro 
qualità è adeguata allo svolgimento del previsto programma di misure e calcoli 
diretto a determinare la sezione d’urto differenziale nello scattering elastico 
del pione su protone. 


2, — Struttura e funzionamento. 


Il volume utile, di dimensioni 7 X 14 X21 em, è ricavato da un blocco di 
ottone massiccio (Fig. 4). Le facce 14x21 cm sono in realtà due grandi finestre 
ognuna delle quali è chiusa con due lastre di vetro temprato spesse 30 mm: 
fra esse è situato un foglio di plexiglas di spessore 4 mm, destinato a pro- 
teggere la lastra esterna dal diretto urto di schegge nel caso, molto impro- 
babile, che quella interna cedesse alla pressione. Il foglio di plexiglas è immerso 
in olio mantenuto a pressione metà della massima raggiunta dal propano, 
cosicchè lo sforzo viene diviso fra le due lastre. 

Questa soluzione ha permesso di rinunciare ad un involucro esterno, dap- 
prima preso in considerazione per raggiungere il grado di sicurezza ritenuto 
necessario. Poichè l'insieme delle lastre ed intercapedine offre grande resi- 
stenza alla propagazione del calore, anche le caratteristiche termiche sono 
soddisfacenti: la continuità ottica è, d’altra parte, assicurata. Il regime di 
temperatura viene raggiunto e facilmente mantenuto con buona stabilità, 
mediante quattro resistenze elettriche indipendenti collocate in apposite cavità 
nel corpo metallico. 

Una delle facce 7x21 presenta tre grandi fori di diametro 45 mm, ester- 
namente ai quali è sistemata la membrana formata (nell’ordine) da un foglio 
di mylar di 0.1 mm di spessore e da uno di gomma di 5 mm. La membrana 
è situata fra due diaframmi di ottone fittamente forati che le consentono una 
corsa massima di circa 2 em: di essi uno è sempre immerso in propano, a 
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distanza tuttavia notevole dalla zona illuminata, e l’altro in aria. Il complesso 
formato dalle lastre forate è facilmente estraibile a mo’ di cassetto, il che 
rende semplice e facile il cambio della membrana. Contrariamente a quanto 
si temeva, il caratteristico ciclo termodinamico dello strumento non è stato 
affatto turbato dalla presenza del diaframma nel propano, che costituisce 
una preziosa barriera rigida di sicurezza. 


Taie, 24, 


Il meccanismo di espansione è tanto rapido che, se non entrassero in gioco 
altre limitazioni, potrebbe compiere fino a 30 cieli al secondo: la grande super- 
ficie della membrana favorisce poi certamente la formazione di onde di pres- 
Slone pressocchè piane nella camera, requisito necessario per avere distorsioni 
piccole. Al di là della membrana, oltre un volume funzionale che costituisce 
il retrocamera, Si trovano, installate su una flangia terminale, due valvole 
elettromagnetiche a tre vie, disposte in parallelo: esse aprono il retrocamera 
o verso due serbatoi ad alta pressione (fase di compressione) o verso altri due 
a bassa pressione (fase di espansione, contro pressione) 


res 
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X21 em opposta al meccanismo di espansione è ricavata 


esternamente una nicchia destinata all’ingresso del beam: una nicchia analoga 
si trova su una delle due facce 7x14 cm, cosicchè vi è possibilità di scelta 


per il percorso delle 
nella camera. 

I fili di una coppia termo- 
elettrica raggiungono il propano 
attraverso un foro nel blocco 
d’ottone. Un altro foro passante 
di diametro 10 mm è chiuso da 
una lastrina di acciaio le cui defor- 
mazioni elastiche, con ottima ap- 
prossimazione proporzionali alla 
differenza di pressione sulle due 
facce, si trasmettono agli anodi 
mobili di due triodi RCA 5734. 
Questi (Fig. 5) sono stati montati 
su due lati adiacenti di un ponte 
di Wheatstone e disposti in modo 
che i segnali di placca dovuti alla 
medesima deformazione della la- 
strina diano origine sulla diagonale 
di uscita ad un segnale composto 
di intensità pressocchè doppia, 
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mentre le derive di caratteristica e di accensione dei triodi si compensano. 
Il dispositivo, la cui frequenza propria è circa 10 000 Hz, è idoneo a funzionare, 
nel nostro caso, per pressioni così statiche come variabili. Equilibrando il ponte 
sul valor massimo della pressione, si hanno 80 V in uscita per il valor mi- 
nimo: il segnale era quindi di intensità sufficiente per essere inviato diretta- 
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mente all’oscillografo nella stanza 
di comando e controllo, situata a 
circa 100 m di distanza. 

Si è rivelato inoltre particolar- 
mente utile controllare spesso, du- 
rante l'esposizione, il funzionamento 
della camera con una macchina foto- 
grafica «Polaroid » che fornisce in 
un minuto positive sviluppate. 

L’inizio della sensibilità del pro- 
pano si manifesta verso i 56 °C con 
la produzione di bollicine isolate al- 
l’atto dell’espansione in presenza di 
radio; il funzionamento diviene sod- 
disfacente fra 58-°C e 60 °C, inter- 
vallo di temperatura cui corrisponde 
la tensione di vapor saturo di circa 


20 atm. In relazione a ciò i serbatoi ad alta pressione vengono mantenuti da 
un compressore a 25 atm, quelli a bassa e l’olio fra le lastre di vetro a 10 atm: 
regolazioni e stabilizzazioni così delle pressioni come delle temperature sono 


affidate a strumenti d’uso corrente 
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Le valvole elettromagnetiche sono azionate scaricando, a mezzo di un 
thyratron, un condensatore attraverso le due bobine in parallelo (induttanza 
complessiva 0.5 H). La scarica cessa dopo un tempo eguale à metà del periodo 
proprio del circuito oscillante; variando la 
capacità del condensatore si regola la du- 
rata dell'impulso di corrente: la tensione 
del condensatore ne determina ovviamente 
l'ampiezza. 

La Fig. 6 è una fotografia dell’oscillo- 
grafo azionato come sopra descritto, in 
condizioni normali di lavoro. La pressione 
del propano comincia a scendere 13 ms 
dopo l’accensione del thyratron: in circa 
5 ms essa raggiunge il valore minimo di 
12 atm e quindi riprende a salire per rag- 
giungere la pressione di vapor saturo. Il 
ciclo può ritenersi concluso a 30 ms 
dall’inizio: il tempo sensibile è di circa 
D ms. 

Nell’esposizione al sincrociclotrone del 
CERN il funzionamento della camera ve- 
niva sincronizzato con quello del sinero- 
ciclotrone in modo che un impulso del 
beam ogni sessantaquattro trovava il pro- 
pano a pressione minima: il flash seguiva 
il beam con un ritardo di (0.6-- 0.3) ms. 
Gli impulsi erano 55 al secondo: le espan- 
sioni si succedevano dunque a circa 1.1 s 
di intervallo. 


3. — Illuminazione e fotografia. 


\ Il sistema ottico è schematicamente 
descritto in Fig. 7. La camera è illuminata 
da un flash PTW, racchiuso in uno schermo 
opaco, bianco diffondente nella parete in- 
terna e dotato, dalla parte opposta alla 
camera, di una finestra rettangolare di 
forma simile alla camera stessa. Uno 
specchio sferico agendo da condensatore 
forma un'immagine reale della finestra, 
di dimensioni eguali a quelle della camera, 
proprio nell’interno di questa. 

Una lente sferica a contatto ottico con 
la camera forma invece un’immagine reale 
dello specchio sul piano degli obiettivi: 
i raggi uscenti dal vertice dello specchio 
| passano quindi per il punto di mezzo del 
LES VE segmento individuato dai due obiettivi. 
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Tali raggi sono assiali per i fasci di luce che (per effetto dello specchio) 
convergono sui punti della camera a bolle: si riconosce che in tal modo si 
ottiene la massima e più uniforme illuminazione media contemporanea degli 
obiettivi da parte di ogni centro diffondente luce dall’interno della camera. 
Per evitare che gli obiettivi vengano anche direttamente illuminati dallo 
specchio basta annerire le due regioni (piccole) ove, sullo specchio, si formano 
le immagini dei due obiettivi date dalla lente. 

I due obiettivi sono Schneider-Angulon f = 65 mm, usati con apertura 1/16. 
I due fotogrammi vengono presi su di un unico film da 35 mm; Vingrandi- 
mento è circa 1/9 e l’angolo stereoscopico 25°. 

Usando film pancromatico di sensibilità 17/10 din si ha luce sufficiente 
con 2J per flash. Uniformità di illuminazione e contrasto possono essere 
giudicati dalla Fig. 8. 


4. — Distorsioni. 


Una valutazione delle distorsioni è stata fatta su tracce di pioni di 230 MeV 
con lo stesso metodo delle differenze seconde usato nelle misure di scattering 
in lastre nucleari. 

Nella seguente tabella sono riportati i risultati delle misure relative a celle 
di lunghezza t eguale rispettivamente a 0.92 em, 1.84 em e 4.6 em. Il A 
rappresenta il disallineamento medio di tre punti, del resto qualsiasi, che si 
succedono a intervallo t sulla traccia e più precisamente, se si segue l’ordine 
di osservazione, la distanza media del terzo punto dalla retta dei primi due; 
ciò evidentemente solo per tracce il cui scattering coulombiano è trascurabile. 


| t ( em ) AS; ASS | ANI 
| 
| 0.92 PRE 20 | 35 ae 3.5 
1.84 75-10 59 45,207 
4.60 190 + 40 245 | = 


Il risultato relativo a celle di 4.6 em è privo di significato: in ogni modo 
è evidente il fatto che il A cresce con la lunghezza della cella molto più 
lentamente che il Au: Si è inoltre trovato che le distribuzioni di Aj, sono, 
entro gli errori statistici, gaussiane con valore medio zero: la camera non 
appare quindi affetta da apprezzabili distorsioni sistematiche. 
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SUMMARY 


A propane bubble chamber of about two liters volume is described: details con - 
cerning the membrane expansion mechanism, the structure of the windows and the 
illuminating system are given. Some features of the use of it, recently made at the 
CERN sincrocyclotron, are indicated. 
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LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Possible Polarization Measurements in K'-H Scattering (*). 


B. J. MALENKA 


Department of Physics, Tufts University, Medford, Mass. 


(ricevuto Pll Gennaio 1959) 


According to present data, the large 
drop in the elastic K°-H differential scat- 
tering cross-section that is apparently 
observed in the center-of-mass angular 
interval — 2 > cos 0 > — 1 is rather dif- 
ficult to explain in terms of a small 
number of phase shifts (1). It is stated 
that the data can be fitted by mainly 
S wave scattering except in this latter 
backward scattering interval which ap- 
pears to be somewhat of an anomaly. 
Hence, additional measurements in the 
region would be most welcome. 

Although there is some question of 
immediate applicability at present inten- 
sities, it is perhaps well to point out 
that further information about the K°-H 
interaction particularly its spin depen- 
dent nature can be obtained from a 
measurement of the possible polarization 
of the recoil proton. 

Assuming a spin zero K*-meson and 
a spin + proton, the scattering amplitude 
takes the comparatively simple form 


A(6) + B(6)o-n , 


where m is a unit vector perpendicular 
to the plane of scattering, o the proton 


(*) Supported in part by the U. S. Atomic 
Energy Commission. 

(2) J. E. LANNUTTI, G. GOLDHABER, S. GOL- 
DHABER, W. W. CHUPP, 8. GIAMBUZZI, C. MARCHI, 
G. QUARENTI and A. WATAGHIN: Phys. Rev., 
109, 2121 (1958). 
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Pauli spin matrix, and A(0) and B(6) 
well known functions of the center-of- 
mass scattering angle 0 and the phase 
shifts. The expression for the diffe- 
rential cross-section and polarization as 
functions of A(#) and (0) are now 
equally well known and conveniently 
obtained (?). 

With the present techniques, a de- 
termination of the proton polarization 
generally involves a measurement of the 
asymmetry resulting from a second elas- 
tie or quasi-elastic scattering of the 
proton at relatively high proton ener- 
gies (*). Presumably, with the advent of 
sufficiently high K* intensities such meas- 
urements could be carried out in emul- 
sions (4?) and bubble chambers as well 


(?) See, for example, K. B. MATHER and 
P. SWANN: Nuclear scattering (Cambridge, 1958), 
p. 269 et seg. 

(3) Experimental references are listed in 
ref. (*) and in L. WOLFENSTEIN: Ann. Rev. of 
Nucl. Sci., 6, 43 (1956). 

(4) J. FRIEDMAN: Phys. Rev., 104, 794 (1956). 

QRBa i HHED and Ba Ga MACIIC-RPAYS: 
vev. Lett., 1, 375 (1958). 

(AC J. Barry and S. J. Gompsack: Proc. 
Phys. Soc., 72, 1130 (1958). 

(7) In emulsion, it is hoped that once the 
usual K*-H scattering event is located and 
measured, it would require relatively little 
additional effort to follow the recoil proton 
track for a nuclear scattering in spite of the 
relatively long mean free path of the proton. 
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as counters (*). However, first we should 
note that, kinematically, obtaining re- 
coil protons of such energies would be 
entirely feasible. 

A simple calculation (8) shows that 
since the rest mass of the K* meson is 
somewhat larger than half the proton 
rest mass, for backward scattering of the 
K*, i.e. forward scattering of the recoil 
proton, better than 90% of the K* ki- 
netic energy is transferred to the recoil 
proton even at low energies. In fact, at 
a recoil scattering angle of 25°, better 
than 75% of the K* energy is transferred 
to the recoil proton. Therefore, at re- 
coil angles less than 25°, for example, a 
200 MeV K* yields a recoil proton whose 


() For formulae containing the non-relativ - 
istic and relativistic kinematics of recoil scat- 
tering, see for example, ref. (*), Appendix A. 


B. J. MALENKA 


kinetic energy is larger than 150 MeV. 
Such energies are ample for the usual 
asymmetry measurements. Finally, we 
note that the center-of-mass scattering 
angle 0 is related to the proton recoil 
angle 0 by the simple relation (non- 
relativistic) 


Thus it is readily seen that kinematically, 
the region of large energy transfer inclu- 
des the region of interest, where 
—2 > cos 0 > —]1 mentioned above (°). 


(*) Similar statements could also be made 
about the recoil proton obtained from K*-neutron 
charge exchange scattering. This case differs 
from K*-H scattering in that it includes the 
effects of both the isotopic spin zero and one 
amplitudes whereas K*-H involves only the pure 
isotopic spin one states. 
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Determination of the Mass of the A° Hyperon. 


J. Bogpanowicz, M. Danysz, A. FILIPROWSKI, E. MARQUIT, 
E. Skrzypczak, A. WROBLEWSKI and J. ZAKRZEWSKI 


Institute of Physics, Warsaw University 
Institute for Nuclear Research - Warsaw 


(ricevuto il 22 Gennaio 1959) 


A stack of 180 pellicles 10 em x 10 em 
of 600 um Ilford G-5 emulsion was ex- 
posed to the Berkeley K~ beam. In the 
80 plates of the stack which were at 
our disposal a search was made for de- 
cays of A° hyperons. The region in 
which the majority of the K~ mesons 
stopped was area scanned for ox- and 
o,- Stars. The pion tracks from the o_- 
stars were traced back to their point of 
origin: o,x- Stars, two-prong stars, scat- 
tered pions from the beam, ete. If no 
unambiguous origin was reached in a 
distance of 20mm from the o,- star, 
the track abandoned. Sixty one 
two-prong stars with a second outgoing 
track were found out of 660 pions that 
were followed in this way. In three 
cases the second track left the stack. 
The remaining 58 with measurable 
ranges were included in the analysis. 


Was 


1 - Calibration 


From an additional scanning per- 
formed on o,- stars, 22 protonic decays 
of the X* hyperons at rest were found. 
Their analysis, using the method devel- 


oped by Fry and WHITE (1), yielded 
for the shrinkage factor of our emulsion 
the value: 


81 = 2.20 + 0.033, 


the value being referred to standard con- 
ditions in our laboratory, and for the 
mean range of the decay protons the 
value: 


ie (V676.2.25 753) um 


Independent information on the shrink - 
age factor was obtained from an anal- 
ysis of 53 two-prong stars classified as 
good examples of the A® hyperon decays 
(see below). 

Taking into account the random er- 
rors of the individual Qa» values and 
also their known dependence on the 
shrinkage factor, we found for the most 
probable value of the shrinkage factor: 


¢) W. F. Fry and G. R. WHITE: Phys. 
Rev., 90, 207 (1953). 
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The weighted mean of the two inde- 
pendent estimates s, and s, 


21 + 0.027 


va) 
I 
bo 


was used in further calculations. 

The density of our emulsion as com- 
municated by Dr. E. J. LOFGREN was 
contained between 3.809 g/em? and 
3.853 g/em*®. Using this information we 
estimated the ratio of ranges in stan- 
dard emulsion ((?), p. 10) to those in 
our emulsion: 


Ry/R = 1.0033 + 0.0045. 


The ratio Ry/R was evaluated inde- 
pendently using our value of the mean 
range of the =* decay protons (*) and 
the corresponding value of this range 
caleulated for standard emulsion from 
the mass value of the X* hyperon from 
reference (3). Calculation yielded: 


Ry/R = 0.9993 + 0.0055 . 
The weighted mean value of these two 
independent estimates of the R/R ratio: 


R/R = 1.0017 + 0.0035 


was used in further calculations. 


2 - Measurements and errors 


As already mentioned, 58 two-prong 
stars tentatively identified as A°>p4 x 
decays were chosen for detailed analysis. 
Ranges of the decay products were 
measured to an accuracy of 1%. Each 
measurement was made independently 
by two observers. A Koristka MS2 
microscope coupled, to a Zeiss Optimeter 
was used for depth measurements. Dur- 


(*) W. H. BARKAS et al.: UCRL-3768 (1957). 

(*) The mean range R=(1679.2+7.3) um 
was obtained by using the weighted value of s. 

(5) W. H. BARKAS and A. H. ROSENFELD: 
UCRL-8030 (1958), VI-41. 


ing each series of measurements, pre- 
viously selected standards of depth in 
the emulsion were measured to allow 
conversion of the data to our standard 
conditions. 

In evaluating random errors of the 
Q-values, errors arising from straggling 
and multiple scattering were compounded 
with errors of measurements of the 
ranges and of projected and dip angles, 
and also with those connected with hu- 
midity fluctuations. 

Systematic errors were split up into 
those resulting from the uncertainty in 
the shrinkage factor, the uncertainty in 
the stopping power of our emulsion, the 
uncertainty in the range-energy relation, 
and the uncertainty in the values of the 
rest mass of protons and pions. 

The values of straggling were taken 
from the work of Barkas and Young (*); 
uncertainty of 0.59% in the Barkas range- 
energy relation (5) was assumed (6); the 
mass values used (3) were 


mp = (938.213 + 0.01) MeV, 


m,—== (139.63 + 0.06) MeV . 


3 - Results 


The histogram in Fig. 1 shows the 
distribution of the Y-values as calculated 
for the 58 two-prong stars. Fifty four 


10 20 30 40 Qin MeV 


Fig. 1. — Histogram of the Q-values for 58 two- 

prong stars with both outgoing particles stopping 

in the emulsion. OQ-values were calculated as- 
suming A" —p+zr decay. 


(*) W. H.  BARKAS and D. M. Youne: 
UCRL-2579 (1954). 
(*) W. H. BARKAS: UCRL-3769 (1957). 


(°) W. H. BARKAS et al.: Padua-Venice Con - 


ference (1957), VI-41. 
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of them are concentrated in the vicinity 
of the mean value; in one of them the 
r -meson is inelastically scattered and its 
energy was determined by ionization 
measurements; therefore the Q-value for 
this event has a rather large error (Q = 
= (39 + 7) MeV). The remaining four 
events are well out of this grouping; of 
these four events, one (Q= 41 MeV) is 
situated 2 mm from the edge of the 
emulsion in a strongly distorted region, 
where measurements are not reliable. 
Three to five events of our sample thus 
appear to be background. Fifty three 


Area corresponding 
to one event 


le) 


40 Q in MeV 


— Constant area histogram for O-values 
of 53 two-prong stars (see text). 


336s 370 380 39 
Fig. 2. 


events yielding closely concentrated Q- 
values were analysed in greater detail. 
A constant area histogram of the 53 Q- 
values where the random errors were 
taken into account is shown in Fig. 2. 


g 47 - Il Nuovo Cimento. 
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The weighted mean value of the Q 
distribution is 


Q = (37.58 + 0.18) MeV, 


in good agreement with the Barkas 
value (5) 


Q = (37.45 + 0-17) MeV. 


The error quoted includes both random 
and systematic errors. The standard 
deviation of the mean of the distribution 
due to random errors is 0.082 MeV, in 
good agreement with the standard de- 
viation of the mean, calculated from 
the spread of the observed @-values 
— 0.084 MeV. 

For the mass of the A° hyperon we 
obtain: 


m = (1115.42 + 0.19) MeV . 


* * * 


We are greatly indebted to Dr. E. 
SEGRÈ for making the stack available, 
to Dr. E. J. LoFGREN for exposure of 
the stack and to Dr. 8. GOLDHABER for 
processing the emulsion. We wish to 
thank Dr. P. ZIELINSKI for many val- 
uable discussions. We also would like 
to express our appreciation to the scan- 
ning staff of our laboratory: Mrs. K. 
BoBINSKA, Mr. R. DABROWSKI, Mrs. M. 
PAZDANOWSKA, Miss W. SANIEWSKA for 
their careful work and especially Mrs. I. 
PRZYPKOWSKA for her efficient help in 
scanning, measurements and calculations. 
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(ricevuto il 23 Gennaio 1959) 


The energy now available in acce- 
lerators is sufficient to produce the & 
particle. Hence the amount of data con- 
cerning the = particle should consider- 
ably increase in the nearest future. It 
is possible that bound systems of & with 
nucleons may be observed. The & par- 
ticle in the vicinity of nucleons enters 
into the fast process H+ N — 2A+30 MeV 
Therefore only those systems can exist 
in which the above reaction is forbidden. 
This forbidness may have two reasons: 
Large binding energy (eg. 5He=, for 
Bs>9.5 MeV, is stable) (1), and the 
conservation of total isotopic spin of the 
system and its projection. A consequence 
of this is that 2p 2°, 2nE E:N and other 
hyperfragments should be stable with 
respect to fast transitions into 2A. If 
we neglect electromagnetic interactions 
the E bound with a nucleon in the 7= 1 
state can decay only in cascade. If one 
takes into account electromagnetic inter- 
actions, then the states 7,= +1 (Ep, 
E n) too can decay only by the cascade 
process E°p>70+A+p=>7°+p+(77+p). 


Il [x] 


(1) P. ZIELINSKI: to be published. 


On the other hand, the state 7= 1, 
= de (1/V2)(E p — Hn), as a re- 
sult of the Coulomb energy between the 
E and p, becomes a superposition of 
the states T—1 and T—0 (the coeffi- 
cient for 7= 0 is proportional to e?) and 
therefore the decay into 2A can take 
place only if there are no additional ef- 
fects which forbid it, but the lifetime 
will be l/e* longer. In this case the 
decay via 2A+y will be a competitive 
process. We have here a ceitain simi- 
larity to the case of the >. 

For the (1/4/2) [= p+ En] hyperfrag- 
ment with 7= 0 there are no selection 
rules resulting from isotopic spin inva- 
riance as for T=1. Here the selection 
rules can arise only from the identity 
of the particles produced. If A has the 
spin 3 then from the condition that the 
function of the two A particles be anti- 
symmetric, it follows that L+S is even, 
where Z is the relative orbital moment 
and S is the total spin. If we take the 
total angular momentum J= 1, then 
we obtain L= S=1 and the system 
has the parity —1. The state P,, = 
= +1 for J=1 does not exist. For 
J = 0 we have two possibilities, L= S=0, 
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Tein a te a Wi Sah ey 
Here there are no restrictions on the 
parity of the system; the same is true 
for J > 1. In Table I we give the parity 
for the EN system with respect to the 
parity of the nucleon for J= 1 under 
the assumption that the E has spin 3. 


Trio Le 
S | II i e SI 
ul | = 
CA | 0 7 = 
OL = | sli 
1 | 1 EAN PRE 
| 1 2 + | = 


If we suppose that the ground state 
is the state with L=0, J=1 (there 
may appear a mixture with L = 2) and 
that Ps = 1, then the rapid decay into 
2 A is forbidden and the hyperfragment 
with 7 =O decaying via 2A+y can 
exist. 

We now present 
which can be drawn as regards the parity 
of the E according to which of the hyper- 
fragments discussed above exist. Let us 
consider several possibilities: 


some conclusions 


1) The charged hyperfragments 
E n, E°p do not exist; at least one de- 
cay of a neutral hyperfragment via 
2A+y is observed. (The non-existence 
of charged hyperfragments is necessary 
to exclude the possibility that the ob- 
served decay into 2A ++ took place from 
the state 4 1. 7o=— 0 which, if Een! 
and E°p did not exist also does not exist 
owing to charge symmetry.) We can 
establish with large probability, that the 
E is bound in the state s= 1 and has 
the parity Pa= 1. 


2) A neutral hyperfragment de- 
caying into 2A+y does not exist. Then 
either the EN forces are too weak for 
a bound state to exist or if the system 
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existed, then there 


lities (*): 


are two  possibi- 


a) the ground state S= 0 witb 
any parity of the &; 


b) ground state S=1, Ps=—1. 


3) There exist charged hyperfrag- 
ments; decay into 2A+y are observed. 
If the frequency of decays into 2A+y 
is greater than the frequency of occu- 
rence of charged hyperfragments, then 
some part of the decay via 2A4y must 
come from the state 7—=0 and we can 
draw conclusions as in 1). In the case of 
the same frequency or of small statistics 
no conclusions can be drawn. 

Experimental confirmation of whether 
any of the above mentioned cases occur 
is difficult, especially in cases 2) and 3) 
which require good statistics for the de- 
cays into neutral particles. On the other 
hand, as regards case 1) it is not excluded 
that with the same accuracy which today 
allows us to state the non-existence of 
the An and Ap hyperfragments we shall 
be able to establish the non-existence 
of the Hn and Hp hyperfragments. 
Then the observation of one decay case 
into 2A+y, or the less probable, but 
easier to discover process or 2A +e+-+e- 
would give us very valuable information 
on the &. 

It is tempting to make an analogy 
between the EN hyperfragment and the 
two nucleon system. The strong inter- 
action schemes suggested at the present 


(*) It seems that even if a fast reaction into 
2A occurs, it is meaningful to speak of a bound 
state. This results from the fact that for the 
reactions in light nuclei or for surface inter- 
actions in heavy nuclei the ratio of the decay 
time to the time necessary to create the hyper- 
fragment is 


Ria (1) |2 0 
(Vra(R))| 0,4 

This expression, for sufficiently large H= is 
distinctly larger than unity and the existence 
of the bound state can be established by inves- 
tigation of the energy of the A-particle. 


time (24) assume identical coupling 
through a pion for the = and N, and 
therefore we can expect that the poten- 
tials Vyx(~) and Vys(m) will be qualita- 
tively similar. The & particle and the 
nucleon are coupled differently with the 
K field, but. also in this case the EN 
hyperfragment such as the deuteron can 
occur only by the exchange of two K 
mesons, which corresponds to a poten- 
tial with a range 1/2Mx. In addition, 
taking into account the fact that gg is 
probably less than g- we see that the 
difference between Vyy(K) and Vys(K) 


@) M. 
(1957). 

(3) B. D’ESPAGNAT, J. PRENTKI and A. SA- 
LAM: Nucl. Phys., 5, 447 (1957). 

(5) W. KROLIKOWSKI: Nucl. Phys., 8, 461 
1957). 


GELL-MANN: Phys. Rev., 106, 1296 


J. SZYMANSKI . 


should not produce any essential chan- 
ges. In such a case the ground state of 
AN® woulde be  unesstate: == 0") S ai, 
The states with T— 1 would not exist 
as bound ones and the establishment of 
the existence or non-existence of EN as 
a bound system would uniquely deter- 
mine the parity of the & particle as 
Pzs=1 or Ps=—1. This analogy can- 
not be applied to the 2nÆ and 2p®° 
systems because of the Pauli exclusion 
principle for nucleons, which forces the 
third nucleon, in 7=4, states to occupy 
a state with J/= 1. 


OK OK 


In conelusion I would like to thank 
Prof. W. KROLIKOWSKI, Prof. J. WERLE, 
and Dr. P. ZieLINsKI for valuable dis- 
cussions. 
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On the Energy Equation for a Gravitating Test Particle. 


C. CATTANEO 


Istituto di Matematica dell Universita - Pisa 


(ricevuto il 31 Gennaio 1959) 


Introduction. 


In a previous paper (1) we have shown that the equations of motion of a test 
particle M in a gravitational field, with respect to a general physical system of 
reference S, can be given a newtonian form 


d'p, 
— = MG, : 
aT 


(1) 
In order to prove this it is sufficient to define suitable standard quantities, relative 
to the system S, and to adopt a particular differential operation d* (covariant trans- 
verse differentiation). 

In addition to (1) another scalar equation holds for M, i.e. (I, (90)). This equation, 
which can be interpreted as the energy equation, will now be examined briefly. 

Let us shortly recall some of the results of I; for all the details the reader is 
referred to the previous paper. 


1. - Summary of the previous results. 


In the space-time manifold V,, let x° (2) be physically admissible co-ordinates 
(xt — var. being time-like lines), ds? = g,,dx°dx* the fundamental quadratic form 
(hyperbolic of the type + + + —), and a a unitary 4-vector tangent to the lines x4 
(vy? = (0, 0, 0, LV ga); Y:= GialV — gu): In every point x of V, let Z, be the flat 
3-space tangent to V, and orthogonal to y. The covariant transverse derivation Ve 
—- introduced in the paper I — operates on a vector field v everywhere orthogonal 


(3) C. CATTANEO: Nuovo Cimento, 10, 318 (1958). This paper will be in the following quoted as I. 
(2) Latin indexes vary from 1 to 4; greek indexes vary from 1 to 3. 
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tensor Viv, still belonging to Z,. Formally such a derivative can be calculated 
as a usual covariant derivative in a riemannian 3-dimensional continuum with a 
sort of Christoffel symbols constructed by means of the tensor 7,6 = 9a8t Va? 
instead of 4,3, and the systematic substitution of the ordinary partial derivation 0, 
by the transverse partial derivation D, — da (Va/Ya)da (see: I, (55), (56), (57)). The 
definition of the absolute transverse differentiation d*v; is then immediate (I, (53), 
(63), (64)). 

For a moving test particle we consider the absolute quantities my (proper mass), 
U= dxi/dr (4-velocity), P=m,U (4-momentum); and the following relative standard 
quantities: a) standard spatial metric tensor: yap= 9x3+Ya?p (8): b) relative standard 
time: dT = — (1/c)y;dai; c) relative standard 3-velocity: v*= da*/dT, v? = Vap vb; 


to y (vy= 0) (or on a tensor field everywhere belonging to X,) and produces a 


d) relative standard mass: m= m,(dT/dt) = mo V1 — v?/c?; e) relative standard 
3-momentum: p%— P* = mv*. All these quantities, as well as the notion of trans- 
verse derivation, are invariant in all internal changes of co-ordinates (I, (4)). 

By means of these definitions and the differential operations recalled above, it 
has been possible to show that if we project orthogonally on £, the absolute 
equation of motion of a freely falling particle, DP/AT = 0 (4), we obtain the three 
equations (1), where the gravitational field G, has the form (I, (101)) 


yl 


G,=G,+@, 


(2) ; Sia Va yi AGO 
Ga 50 WAI a Gg = €2,p0". 
VA y 


Va 


The antisymmetric spatial tensor (5) Q,g= valda(Y 5/74) — dglvalva)] is the projection 
on XY, of the tensor 4-curl Q,;= d;y; — 0,y;. We add here that the last term in (2), 
Gi, is invariant with respect to every internal change of co-ordinates (I, (4)) whereas 
the two former terms are separately invariant only for a spatial change of co-ordi- 
nates (I, (5)). However their sum, Gi, is invariant for a general internal change, 
as we can recognize, for instance, by putting it in the form G,=e2Q,;7°. 


2. - Energy equation. 


Let us now consider the fourth equation, (I, (90)), which completes the absolute 
equation (I, (85)): 


(3) Sage: 3950 = 


= 0 
ar 


If in its left-hand side we separate the terms which have one or two indexes equal 
to 4, and we remember the expression P!— m(dxt/4T) = — me(1/ys+ypv3/ey4) 


(*) It is the same tensor employed in the transverse covariant derivation. This tensor had 
already been considered in general relativity; see, for instance: A. LICHNEROWICZ: Théories relati- 
vistes de la gravitation et de l’électromagnétisme (Paris, 1955) : = ivi 

8, » P. 9; C. MOLLER: The theory of relat 
(Oxford, 1952), p. 238. Re ne 


(*) D is the symbol of absolute differentiation in V,. 
() We have changed the previous notation writing Qxp instead of 74 9,8 as appears in [I, (101)]. 
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(see: I, (72)), we can give the equation (3) the following form 


dP Il vo day 

4 9 454 

(4) -——=m|-daYag:®"v8 + cy, da (È A AE 
d7 2 Va Va 


In order to give to (4) a physical interpretation, we shall add the following 
definition: standard effective relative energy 


(5) H= — ey°P,=e—='mc3. 


As well as the other standard quantities (47, m, v*, p*), E has an invariant 
meaning with respect to any internal change of co-ordinates and has the same form 
as the corresponding quantity in the Special Theory (°). 

If now we substitute to P, its equivalent expression P,= mey, — (1/c)y,E, (4) can 
be written 


dE dx) y nr 
6 —— = — mo? | (E) met 098, 
ci ar | = \Va) ai 2 = Va 


The first two terms in the right-hand side are simply the instantaneous power, 
AL/AT=MmG0%= mG,v*, of the gravitational force acting on the particle. The 
last term may be interpreted as the rate of change of the energy of the particle 
due to the dependence on time of the system of reference. 

If the physical fluid of reference is stationary (g,. independent of 4) the last 
term vanishes and (6) assumes the classical form 


(7) dE — dL: 


the increase of the energy of the particle equals the amount of the work done by the force. 

In the same case of a stationary system of reference one can also introduce 
the definition of the total relative energy of the particle, recalled in (I, (84)) and 
already employed in general relativity (see, for instance, MOLLER: op. cit., p. 294): 


(8) H =—cP,=—me’*y, =— Ey, . 


H is not invariant with respect to a general internal change of co-ordinates (it has 
an invariant character only by respect to a spatial change of co-ordinates [I, (5)]). 
However it can be noted that also the stationary character of the system of refer- 
ence is not conserved in a general internal change. So the physical meaning of H 
is linked to the fluid of reference with a particular system of clocks which makes 
it stationary. 

In such a system (7), or more immediately (4), gives the well known integral 
of motion (see MOLLER: op. cit., p. 294): 


(9) 11 = COMA 


which justifies for H the name of total energy. 


(*) By means of the standard definitions, the gravitational field does not enter explicitely in 
the definition of the dynamical quantities: it has only a direct influence on the space and time 
measurements. Occasionally we remember that the adoption of the standard spatial tensor and 
the standard relative time gives to ds the form [I, (68)1: any optic anisotropy of the space disappears. 
the speed of the light being equal to c in every direction. 
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On the 34 Wave Functions in Atoms of the Iron Group. 


I. FIDONE 


Istituto di Fisica dell’ Università - Palermo 


(ricevuto il 4 Febbraio 1959) 


Among tbe parameters which appear 
in the theory of the hyperfine structure 
and quadrupole interaction of the iron 
transition group (1) there is the average 
value of 7-3 for the 3d shell. Since the 
experimental values of <r-*) for a certain 
number of elements have been estab- 
lished (2) it seems interesting to try to 
calculate theoretically this parameter by 
choosing the appropriate wave functions 
for the elements with atomic number, Z, 
from 21 to 29. A theoretical calculation 
of <r*> for Cu has been performed 
using Hartree self-consistent wave func- 
tions (3). Numerical tables for all 34 wave 
functions are not available and in every 
case require numerical integration. For 
this reason the use of simple analytical 
forms is advisable. 

The most simple analytical form of 
wave functions has been given by SLA- 
TER (4) and analytical expressions less 


(*) A. ABRAGAM and M. H. L. PRYCE: Proc. 
Roy. Soc., 205, 135 (1951). 

() A. ABRAGAM, J. HOROWITZ and M. H. 
L. PRYCE: Proc. Roy. Soc., 230, 169 (1955). 

(5) R. M. STERNHEIMER: Phys. Rev., 405, 
158 (1957). 

(4) J. C. SLATER: Phys. Rev., 36, 57 (1930) 


simple than those of Slater have been 
suggested by LòwDpIin and APPEL (°) 
(Table V of ref. (°)). However, the 
functions of Slater and of Lowdin and 
Appel cannot be used in the calculation 
of <->) since the main contribution 
comes from the part of the wave func- 
tion within 0.5 atomic units (this region 
is smaller for larger Z; for example, 
for Copper it is ~ 0.3 a.u.) of the nucleus 
and in this region these functions are 
inaccurate. What we need, therefore, 
is a wave function accurate within this 
distance of the nucleus. 

In the paper of ref. (5) L6wp1n and 
APPEL suggest simple expressions for 
3d wave functions in a region very close 
to the nucleus. These functions are given 
by the following formula: 


VA 
fsa(Z, 7) = O3(Z) r3 exp (- 3 | > 


with the coefficients C,(Z) given in 
Table VI of ref. (5). These wave func- 
tions are of a conveniently simple form 


(5) P. LÜWDIN and K. APPEL: Phys. Rev. 
403, 1746 (1956). 
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but agree with numerical tables deduced 
by the self-consistent method for V2* (5), 
Mn?+ and Cut (7) only for very small 
values of r. Nevertheless it is always 
possible to use an equally simple form 
with the same coefficients 0, by modi- 
fying properly the exponential, e.g. an 
expression of the type: 


CNE Zs 2) — Csr exp (- 2 ) È 
It would perhaps be more appropriate 
to introduce a Z,,; but for ease of cal- 
culation it is simpler to vary the denom- 
inator. Using the numerical self-con- 
sistent wave function for V?+, Mn?+, Cut 
we attempt to fix # in such a way that (1) 
is in accordance with these wave fune- 
tions. When this is done for the 
three preceding elements we find that 
the best values for x satisfy the fol- 
lowing linear equation: 


(2) x = 7.50125 —0.12575Z . 


We assume that (2) is valid for all the 
elements in the iron transition group. 
The value of 0, for Cu has been taken 
by extrapolation and is approximately 
equal to 251.5. 

The value of «<-> has been cal- 
culated using (1) and (2) in the usual 
manner: 


Cr) = car) r2dr , 


in which R(r) is the normalized radial 
part of the 3d wave function. Since the 
part of rR(r) outside 0.5 does not con- 
tribute to the integral and since [rR(r)], 
near the origin, is well approximated 
by fsa we have 


2 
603 , 
r= — 7. 
73 (22)4 


HN 


() B. H. WorsLEY: Proc. Roy. Soc., 247, 
390 (1958). 

(7) D. R. HARTREE: The Calculation of 
Atomic Structures (New York, 1956). 


The numerical values obtained using (3) 
are in very good agreement with the 
experimental values (?) as is shown in 
Table I. 


TABLE I. — Experimental and calculated 
values of <r=*»> (a.u.). 
Element x ue or? 

| exp. cale. 

| = 
Se | 4.86050 | 1.2 1.20 
Ti | 4.73475 = 1.70 
W 4.60900 | 2.3 2.30 
Cr | 4.48325 — 2.98 
Mn 4.357 50 3.4; 3.74 
Fe | 4.23175 Dr 4.54 
Co | 4.10600 5.2 5.37 
Ni 3.98025 — 6.20 
Cu 3.854 50 7.4, 7.40 


An interesting question is raised by 
the remark that the observed hfs from 
which the experimental values have been 
obtained has been interpreted in terms 
of an isotropic part from unpaired 
s-electrons and an anisotropic part from 
d-electrons. The explanation of the 
isotropic part is not a one-electron 
process and one may ask whether the 
anisotropic part is a one-electron process. 
The agreement between the experimental 
values of <7-%> and those calculated by 
using simple analytical one-electron wave 
functions leads us to presume that one 
can understand the hfs in terms of 
reasonable wave functions, and that, 
for example, one does not have shielding 
and anti-shielding effects such as are 
obtained in nuclear quadrupole spec- 
troscopy. 


* OK OK 


The author would lke to thank 
Prof. STEVENS, Prof. Fumi and Prof. Pa- 
LUMBO for useful discussions. He wishes 
to thank Prof. SANTANGELO for his 
continuing interest during the course of 
this work. 
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A Possible Experimental Test of the Relative Parity 
of Charged and Neutral K's. 


G. MoRPURGO 


Istituto di Fisica dell Universita - Parma 
Istituto Nazionale di Fisica Nucleare - Laboratori di Frascati 


(ricevuto il 6 Febbraio 1959) 


In some recent attempts by GELL- 
MANN (1), SCHWINGER (2) and particu- 
larly by Parts (3) the idea has been 
developed of examining the consequences 
of introducing in the Lagrangian of the 
strong interactions rather high symme- 
tries. In particular Pats (3) has noted 
that a too high degree of symmetry is 
incompatible with the experimental data, 
because it implies some incorrect rela- 
tions between the observed rates of 
some processes and also the vanishing 
of certain experimentally non zero cross 
sections. 

To avoid such unpleasant conse- 
quences of this excess of symmetry one 
has to destroy part of it. PAIS (4) inves- 
tigates whether it is possible to do this 
by assuming opposite parity for the 
charged and neutral K and introducing 
an appropriate KKr interaction (the 
strong interactions are still assumed to 
conserve parity). 


(1) M. GELL-MANN: Phys. Rev., 106, 1296 
(1957). 

(2) J. SCHWINGER: Ann. Phys., 2, 407 (1957). 

(*) A. Pars: Phys, Rev., 110, 574 (1958). 

(*) A. Pats: Phys. Rev., 112, 625 (1958); 
we shall use the notation of references (*) and (4). 


Notice that if the charged and neu- 
tral K have opposite parity the [ K] inter- 
actions destroy the charge symmetry. 

The purpose of the present note is 
simply that of pointing out an experi- 
ment which is practically feasible and 
which should be appropriate to decide: 
a) whether charge symmetry is violated 
in the [K] interactions; b) whether such 
violation, if it exists, has to be attri- 
buted to a parity difference. 

We start with some remarks: ob- 
viously the first possibility to which one 
may think for examining a parity dif- 
ference between the charged and neu- 
tral K consists, as Pais has noted, in 
exploring the reaction K-+ Nucleon + 
— K-+Nucleon, with and without charge 
exchange. A different low energy behav- 
iour of the charge exchange versus the 
non charge exchange cross section should 
be indicative of opposite parities of 
charged and neutral K. However for 
incident K there are in such reactions, 
too many channels open; though events 
both of the kind K-+p-> K +p and 
of the kind K +p > K°+n, have been 
observed, à quantitative comparison of 
the energy dependence of the above 
reactions seems in practice to be dif- 
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ficult. The situation would be much 
better if it might be possible to perform 
the charge exchange and non exchange 
reactions of K* on free neutrons; there 
infact only the two channels K*+n + 
> K'+n and K'+n + K°+p are open; 
however the most free neutrons of which 
we can dispose are to be found in the 
deuteron and there the effect of the 
Pauli principle introduces a spurious 
inhibition of the exchange reaction which 
has nothing to do with our problem. 

The situation is, on the contrary, 
experimentally much more favourable 
and clean if we consider the reactions 
induced by pions. The obvious method 
for investigating the charge symmetry 
is that of bombarding a self conjugate 
nucleus with z+ and x and compare 
the rates of the «mirror » (5) processes 
for equal energy of the x+ or 77; in our 
case the most favourable nuclei are 
probably D or 4He. In this way it 
should be possible, very simply, to give 
an answer to question a). As to ques- 
tion b) the answer can be found simply 
in comparing the angular distributions 
and energy dependence of the various 
«mirror » processes induced by the z+ 
and =~; a difference in parity would 
certainly imply strong differences in these 
angular distributions and energy depend- 
ence near threshold. In particular, con- 
sidering the reactions in 4He the most 
appropriate and probably most frequent 
reactions of associated production in- 
duced by 7 and respectively x+ below 
the Y threshold should be: 


n-+4He—> *H +A°+R°, 
(1) 
n++4He > *He+A°+K*, 


and also: 
n “He *H, --K°, 
n++4He— ‘He, +K*: 


(8) The word «mirror » refers here to the 
simultaneous change of rt in x, p in n, 
Tesero ina 
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Both the reactions (1) and (2) should 
be very clearly interpreted in a ‘He 
bubble chamber. 

We conclude adding a few points: 


1) It is clearly not very « aesthet- 
ical» to gain the high symmetry of the 
doublet approximation, (f) loosing the 
old good charge symmetry; so we cannot 
help from expressing our hope that the 
reactions considered above really behave 
as mirror reactions. 


2) One way to conserve the ad- 
vantages of the doublet approxima- 
tion without renouncing to charge 
symmetry might consist in renouncing 
to parity conservation in some « moder- 
ately » strong interactions; writing the 
doublet approximation as it is written 
for even p(K), but adding, for instance, 
the KKz interaction. 


€ 


3) If it should however happen 
that in the reaction Z+N > A+N 
parity is too strongly (7) violated, one 
would probably have to renounce also 
to the doublet approximation (at least 
without derivative couplings); infact if 
we insist on having charge independence 
and time reversal for the [Xxz] inter- 
action, such interaction must conserve 
parity, and the same should be true, in 
the spirit of the doublet approximation, 
for the [A7] interaction. 


Note added in proof. 


The same experiment has been pro- 
posed by A. Pars in the December 1, 
1958 issue of Phys. Rev. Letters, which 
was available here only after this paper 
had been submitted for publication. 


(5) Compare, in addition to ref. (*), J. TIOM- 
MO; Nuovo Cimento, 6. 69, (1957); N. DALLA- 
PORTA; International Conference on Mesons and 
Recently Discovered Particles, V-3 and Nuovo 
Cimento, 7, 200 (1958). 

(7) This means: if it should be necessary 
to attribute such violation to the strong pion- 
baryon interaction. 
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Note on Moller’s Energy-Momentum Pseudo-Tensor (*). 


B. E. LAURENT 


University of Stockholm - Stockholm 


(ricevuto il 20 Febbraio 1959) 


In a recent paper (1) C. MoLLeR has shown that the energy-momentum pseudo- 
tensor of matter and gravitation introduced by him (?) may be obtained by varying 
the curvature-scalar density in contrast to the Einstein energy-momentum pseudo- 
tensor that is obtained by varying the pseudo-scalar density (3) 


1) «ee Re, 1e)" 


Mollers pseudo-tensor has the advantage over Finstein’s pseudo-tensor that it 
gives an energy-density and energy-current transforming as a scalar- and a vector- 
density respectively when purely spatial transformations are performed. 

In (M2) the variation in question is performed with the curvature-scalar den- 
sity 3 regarded as a function of g‘” and their first and second order derivatives. 
It is shown in this note that a simpler procedure is to regard 


(2) Re gl Pu 
as a function of the variables g!”, 7%, and I, 5 (4). R,, 
position of /%;, and their derivatives. 

The local variations in g#” and 7%;, induced by the coordinate transformation 


being the ordinary com- 


(3) ao! = al 4 En(x), 


are in our case 


(4) OYA Use EP Pure Te Er 


Oy 


(*) This investigation was sponsored by the U.S. Air Force. 

(1) C. MOLLER: The energy-momentum complex in the general theory of relativity. To be published 
in Dan. Mat. Fys. Medd. In the following denoted (M2). 

(?) C. MoLLER: Ann. Phys., 4, 347 (1958). In the following denoted (M1). 

(3) gi” — =. gt”. 

(*) I owe many thanks to prof. O. KLEIN for proposing the use of this method, the Palatini 
method, in the present case. 
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The Y4 can stand for the g“” as well as for the 7%,,. In the former case 


(5) wir = gol dp + ad —g OE, 
(6) ee ey 


while in the latter case 


‘ 


(7) wi uv a 1°, Of hà 17, On CE den, Oy > 
(8) Van = — 2 (05 0p 0) + 050, 07) . 


The variation of N is very easily obtained. Using g/”,,= 0 one finds that (5) 
(9) OU a ogg Rol) CAC OL) gue 


From (9) and (6) it follows that 
{10) Se, ee 


Using this formula, equation (30) in (M2): 


où où 
(11) ¥4 b| ut) =0, 
E È 
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is shown to be true also in our case. 
The fact that 9 is a scalar density yields the following local variation 


(12) ON = — (RE), . 
Using (9) and (4) in the left hand side of this and putting all the coefficients 


of the different derivatives of £* equal to zero gives the formulae corresponding 
to (33)-(39) in (M2): 


(13) S Page Oe 
(14) SP —=V.fr,, 
(15) VA TRE 10, 
(16) Va =i pr sla Gale — 0 ; V Pv = yore ’ 
where 
Om on 
17 SP = —— uf -Y4,—05R, 
(17) spa Wet aya, L'an 
. où 
(18) V2? = —— u#B, 
71, 
on 
vo == m4By 
(19) V bv) = A, vABy . 


(5) See f.i. C. MoLLER: The Theory of Relativity (Oxford, 1952), p. 334. 
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Notice that our Vf have not exactly the same symmetry properties as the 
corresponding quantities in (M2), but still an antisymmetric quantity U,”? fulfilling 


(20) 8,2 = UP, 
(21) Ups Ur = 0, 
can be constructed. It is easy to see from (14)-(16) that 
(22) U,8Y = LV + 4V,,79F 5 — (BY) > 
fulfils (20) and (21). This formula corresponds to (46) in (M2). 
It is now a simple matter to calculate S,f and U,°” explicitly from (17), (22), 


(9) and (5)-(8). 
We obtain 


(23) SARL = 28,8 e CH OSIO) 5 


Using the formulae (8) 


) ox x ! ja Dx 1 x Dx 
(24) da Eee 5 On via CEE Pr 
and 

= 1 CL 8 

25 == - nm, — OP 2, 

( ) V I € Dee ET i + x 


S.? in (23) is easily seen to be identical with ©, in formula (106) of (M1). 
Further we obtain after a short and straightforward calculation 


(26) UE ET ET hee 
By expressing the 7,, in terms of the g,,, and their first derivatives, eq. (26) gives 
(27) Ue “i V— 994° 9?" (Gxc.6 i doro) 4 


This is precisely the expression for 7,8” given in formula (65) of (M1). 


The author wishes to thank professor O. KLEIN and professor C. MoLLER for 
the discussion in which this investigation was suggested. I am also indebted to 
Professor Moller for his hospitality at Nordita, where this note was written. 


(°) See 7.7. TOLMAN: Relativity, Thermodynamics and Cosmology (Oxford, 1934), pp. 222-224. 
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F. G. FRIEDLANDER — Sound Pulses, 
volume in-16° della Cambridge 
University Press, pubblicato il 
27-6-58; 40 scellini. Fa parte della 
raccolta « Cambridge Monographs 
on Mechanics and Applied Mathe- 
matics ». 


Il volume Sound Pulses di FRIED- 
LANDER tratta la propagazione di pertur- 
bazioni non periodiche cd impulsi sonori, 
fondandosi sulla teoria delle equazioni 
differenziali lineari a derivate parziali di 
tipo iperbolico. 

Nel primo capitolo si stabilisce la 
forma generale dell'equazione di propa- 
gazione di piccole perturbazioni in un 
mezzo omogeneo e si dimostra che la 
soluzione di Poisson dell’equazione gode 
della fondamentale proprietà di propaga- 
zione. 

Nel secondo capitolo si estende tale 
proprietà al caso più generale del mezzo 
non omogeneo e si discutono i problemi 
al contorno che sorgono. 

Nel terzo capitolo si estendono i 
risultati precedenti ai fronti d’onda 
d’urto, dove la pressione è discontinua e 
si pongono le basi dell’acustica delle 
onde d’urto, che per la sua analogia con 
l’ottica geometrica viene battezzata 
« acustica geometrica ». 

I capitoli seguenti, quarto, quinto e 
sesto, si occupano di problemi di rifles- 
sione e rifrazione ed infine di problemi di 
diffrazione delle perturbazioni non perio- 
diche, sulla base dei risultati fondamen- 
tali acquisiti nei primi tre capitoli del 
volume. 

Di notevole interesse per il fisico- 
matematico, l’opera sarebbe di proficua 
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lettura per molti specialisti di altre 
branche e precisamente per chi si inte- 
ressa di aereodinamica, di acustica teorica 
o applicata, di sismologia o di prospe- 
zione geologica. 

Purtroppo però il volume riesce di 
lettura abbastanza difficile per i non ini- 
ziati e la trattazione resta troppo spesso 
puramente matematica cioè assai lon- 
tana dalla mentalità dello sperimentale. 


F. S. GAETA 


Dislocations and mechanical properties 
of crystals, edito da J. C. FISHER, 
W. G. JOHNSTON, R. THOMSON e 
T. VREELAND Jr. (J. Wiley, New 
York, 1957), pag. xIV +634. 


Il libro raccoglie i lavori discussi alla 
Conferenza di Lake Placid, New York, 
dedicata alle osservazioni sperimentali 
delle dislocazioni ed alla discussione del- 
l'influenza di queste imperfezioni lineari 
sulle proprietà meccaniche dei solidi. 
I lavori presentati appaiono distribuiti 
in otto sezioni principali, riguardanti le 
osservazioni dirette di dislocazioni, la 
deformazione plastica di monocristalli 
« puri», l’indurimento da lavorazione a 
freddo e la ricottura, le proprietà mec- 
caniche di leghe metalliche e cristalli 
impuri, il comportamento di solidi sotto 
sforzi alternati, argomenti particolari 
della teoria delle dislocazioni, i « baffi di 
gatto » e cristalli sottili, e il danneggia- 
mento da radiazione. Di particolare inte- 
resse è la sezione dedicata alle osserva- 
zioni dirette di dislocazioni, che sono 
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recentemente venute a comprovare la 
vasta massa di lavoro teorico compiuto 
dal 1940 nel campo delle dislocazioni. 
Prove dirette dell’esistenza di disloca- 
zioni nei solidi, e di varie loro proprietà 
quali l'interazione con difetti puntiformi 
e moti di scorrimento sono state otte- 
nute con tecniche varie, quali decorazione 
con particelle colloidali di precipitato, 
studio di figure d’attacco chimico e ter- 
mico, microscopia elettronica, e sono 
descritte per NaCl da AMELINCKX, per Si 
da Dasx, per gli alogenuri d’argento 
da MircHELL, per Al da HirscH e colla- 
boratori, per LiF da GILMAN e JOHNSTON, 
per Ag da Macuiin. Di notevole inte- 
resse sono pure gli studi del comporta- 
mento di solidi sotto sforzi alternati, che 
mirano a dare una spiegazione del feno- 
meno di affaticamento dei metalli, e che 
forniranno probabilmente nei prossimi 
anni valide informazioni sulle interazioni 
fra dislocazioni e difetti puntiformi me- 
diante tecniche di anelasticità. Degne di 
menzione appaiono infine le osservazioni 
sperimentali di fenomeni plastici de- 
scritte per Cu da BLrewIrT, per KCl da 
T. Suzuki, e per leghe metalliche da 
H. SUZUKI e da HONEYCOMBE, la discus- 


LIBRI RICEVUTI E RECENSIONI 


sione delle proprietà plastiche dei me- 
talli a struttura compatta ripresa da 
SEEGER, e le note di LEIBFRIED sulla 
interazione di dislocazioni con vibrazioni 
reticolari. 

L’opera degli editori è stata partico- 
larmente brillante nel corredare ogni 
lavoro di un resoconto dettagliato delle 
discussioni che l’hanno seguito: è così 
possibile rendersi conto dei punti su cui 
esiste ancora disaccordo o che non sono 
ancora ben compresi. L’effettiva vitalità 
ed attualità degli argomenti trattati, su 
cui è riportata l’opinione dei maggiori 
esperti del campo, rendono il libro racco- 
mandabile sia ai ricercatori in fisica 
dello stato solido che ai metallurghi. 
Naturalmente, poichè il libro fissa un 
particolare momento di un campo di 
ricerca in continuo e attivissimo sviluppo, 
alcuni punti di vista in esso espressi 
appaiono già superati a due anni di 
distanza: ricordiamo a questo proposito 
il modello proposto da GILMAN e JOHNS- 
TON per la generazione di anelli di dislo- 
cazione in cristalli temprati, ed il mec- 
canismo discusso da FRIEDEL perla ricot- 
tura di cristalli contenenti posti vacanti 
in soprassaturazione. 

Wh, 125 SMOKE 
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